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MEASURING CARDIAff* DUTPUT 



1. 'the technique of dye dilution' 

The volume? of blocyi a -p£grsoii)s heart pympS per unit 
time (that rs, the rate at which 'it pumps brood) is ' 
called the person's cardiac output. - Normally Irt a per- * 
son at rest this rate is about -5 liters pef minute. But 
aftej strenuous exercise it can rise to more than 30 
liters per minute. ItLcan also be raised or JLower.ed • 
significantly by certain diseases o-f the blooS vessels, 
heart, and n^rvpus system. 

~ In this^unit we shall discuss a technique for 
measuring.. cardiac output known as dye dilution. .The 
technique works as follows.. At time t = 0. a known amount 
D of a dye is injected into a main vein near the heart. 
The dyed blood circulates' through 'the right side of tffe 
.heart, the lungs,* then the left side of the. heart, and 
finally appears, in^the arterial system. The concentration 
^of the 'dye is monittfred at fixed time interyals ^t at \ 
some convenient point ih' the arterial system. Typically, 
At might equal one second... For purposes of the mathe ; 
matical deveiopmentjji 'this ' unit . we shall assume |he 
monitoring i's done in the .aoSrta near the he'art. In 
Exercise 1 it wirll be assumed that "the dye concentration 
is n&nitorecf in a branch artery ^instead, and you will be 
aslced to make appropriate changes in the analysis that 
fallows. . * ** * . 

♦•Normally it will take onl^ a few seconds ftfr the dye, 
td^Tass through* the heart and lungs once and begin to 
appear in ther aorta. 'A typical 'set of readings will be * 
as in Table* L, where -we see one result of* injecting 
D = 5 mg of dye in 9/maiir vein' ne^r the heart at time 1 . 
t = 0 seconds* 1^ we plot these readings .on ..graph gaj>er- 
we. get £?ie points shown'in Figure 1." % 



Our question 15: How may we use the empirical data 
given in Table 1 to 'determine the cardiac output? 

» 

TABLE I 

Typical' Data for the Dye Dilution Technique 



Time (seconds) : 


0 1 


2 


3 


4 


5 


6 


7 


8 


^ Concent ratTbn: 
(mg/1 i ter) 


0 0 


0 


0.1 


0.6 


0,9 


J. 4 


J. 9 


2.7 


Time (seconds) :< 


9 V 10- 


11 


12 


13 


14 


15 


16 


17- 


Concent ra t«ion: 
(mg/1 iter) 


3.0 3.7 


4.0 




4.0 


3.8 


3-7 


a. 9 


2.2 

k . 


Time (seconds) : 

Concentrat t on : 
(mg/1 i ter) 


* 18-, 19 


20 


21 


22 


23 


2k 






1.5 l.l 


0.9 


0.8 


0.9 


d.9 


0.9 


s 





■A ' 



r-4-- \ 



3-- 



2- 



' c 
o 

3 

(0 



c 

. 0) 

o 
c 
•o 



1 



0' 



10 1 15 
Time (seconds) 



20 



25 



Figure 1. Typical readings in the dye dilution technique, 
when D « ^ mg of dye are injected^at time t = 0 seconds'. 



• 2. THE FORMULA «FOR CARDIAC OUTPUT 

2»1 Preliminary Illustration ^ 

* Let, us Set v the stage by considering a somewhat" . 
artificial simplified version of the question. Suppose 
it were possible to set things jup so dlZ of t^ dye 
flowed through the heart exactly once in a time interval 
oi length T seconds, tat a constant concentration of C mg/£ 
{A record of our observations would look* something like 
Figure 2.) Then we" could express the amount of dye by 
the formula D = CV, where V is the volume (dn latefs) of, 
blood flowing through the hearten this time interval, or 
V = D/C. The cardiac output R (the ratt?) would then be 
^ven by**the formula V = RT (volume .= rate * time), whi<;h 
can be written in the fprm * % 

(11 • R tf/T = D/CT> ' • * 

where D, C, and T are"^all known.* Notice t$at CT is the j 
area -of the .rectangle in Figure 2. 




Time (seconds) 



Figure 2. Idealized observations of dye concentration in 
the aorta. / O 4Kr 



tye cannot achieve this situation. Even if we Were to 
wait a very extended period of tiine to achieve a constant 
concentration of dye in ( the bloodstream, this would be 
useless, ^since we would have no w,ay of know.ing how* long 
it \pok all the dye to pass tire monitoring point .exactly 
once „ A # 

How can we modify this simple algebraic commutation 
to analyze \he data of Figure 1, where the dye concentra- 
tion is not constant? 

2.2 Rectangular Approximation * „ 

• There are two essential differences- between the 
idealized observations of Figure 2 and the "realistic * 
observations of Figure 1. One is that in Figure 2 the 
dye concentration is constant. x *The other is less striking 
but *eqaal]V important — in Figured we can identi'fy a time 
intervals during which we know exactly how much dye has 
passed by our monitoring point. „ f 

Let us consider this second difference first. In 
Figure 1. we se*e that the dye concentration rises sharply, 
then falls sharply, and then, jusfwhen we think it is 
going to fall off tp zero, it rises a^ain. This second 
rise occurs at about t = 2oJseconds^ Now, physiologists 
know that 20 seconds is just about long enough for some . 
blood passing through the aorta to make a round trip of 
the body and the lungs, and reappear in'the aorta. ^ 
Apparently what is happening is that most of the dye 
passes through ttfie aorta in the first 14 or 15 seconds. 
, The dye concentration then falls off rapidly (from t = 15 > 
tQ t = 21) 'as ithe rest of the dye trickles through. Then, 
at about t = 21^ a little dye, having completed its round 
trip, appears f o,r *tbe second time and -mingles with what 

♦ 4 W 

is left of the "first-time-through" dye to cause the jump 

* *■ ' 

in the graph. * , 

V * 
We must attempt * to pick out what part of the dye * 

concentration after V = 21 is due to rt f irst-time-thraugfi 11 ^ 



dye*. Notice that jtfst before t = 21 (and especially from 
t ■ 19 to t = 21) the- dye concentration is decreasing at 
a pretty steady^rate. Let us assume that "first-time- 
through" concentration continues to decrease at this rate. 
Then the graph of "first- time- through" concentration, 

Z stead of rising at t « 217 will pass through the points 
, B, and C as shown^in Figure 3. 

/ Jn Figure 3 we simply drew in A, B, and C by eye. 
They are approximately: A £22,0.5), B* (23 , 0 . 3) ,^and 
C (24,0).. They- represent,, at best, a*shrewd gu^QL^pL 
there is no point in agonizing over their ,exact location. 
By the end of this section we shall see that the portion 



21 has only a small effect on our 



of the graph after t 
fina^l result. • . 

. Now we are ready to confront the "first of the two 
essfential differences^ mentioned at the beginning of this 




Time (seconds). 



Figure 3. Rectangular approximate of the dye flow." 



to 




section. In Figure 3 we have, drawn a Succession of 
rectangles^ Each rectangle has base At (the interval 
between observations) and height c(t i ) C*he observed 
"concentration^ at time t^. In our example At = 1*. In 
Figure 3 we have illustrated, for i = 6 » / t i = t 6 = 6 * r 
and e(t i ) = c(6) a 1.4. ' ' / 

Now let us consider the time interval from t. to 
t i+l» of ^eng** 1 At. At the beginning "of that interval 
the dye concentration is observed to be cft^. The volume 
of blood flowing past our observation. point during the 
time- interval risHVAt. Recall that R is ^ rate. If. the 
dye concentration were £on£tant for this time interval, 
the total ^amotrtit of dye flowing 'through in this interval 
would be*c(t i )RAt, of R times the. area of rectangle num- 
ber' i in Figure 3. 

i 

The time interval At is rather small compared to' the 
total time^ involved, and' the, dye concentration nearer * 
changes abruptly, so the error introduced by making this 
approximation is not great. 

Since we have assumed that the monitoring is done in' 
the aorta near the heart, 'all the dye must flow by our^ ^ 
monitoring point between t = 0 and t - Tq. If we add all 
the approximations corresponding to the rectangles from * 
't = 0 to t = T Q we must account approximately for the 
total amount. of dye D: N 



(2) D * I c(t.)RAt = R I c(t.)At, 

i=l y i=l 1 

where n is the number of rectangles. In our example - , 

n = 23 if we count €he 'first two "rectangles , 11 from t - \ 

to t ■ 3, each of which has "height" zero. Thus, 



(3) R 



D. 



\^ c(t r )At 



where the denominator is the* total area 'of the rectangle!; 
in Figure 3. ^ 



ii 



abased* sjbHctly upon the 




^V:^t;t .••c(t)r defined, (1>ut not observed) -for all t between t> - 0 

and t */T A (see Figure 4). This 'function may be approxi- 
^ir'Kv: »*ted by.iitting, a saooth curve.to the observed points 
&f A, B, and C at the 
SS^^-Md*.' Th^ denominator' in (3) is then an* estimate qf the 
%-4:'' : iria under this curve* In fact it is orie of the appfroxi- 



^A^v * Mating sums* used in defining tfye definite .integral' 



l 0 - V n - 

• c(t)dt ="lim-J c(t.)At,,*r 



jVSgi V^ i^Qh we think of a and T Q as given and set\At - T Q /n. 





We can noj& write 
D 



.(40 • R s 7T — ■> . 



c(t)dt 



0 



We must use an approximation sign because our curve c(t) 
is at best a .curve which fits the data well. We have no 
way of knowing if it is exact. 

3. COMPUTATION OF CARDIAC OUTPUT 

3 , 1 Antidif f erentiation * 

How we use Equation (4) depends on the nature of the 
function c(t)-. It may be that a curve can be fitted to 
the data points in Figure 1 which is the graph of a func- 
tion c(t) whose antiderivative C(t) is known. In that 
case we would use/ the" fundamental theorem of calculus to 
compute 4 

c(t)dt « C(T 0 ) - C(0f^ 



0 ' % r 



and then 



-3>2^ Numerical' Methods 

; .More likely, however, there will be no expHcit 
formula forc(t), let alone for its antiderivative. In 
t this case we use one, of. a variety, of ways to estimate 0 
the denominator of Equation^ (4) , and thus obtain an * 
approximation of R. W*e shall list several, and illustrate 
some:iof -them ^ with the data of . Table 1. Hecall that these 
data w/ere\obtained with a dye dosage of D a 5 ftg. t C 

- (a). We can use the denominator of Equation (3) . 
/.This uses the areas of the rectangles A in Figure 3, rather 
than the area ^nder the curve in Figure 4. , 



» In our example, ^ , ^ 

v 2?* 23 * * * 

'c(tOAt = I c(t.)'. = 0 + 0 + 0.1 + .0.6 +.., 
' i»l 1 i«l • 1 



Then 



+* 1..1 +*0.9'\+ 0.8*+ 0.5 + 0\3 = 44.1. 



R * 44^T = 4^pX Oj: 113 liters/second 



\ 

* 6.8 liters/minute. 

(Although the concentration Measurements were taken by the 
second, output is usually measured in liters per minute.) 
Notice that in making this computation ,we replaced the 
last" three experimental points of Table i witfi the points 
A*(22,J3.5), B (23,0.3), and C,(24,0). fhe^ reason for * 
doing thi,s was discussed in JSection 5.2. 

(b) More laboriously, but also more accurately, we 
could ske* tch, Figure 4 on a large sheet of graph paper and 
count the number of squares' that fall between 'c(t) and 
the horizontal axis.^ We *wouTd then multiply ' this total 
by the unit of area i^pre sen-ted «Jby a single Square/ „ 
There .are- also mechanical devices, c^led planijneters^ » 
with which it is possible to trfcce the boundary % of' a * 
region $nd then read an estimate of the area of the region 
from a meter. We^could use one of these instead of count- 
ing squares. > 

* (c) If the interval fa,b] is* divided into n equal 
parts [a = t Q < t 1 < « . .. < t ^ < t ■ b] , then the 
trapezoidal' vule says': . , > ■ ' . ~ 

• .V * ' ■ • 

c(t)dt - L^l ( y o + 2Yi + 2y2 + /. . + ^ , + j . 
a ♦ * ( * ' 

where we have written y 1 = cft^) for,i'= 0,l,2,...,n. 



/ 



In our example, 

24- " 

c(t)dt a 2 A s ; 4 j° (0 + 0 + 0*0.2 + 1.2 + 1.8+... + 2.2* 
>0 , \ ^ ^ 

+ 1:8 + 1.6 + IrO + 0.6 + 0) * 



= |(88.2) = 44.1, 



As in part (a)', R * 6.8 liters/minute. ' * 

v - 

(d) If the interval [a,b] isdivided into n equal 
parts [a = t Q . < t^< ... < t R-1 < t =* b] , where n is an 
even number, then the parabolic rule , 'aJtso known as 
Simpson's ruie^ says: ^ * 



c(t)dt ; * (y Q + 4y x + 2y 2 + 4y 3 + 2y 4 + 



2y 



n-2 + 4 vi + yJ> 



with theu notation of part (c) . 
In our example, , ^ 



24' 
c 



6) 



+ 4(0.9) + .... + 2(0.9) + 4(0.8) + 2(0.5) 
+ 4(0.3) + oj ' 

=4(132.2) 1 



* 44.1 (to the. nearest tenth) 



Again,- R- * 6. 8^ liters/minute. 



10 



15 



] 



EXERCISES* 



H .Assume the dye monitoring takesTplace at a b'ranch artery which 
» » receives only 1/10 of the .blood comln$*from the heart. What 
( -changes are necessary In the analysis contained in* Sections 2.2 
and 2.3? How does this* affect Equations (J) and (4)? v 

*/ - * * ■ . 

' . * ' / 

2, Suppose eft)- fs measured In mi 11 igrams/1 iter,, 1? in seconds, and 

*< 

0 in milligrams. In .what uni ts should # ' „ 



c(t)dt 



be expressed. 



3. Suppose that at time t the dye, concentration is , 
C(t) = -bt(t - T Q ) = -bt 2 + bttjj, 

where *b and Tp'are positive constants. * ^ 

i. Graph c(t)4 , 

b« ,«£lnd R In terms^ of b, T^.vand the total amount D of dye 
Injected. v . , v 



4. Suppose c(t) is as*shown in Figure 5 



SERIC, 




4 u \ * # Time (seconds) 

l 

Figure 5* A hypothetical concentration curve. 

■". . ; .. . re ... ■■ 



\ 



11 



a; Find R in terms of Tq» Cq» and^the total amount^) of dye 



injected. y/ 



b. How is"* affected if (1) T Q is doubled and C Q is kept 
-constant? (2) T Q is* halved and C Q is doubled? 



Find R in terms of Tq, C qi and D '(the total amount of dye 



injected) if c(t) is as' shown i« Figure 6. 




Time (seconds) 



Figure 6. A hypothetical concentration curve. 

» t 

In^an attempt to determine cardiac output, 10 milligrams of dye 
are Jnjected into a main vein near the heaTF! The dye .concen- 
tration Is monitored at the aorta. The following observations 
are made: • K 

" z 5 7 

Time (seconds): 0 1 2 3' 4" 5 6 7 8 

Concentration: a 0*1 0.2 0.6 1.2^2.0 3.0 4.2 5.5 
(mg/liter^ * 9 

Time (seconds): 9 10 IV 12 13 14 15 16 17 

Concentration: v 6. 3 7.0 7.5 7.8 7.9 7.9 7.8 6.9 6.1 
(mg/llter) 

Time % (seconds): 18 19 20 21 22 23 % 24 25 26 

Concentration: 5.4 4.7 4.1 3.5 2.6 2.-1 2.2 2.1 2.2 
• . (mg/liter) 



17 



/ 



a/ Plot these, observations on graph paper. 

b. , At vrf^at time does recirculation begin? 4 

e. What points would yoj^Viiid to the graph corresponding "to A, 
6, and C in Figure 3? 

Calculate the cardiac output R from the data in Exercise 6: 
* 

a* using EqOation (3) directly. 

b. us^ing the trapezoidal rule. 

c. using the parabolic rule. 

/ 

5, ANSWERS TQ EXERCISES 

Throughout Section 2.2, D and'R must be replaced by D/JO and 
R/10, respectively. Equations (3) and (4) are unaffected. 

Milligram-seconds per liter. 




b. (I) halved. (2) unchanged. 



-2* Vsec- J2L */min * [ . J> ****** is ^e area ) 
^ C 0 T 0, Vo I A of the trapezoid - J* 



6. a. 



c — 

g ? 

o ^ 



8-- 

7-- 
6-- 

5-- 

3 — 
2 

I 



.• • • • 



+ 



h- 

10 15 ■ 20 

Time (seconds) 



25 



30 



b. just after 23 seconds. 

c A (2*4,1.2), B (25,0.7), C (26,0) Is one possible answer. 



7. Did you lumber to replace the last three data points irV>the 
table 'bv^ three faints approaching the t-axis? A (24 J. 2) • 



'bvth 

1 B ^25,%^, C"(26,p) will do. 

a. The denominator of Equation '(l) is 
26 26" , 

• rc(t.)At - i C (t.) 

i»1 ' t~1 1 



i»1 

„ Therefore^ 
R 



l«1 



0 + 0.1 + 0.2 + ... + 2,8 + 2x1 + 1 .If + 0.7 + 0 



10 



T06T7 T0ST7 



- 0.094 liters/second 

- 5.^ I iters/minute. 



c( * )dt = tttzt (0 + 0.2 + o.£ + r.2 + ... + 2.8 /n + 0) 

-7 



4) .• 



1T26T 

- ^(213.4) 
=106.7. v 



1.9 



. 14 



' i^Vh * * *%• « 



A^ln.part 7a, ^ 5.6 1 Sprs/minute? 



r26 



c(Odt * (o +.4(0.!) + 2(0,.2)' + 4(0.6f+ ' 

* +.M2.M) + 2(1.4)% 4(0.7,) *"0 

* > » 

-|(320.4) 



- 106.8. 
Again, R = 5.6 liters/minute. 
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. 1. DRUG DOSAGE PROBLEMS 

1«H Gradual Disappearance of a Drug from the. Body , 

The concentration in the b^ood resulting from a 
single^dose o£ a drug normally decreases with time as 
tlie drug is eliminated from the body. (See Figure 1.) 




(hours) 

Figure 1. Jhe concentration 'of a drug In' the blood stream 
^decreases with time. 

t - 

1.2 'What is the Effect of Repeated Doses of a Drug?- 

If doses of a drug were given at regular intervals", 
'what would happen to the concentration of the drug in 
the blood? Would it behave as- shown in* Figures 2 or 3, 
or in some other way? 




. * ^Figure 2. One 'possible effect of successive doses 

* of d drugiv »' * * i s \ 



time 

0 • * 

Figure 3. Another possible effect oi successive doses 
of* a drug . * 

1.3 How to Schedule for a Safe but Effective Drug* Concentration 

For most drugs there is a concentration below which 
the drug is ineffective and a concentration above which 
the drug is dangerous. How can the dose and the time 
between doses be adjusted to maintain a safe but effective 
concentration? „' % v ' 




0 

Figure 4. Safe but effective levels. 
>•* C Q * change in concentration produced by one dosfe 

t A « time between doses 



2.*A MATHEMATICAL MODEL OF DRUG CONCENTRATION 
To give a reasonabTe^nswel^to the two questions 

✓ 

above, we develop formulas from which we can compute drug 



concentration as a function of time. The development 
depends on two assumptions. THe first assumption is 
quite reasonable. The second assumption is reasonable 
in some circumstances but not reasonable in Others, and 
limit s the application of the model we are about to 
describe. * ^ • K 

2.1 The FiEst Assumption 

The first assumption, one that is borne out 
by clinical -evidence, is this: Whatever the mode of 
elimination, the decrease in 'the concentration of the 
drug* in the blood stream will be proportional to the 
concentration itself. If the concentration were doubled, 
the* ra£e of elimination i$ doubled also. If the concen- 
tratidn is reduced by a third, the rate of elimination 
is reduced by a third. The amount being eliminated at 
any given instaat is a fixed fraction of the amount still 
present. 

To model this assumption mathematically, we assume 
that the concentration of dryg in the blood at time t ' 
is a function C(t) whose derivative C (t) is given by* 
^the formula 

Q) * C'(t) = -kC(t) ' 0 

In this formula k is a positive constant, ^called the 
elimination constant of- the, drug. Notice tha,t C'(t) is 
negative, as it should be if it is to describe a 
decreasing concentration. 

i 

2.2 ' Units of Measurement » 

We usually measure the quantities in Equation (1) 
in the following units: * 
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t hours (hr) 

C(t) • milligrams per milliliter of blood (mg/ml) 

*C'(t) 2^!i or mgml^hr- 1 

k hr" 1 ' " 

2.3 Drug Concentration Decay, as a Function of Time 

If we happen to know the concentration of a* drug at 
c a particular time, then we can predict the concentration 
at any later time by integrating both sides of Equation 
(1). Specifically, if C Q 'is the concentration at t=0, 
then we" can calculate C(t) for every t > 0 in the following 
way^ - ^ . 

First rewrite Equation (1) to j»et 

tut) . . 

VTtftf k • 
Then integrate from 0 to t; 

-kt 



in C < 



■*(2) C(t) « C Q e 



Exercise 1, Starting with Equation (l)* carry out In detail the 
steps that lead to Equation (2). [S-i]* . 

To obtain A the concentration at time t>0, we multiply the 

kt -kt 
initial concentration C Q by e . The graph >of Cft) * C Q e 

looks like the one in Figure 5... 



* This reference means that there is additional explanation material 
available in the Special Assistance Supplement a*t the back of the 
unit. , « 
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Figure 5- Exponential model for decay of drug 
concentration with time. 



Exercise 2. Suppose that the elimination constant of drug A is 
k»0.2 hr \ and that of drug B is k=0.1 hr" 1 . Given the same 
ini tiaKcojjfjsntration, which drug will have the lower concentration 
k hours later? 



2J* The Second Assumption 

\ Having -made an assumption about how 'drug concentrations 
decrease with time', we need a companion assumption about 
how they increase' again when drugs are administered. 
What we shall assume is that when a drug is taken, it 
i£ diffused so rapidly throughout the blood that 'the 
graph of the concentration for the* absorption period is, 
for all practical purpose's^ vertical. That is, we assume 
an instantaneous rise in concentration whenever a drug 
is administered. This assumption may not be as reasonable 
for a drug taken by mouth as it is for a drug that is 
injected directly into the blood stream. 

By combining Assumptions 1 and 2, we arrdve at the 
graphs 'in* Figures 2 through 4. 
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3. DRUG ACCUMULATION WIT H REP EATED DOSES 
: — ~ \ — 

3.1^Quantitfes to be Calculated 

v 

What happens to the concentration C(t) if a dose 
capable of raising the concentration /by C Q mg/ml each 
time it is given is administered at fixed time intervals 
of length"t 0 v ? Does the drug accumulate? If so, to what 
level? The next graph shows one possibility, and suggests 
a number of quantities that one should know how^jo 
calculate . • [5-3] 




FTgure 6. One possible effect of 'repeating equal doses. 

3.2 Calculation of Residual Concentration 

- If we let C^ M be the concentration at the beginning 
of the i-th interval and R. the, residual concentration at 
the end of it, we can easily obtain the following table. 
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TABLE i 

CALCULATION OF RESIDUAL CONCENTRATION OF DRUG* 



i c i-i . : R i 



i • C Q multiply > C 0 e" kt 0 

by e" kt 0 J* 



•* add C 

2 t C 0 + C 0 e- kt 0 0 C oS - kt ,0 +V - 2k t 0 

1 \\ + C 0 e- kt 0 + C 0 e- 2kt 0 C 0 e- kt 0 ♦ \ C^O 



> 



^C 0 e- kt O ♦ . ... + C 0 e- nkt 0 



From the table we see that 

4 

(3) R n = C 0 e- kt o ♦ . , . +C 0 e" nkt o 



is the sum of the first n terms of a geometric series. 
The first tei 
Accordingly, 



The first term is C Q e .0 and the common ratio is ef 0 . 



(4) R n = C 0 e- kt 0(i- 




Exercise 5. Calculate Rj and/R 10 for C 0 ■ 1 mg/ml , k » 0.1 hr" 1 
and to ■ 10 nr. (To compare Ri 0 with the result of Exercise k 9 
assume that the data are given to unlimited accuracy.) 
* — ^ ■ — ; ; 

To re'turn^ to Equation 04) , notice that the number 
e -nkt 0 is close t0 o w h en n i s i ar ge. In^ fact, the 
larger n'becomes, .the closer e" nkt ° gets to 0. [5-4] 
As a results the sequence of R *s has a limiting value, 
which we call R: 

a .32 * • 



0 



lim R 
n-H» n 




[5-5] 



If a dose that is capable of raising the 
concentration by C 0 mg/ml is repeated at 
intervals of t 0 hours, then t&e limiting 
value R of the residual concentrations 
is given by the formula 

(5) R = C P . 

- e Kt o - 1 

The number k is the elimination constant^ 
iJf^the drug. £ 



Exercise 4. UsetEq 



t Q given in Exer 



) to find R for the values of C 0 , k, and 
3. How good an estimate of R is R 1Q ? 




3.3 .Results for Lortq intervals Between Doses 

Yrhe only meanii?gfuj~ way to examine what happens to 
the^resi^iual concentration, R, for different intervals, t 0 
between doses is to look at R in comparison with C Q , the 
change in concentration due to each dose. [S-j>] To 
make th^Ls comparison, we form the dimensionless ratio 
R/C 0 by dividing both sides of Equation (5) by C 0 : 



(6T 



C 0 



Kt 0 - Z 

e 0 - 1 



" EquaXion (6) tells us that R/C 0 will, be close to 0 
whenever the time t 0 between doses is long enough to make 
e kt ^ >> 1. As for the intermediate values of R n , we can 
see from fable I that* each R n is obtained from R n -j by 
adding a positive quantity (C 0 e' nkt °)- This means that 
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Ill the R n 's are positive, because R<j is positive. It 
alsomeans that R is larger than each of the R n 's. "In- 
symbols*, , 

(7} 0 < < 



<R 



for all n. 



The implication of this for drug dosage is that 



In 
Kt f 



whenever R is small', the R n 's are even smaller, 
particular, whenever t 0 is long enough to make e K *° » 1, 
the residual concentration from e.ach dose is a^rtipst nil. 
The various administrations of the drug are then 
essentially independent, and .the graph o^C(t) looks 
lilCe the one in Figure 7. < 




figure 7-. Drug concentration Yor lc5ng intervals \ * >. 
1 * between doses . 



3.fr Results for Short Intervals Between Doses- 



between doses is 



rf, however, the length of* time t 0 
so short/that e^J^L * s not very much larger than 1, then 
Equation (6) shows thpt R/C Q . will significantly graater 
than -L. Tne concentrjation will buildup with repeated^ 
do s^es^irtrri^4^*-s-taWJ^ an oj 




R and R # +C Q . fc [5-7] 



cillation between 



Figure Jpx>n page 10. 




St 



10t„ 



t 0 % ^ ^ 

Figure 8. Buildup of drug concentration when interval 

4 \be twee n. doses is short, 9 " * 



4. DETERMINING A DOSE SCHEDULEpFOR SAFE. -BUT EFFECTIVE 
DRUG "CONCENTRATION * \ 

kA Calculating Dose and Interval* . • 

Suppose that a drug^is known to be "ineffective* below 
a concentration C T and harmful above some higher concen- 
tration Cjj. Is -it possible to find, values of C 0 and t 0 
that will*produce a concentration C(t) that is safe (not 
above* Cu) but .still effective (not bVlow C, ) ? To whatever 
extent "the model & valid the answer is -YES, and figure 8 
gives us the clue for how to start. • 

We begin bv«Jooking for values of Cjf.and t 0 that 



make 

(8) , 

Subtraction then yyilds 



R = .C L and C„ + R « . 



'(9) 



'H 



1 ' 



When these values of R and, C„ are substituted in Eqtfation 
CS) , we find that ~ . / . 



(10) 



e kt o 
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* We- then solve, for e kt ° to obtain 

n \ ' • • 

- ' - * C 

(11) • >e kt o- » ^ . ' . 

'«? ... * * 

-When we talce the logarithm of both sides of (11) and 
^ -'divide both, sides of the resulting equation by k, we 
+ learn, that 

IS (12) t tj^ln^ . ' 

H ' \ • 

\ Exercise £ N Solve Equation (10) for e kt 0 to obtain Equation (11). 

; Exercise 6.^$olve Equation^ (1 1) for t Q to obtain^quation (12). 

- f * * u : 

^ 4.2 Reaching an Effective Level Rapidly 

>; " — 7* — — ^ : — — . " 

/ "Jto reach an effective level ^rapidly, administer a 
^V-dose, often called a loading dose, that will immediately 

produce a blood concentration of m$#nl . *This can be s 
\'' foXlrrwed every t Q = £ In ^ hours by a dose that raises 

the Concentration by C 0 = * c h~ c l m £/ ml * • 



^ ^ v 5. EXERCISES 



7 # State two reasons why the model suggested In this unit seems to be 
a' good one. * 

>;> '> "& * Suggest.pther phenomena for which the model described in the^Jext 
v/. Infght be used- 

i: ,9. a) If k » 0.0? hr" 1 , and the highest safe concentration \s e , 
k - * times the ^owes ^.effective concentration-, find the length 
**of time between repeated 'doses that will assure safe^but 
C effective concentrations. > ' 

■? ■ 36 u 

^ ERJC \ \ \ ^ ,7 /\ : . ■ ,\ 



b) f Does (a) give enough information to determine the size of 
each dose? 

HO. Suppose that k « 0.01 hr'Nand t ft * 10 hrs. Find the smallest 
1 

n such* that R >?R . 

n z 

/ / 

11, 



6iven^Cjj^/2 mg/ml , ■ 0.5 mg/ml , and k - 0.02 hr" 1 , suppose 
that concentrations below are not only ineffective but also 
harmful. Determine a scheme /or administering this drug (in 
terms o^f concentration, and times Of dosage.) ^ 

1^. Suppose that k ■ 0.2 hr" 1 and that the smallest effective 
concentration is 0.03 mg/ml. A single dose that produces 
a concentration of 0.1 mg/ml is administered. Approximately 
how many.hours wIL) the drug remain effective? 



6. ANSWERS TO ExkRCISES 

For detailed solutions, see the sections of the 
Special Assistance Supplement referred to in the brackets 
after each answer^ * 

1. See [5-j]. , * , 

. * ■> « • 

2. ■ C A < [5-5] " - * % 



3. R r - 0.36788; R 1Q - 0.58195 

4. R » 0.58l9C(; R and R 1Q agree to four decimal places. [S-10] 

5. See [5-Jj]. 



^-6. See [S-J2]. 



r 7- See [5-JJ]. 

8. See [5-1?]. N ^ 

9. a) t- Q * 20 hours [5-25] 

b) No; but the first dose -could be as large as 2.72 times the 
* minimum effective dose. [5-Jfl] 

10. n - 7 [S-17] 

12 
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11* Give an Initial dose that will bring the concentration to 

2 mg/rnl. Follow this every 69 hours by a dose that will 

raise the concentration by "1.5 mg/ml . [s-20] 

12. About 6 hours f [S-iP] 
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' :V " MODEL "EXAM * . 



^kl\-; : v^^Tboii9h* you Will seldom be asked to take an exam. on * srngle 
#Offc .Mnl'tr' an< exM a clutter of 4 units Is usually made up from a pool 
&. Lof questions similar to those below. ^ - 



1. Assume that^the decay in concentration of a drug 
injected into the blood stream is given by C 3 Cne" , 
/fnd,$$it the drug is .given in such* a way that each dose 
'^U ^te change in the level* of concentra- 

jV;v^:ftioh : of Co* Write, an. expression that gives the residual 
concentration after 3 doses spaced to hours " ap art; i.e., 
find' the concentration at .time 3to. 

State least one deficiency, of the model described 
%n this unit* 

'ir ,v fc: : ; ■ / * 

*jc:.v ;3. Suggest a situation-, different from that described in 
>^J^\ the text, to which this model might be. applied. ^ 

i y x\- 1 ; , - ■+ : : • . W ••: ' 

'j^V**:- - . - 

v ; 1; 4. A certain dose of a drug lis capable of raising the 
>*OV * * J ~ "nuiL » • 

blood concentration of- the drug^byuO.S mg/ml each time 

yK:v:< ^it ( jOT$aKen.'; The decay constant for the drug is 0.1 hr" 1 ; 

dosfes are given every <four~ hours .7 "1 

a^ Find the concentration 'of the* drug Just before* 
the third dose. + % *. 

fe faT'y . Find the concentration just after^the third dose. 

j£ >0im^' ' Given # the drug above and the JcnowleSge that the 

.^Highest safe level of concentration Is 0.9 mg/ml and 
v : ; A the lowpst effective level is 0.6 mg/ml, Revise a 

reasonable schedule (dose size and time interval) for 

administering the drug. " 

iSSbc^.-^.^*^ «'.'.' '• 
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8. ANSWERS TO MODEL EXAM 



1. See 1*able I, page 7 of text. 

2. A drug taken orally, such as aspirin, certainly takes 
a finite time to diffuse into/ the blood stream. Thi^s, 
the assumption of an instantaneous rise in the level 
pf concentration is not realistic for such drugs . 

3. The concentration of active developer in a photographic 
* developing solution might vary in a similar way each 

time replenisher is added to the solution. See 
for other examples. 

4a. 0.SS98 mg/ml 
4b. 1.0598 mg/ml 

5. t 0 * 4.05 hr ; C 0 ■ 0.3 mg/ml . 

- ' . f 

The first dose could be three times this amount. 




9. SPECIAL ASSISTANCE SUPPLEMENT 



Answer to Exercise 1: 



Integration of 



yields 



and, letting C(0) » C Q , 



or 



and finally, 



In C(t) - In C(0) 
In^.-kt 

£iil m e -kt 



-kt 



C(t) 



•kt 



[5-2] 

If the time for the drug to diffuse through the body 
sufficiently to affect the desired organ Is appreciable compared 
to the time between doses, then the assumption of a vertical 
flse In the graph of concentration Is a poor approximation. 
Under these conditions, the graph of concentration versus time *" 
for a srngle ddse might resemble the graph below: 

C 




t ' 

After completing this unit, try to sketch how a series of 
stjch doses might accumulate. If you would like to pursue this 
further, the equation of the graph above Is * 



This equation Is plotted at the top of page SA-2 for two dif- 
ferent values of The diffusion constant k\* The elimination 
constant k 2 lsO;i hr" 1 for both curves. 



SA-1 




0 10 20 30 ' ft 50 hrs. 

Rise and fall of concent rat ion <when diffusion time is significant. 



[5-5] 



Looking at the first two steps of the diagram: 
►C 0 




we see that Ci * Cq +Ri , but Rj «■ C(je 

Therefore, Cj « C 0 +C 0 eT ktd 

• « Looking at the third step- ~ 























n 

















°, *0 2t 0 

we see that C 2 *.C 0 +R 2 > but R 2 « C^e" kt ° 

x . (c 0+ C 0 e" kt 0)e- kt 0 
-C 0 V kt 0 + c 0 e- 2kt 0. 
Therefore, C 2 « C 0 + {c 0 e~ k *° + C 0 e~ 2kt °) m * 
- C 0 +C 0 e" kt0 + C 0 e" 2ktf 0 . 
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♦ • 

Ate reach the results given In Table I (page 7) by continuing 
this process. ^ 



% The term Cge Is the Increase in the residual 
value at the beginning of step n. t 





Note that at. the end of each dose period the residual 
£d~hcentration Is greater than the' last residual amount 
by a smaller and smaller Increment. 



A : 1 


[5-5]' 

Beginning with Equation (k) , 


«n " C.e- kt 0(J 




Use the fact that » ' 
.11. «-«""• -b. R " 

IT*** # 


ll2 R n " <oe kt °(- 


• 

-e" kt o) 


Eliminate parentheses* 


1 -e- kt 0 




Multiply numerator . 
amf denominator, by ft °. 


e Kt 0 . , 





, . There are two pitfalls in looking at a value of R of 

.001 and concluding that it is small. First of all, we do not 
~ know what .001 means physically. It mlgl\t mean .001 kg'/ml , 
which could" be a lethal concentration of many drugs, or It could 
mean .001 mg/ml, which might be an insignificant concentration. 
The number X001 by itself is devoid of physical meaning or 
magnitude. The second pi tfafl Is t that while .001 mg/ml might 

*be an insignificant concentration of one drug, it might be a 
very high ddse^of another drug. * * SA- 
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We can avoid both these pitfalls by not looking at the' 
absolute values of R but only at its size In comparison to Cn 
by taking the ratio of R to C 0 . Thus, if R is .001 g/ml and 
C 0 Is .0002 g/ml, then the. ratio V 

JL » -001 g/ml . - * 
On .0002 g/ml ' * • 

and we see that R Is several times larger than C<). 



As R n becomes larger, the concentration C n after each 
dose becomes larger. The loss jduring the time period 'after each 
dose increases^with larger C n (assumption 1, page 3). Finally, 
the drop rn-xoncent rat ion after each dose becomes imperceptibly 
close tp^the rise in concentration C 0 due to each dose. When 
this^londitlon prevails (the' loss in concentration equalling the 
gain) the concentration will oscillate betweefflfat the end of 
each period and R+C 0 at the start of each period. 



[5-5] * AnBwJr to Exercise 2; . 



-0.8 -0.4 ^ * 

e < e ; "therefore, <C A < Cg 



l' -p -ktn fl-e-" kt 0] 



[S-P] Answer to Exercise 3: 

-nkt 0 ^| 

1 - e 

Co-**! mg/ml ; k »0/1 hr" 1 ; t(j - 10 hr . 
e -kt 0 . e -(0.1 hr-»)M0 hr) m% . t .0.36788 

Rl => C 0 (0. 36788) (1) » 0.36788 mg/ml 

Rio » C 0 (0. 36788) ( 1 1 '. e e '.\°) « c 0 (O.36786)( ]: ;°g°g§ ). 

= C 0 (o. 36788) (-i|22^) =. c 0 (0. 36788) (1.58190) * 0.58195 mg/ml 
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Answer to Exercise 4:' 

* 

C 0 - 1 mg/ml ; k » 0.1 hr' 1 ; t 0 - 10 hr . 

e kt 0 . e (0.1 hr".l)(10 hr) . ftl „ 2 ?1828 ( " f 

R * 2.71628 - 1 " TJvm = (°-58'9«)C 0 - 0-.58198 mg/nl 

\ 

» N 

[S-JJ] t 'Answer to Exercise S: 

Given C. - u" " C W ; solve for e kt °. 

C L 

e kto^ C H- C t + 1 „£tl 
* « r r r 

L C L ■ C L 




Anauer to Exercise 6: 

Given e 1 *' 0 -^, solve for t 0 . 

L * C • 

Take t«%Haa*»rfthm of each side: 1n(e kt °) - 1n|^j 

kt 0 - ln&] 
' L 



Answer to Exercise 7: 

The model appears to be a good one because It Is I 
accord with several common practices of. prescribing .drugs; I 
It accounts for the practice of prescribing an Initial dose 
several times larger than the succeeding periodic doses. 

, .)Th«saodeI aJso provides quantitatively for the pre 
diction of concentration levels under varying conditions- of 
rates In terms of a single easily measured parameter, k. 

What else would you need to know before you could 
actually prescribe a particular' dose rate? 4 . , v 



: C 

[5-24] Answer to Exercise 8: 

Another phenomenon to which the model could be applied 
Is the consumption of alcohol. How often could a can of beer or 
a cocktail be consumed and still not produce a concentration of 
alcohol In the blood at which a person Is legally drunk? 

A very different phenomenon to which thl*s model might 
also be applied Is the burning .of an oH-fashloned woofc stove. 
Here *the rate of burning onheat output Is proportional to the 
charge of wood placed in the stove. There is a maximum safe 
level of burning to be reached as soon as possible, and a lower 
level required to keep the cabin up to minimum comfort. As the 
wood charge Is consumed, the rate of burning, heat output, and 
consumption of wood decrease. 

Sketch possible graphs of heat output versus time through 
several charges of woocL (See Figures 7-8, "Heat Output of a 
Franklin Stove", p^ SS-if Jay W« Shelton, Vfoodburners 1 Encyclo- 
pedia , Vermont Crossroads P>$ss, Waitsfield, VT 05673. 19760 
1 « t 



[5-JS] Answer to Exerbise 9a: 

/• . * 9* 

1 H v * 

given * e and k * 0!050 hr" 1 



t 0 e l c T =1 ln(eV^(20hr)(l) » 20 hr 

1 3 jl± 



[5-J6*] Answer to .Exercise 9b: - % N 

No, not enough InformaVlofl Is given to determine the 
actual size of each dose. We'have only the* ratio of the highest 
safe concentration to th£ lowest effective concentration* lf.» 
the value of one of these 1 1 ml ts were , known, the other could be 
calculated and the difference In concentration »to be produced hy # 
one dose determined. However, the actual dose' required,, to pro- 
duce this change In concentration would depend on the vojujme of 
bloot^In the patient and rtov^qulckly the drug would spread * 
through the entire blood system* * 
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[S-I?] . AnatMr to 10: 



and 



e Kt ° - 1 



j 1 



find n for R > r- R , 
n i 



•The above hnpl.e's c^^S > fc^^-J . 
Some algebra leads to Vf** 0 > (j- 0. / 



Then- 



and 



•" nkt0 <tl4), 



-nkto < 



2 • 



but also "given were k ■ 0*010 hr" 1 ■ * to" ■ 10 for 
so e " nkt ° - e" n( * 01 hr ~ l >< 10 hr > - c "°* 1n - 
7 Therefore/, e"°- ln <{ and e°* 1n >2. 
Taking the 4ogar^hro of each side: 0,1n In e > In 2 
^ 4^ ,* or °* ln > ^ n 2 

' ' % *■ n > 10 In 2 

• " . * n > 6.9. 

Therefore, the smallest nmust be 7». 



4 Anauer to Exercise 11: 
C 4J 



Given " to ■ f 



and, ^ k - 0*020 hr v * £ H ¥ 2l0 jog/mr ; C L - 0.50 mg/ml. 
(50>r)(1>39), - 69hr-. * 




2.0 i*g/i»l 



Cq ■ t*5 mg/ml 



' ' $8 



69 hr' 
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EPIDEMICS 



1 . STATEMENT OF THE PROBLEM 



An epidemic is the spread of an infectious disease 
through a community, affecting a significant flection, of 
the population of the • community . Typically, the ^number 
of infective persons might rise sharply at first, and 
then taper off as the epi4emic runs its course or is 
brought under control. Figures 1 and 2 illustrate this. 

There are two kinds of steps health authorities 
can take to control an epidemic. They can attempt to 
cure .those who are sick, and they can attempt to prevent 
the disease from spreading. Usually^ tfiey will attempt 
both. 

Since the disease is infectious, it seems reasonable 
that reducing contact between those who have or, catry it 
and those who are susceptible to it will help prevent its 
spread. Another means of controlling some epidemics is 
to eradicate the source of infection, £or example, r/Sfent 
populations or mosquito breeding grounds. However, this 
will be of no relevance in^the model we shall consider. 

Reducing contact may be accomplished by reducing 
the number of infective persons in any of several ways 
depending on the nature of the disease and of the commun- 
ity. For example, they may be quarantined, they^may pe 
cured^ ass#ming\jecovery brings immunity and doesnot 
Jeave them as carrWrs, or, in case of their death, their 



bodies may be quickly otemoved . 

Q 

\t what rate will this reduction have to be accom- 
to keep the epidemic under control? Can we pre- 
dict wfj£t jpoftion of the community will eventually catch 
the disease before the epidemic is over? 



plishea 




2. THE .MODEL t 

2„.l Basic Assumptions * ^ 

* . * We shall make the following assumptions about the 
•epidemic we are .modelling: * , 

(a) . The epidemic begins when a small number of in- 
fected persons (perhaps returning from a trip abroad) are 
introduced into a community. <. f 

(b) No one in the community has had the disease 
before, and no one is immune. 

(c) The epidemic is spread only by direct contact 
between a diseased person, or a carrier, and a' susceptible 
person. 

'(d) All persons who have had the disease and re- 
covered are immune. However, some recovered persons may 
be carriers. 

A simplified .description of the progress of the 
eqidemic is shown schematically in Figure 3. ■ In that 
figure we assume that each person is in ex,actly one group 
at a time, and that changes are in the direction' of the 
arrows only* For example, a quarantined person will not 
be released if he is still *a carrier. 

We shall also "assume, for simplicity, that the total 
population of groups, S, I,* and P does not change during 
or shortly after the epidemic. This means, for example, 
that fhere are no births, no deaths^ from other causes, 
and no new people moving into the community. This assump- 
tion is never realized, of course, but it is a reasonable 
approximation to the truth if the epidemic is short. 




"Susceptibles" 



Person^ who 
have never had 
disease anE are 
not immune 



size at tine t 
$(t) 



"Infectives 11 



Diseased 
persons still 
at large 



Recdvered 
persons who are 
carriers 



■size-at time t 
Kt) 



"Post-infective^" 



Recovered 
persons now 
immune 0 



Quarantined 
persons 



Removed bodies 



size at time t 
P(t) 



Figure 3. Progress of an epidemic, 



2.2 Definition o£ the Variables 

Le't us call t 35 0;the time at which* the ^epidemic' 
begins, and let N = the total population. Let S(t), I<t) , 
and P(t) be the number of persons in groups S, jrf^and P . 
respectively at any time t . * Depending on the nature of 
t?he epidemic, t might be measured in hours, days, weeks, ^ 
or r even months. Our basic assumptions tell us amon^g 
other things, that S(0) = N (the total population), that 
P(0) = 0,' and that during and shortly after the epidemic^ — 

(1) S(t) ♦ I(t) + P(t) =? N . 

t; > « * * „ • 

The' number "who have caught the disease by time t is- 
I(t),+ PftJ,* or N iuS(t).. „ . ' . _ • 



2,3 The Spread of the Disease " ; j> 

** ' Each time a person catches the disease - S (O decreases 
by one and I(t) increases by ojie>. ' How frequently this 
happens is' determine^ by how frequently a person, in group 
S- comes* in contact with* 'one £yi group '\\ 

What is a reasonable formula for the frequency of 
thes N e contacts? We would expect it to vary directly with 
S(t) and. also* with I(t). For^ejatirpie , we would expect 
that tripling *the dumber of infectives whiLe holding the 
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number of susceptibles , fixed would triple the contact 
frequency, Sinfivtarly, we would expect that tripling the* 
number of susceptibles while holding fixed the number of 
infectives would also triple the contact frequency. The 
^simplest formula which varies directly with S(t) and with 
I(t) is kS(t)I(t), where k is a positive constant. 

We shall assume that a fixed 'fraction of these con- 
tacts results in "the disease beirtg transmitted from the 
infective to the susceptible* Then the frequency/ with 
which S(t) decreases by one is k,S(t)I(t) for some new^ 
constant k, (0 < k, < k) . In other words, the rate at 
-w hi ch S(t) is- changing -rrs - k, S ( t j i^y-— — § ^ - 

2.4 A Smooth Approximation 

The l,a"s-t sentence of Section 2.3 seems to be a state 
-ment about the derivative ("rate of change") of S(t). 
Strictly speaking, S(t) cannot haVe a derivative, since 
its graph is ndt smooth. It must be a step function* 
(Figure 4)-, with each step being of height one. But it 
is easy to draw a smooth curve, as shown, which is an 
excellent approximation to S(t) ^ It" will never differ 
from the true value /1>y more th'a'n one, 1 which is assumed to 
be a -tiny error compared to %he total p\>pul%tit)n This 
smooth curve has a derivative, and for it we have 

(2) S»(Jt) - -k,S(t)r(t) 

from so me co n>»ant k, > 0.. 

p 

2.5 Reffloyal of Infectives 

It^seems reasonable that the rate a"t which victims/ 
tfie from the disease, and thus enter group P, is propor- 
tional to the number of infectives at any given time. We 
shall extend this to an assumption that the rate of trans 
fer from group I to group P,- for any reason, is propor- 
tional to the s*ize of group I. That is, after 
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S "susceptibles" v 




Smooth approximation 



j t 



, Figure 4. Approximation of 3(t) by a smooth curve, 

"smoothing" as before, 
(3) P*(t) ».k 2 I(t)* 

for some constant k 2 >>0. „ * 



Exercise: 



1. Criticize' this models For example, ate the* assumptions realis- 
tic? Are £hey reasonably translated into mathematical terms? What, 
if any, important^aspects* of the situation are not represented? 



3. CONTROLLING THE EPIDEMIC 

' * » 

, 3.1 Dei inxtion of . "Control^ - * — 

Recall that one question we asked was at*what'rate 
must persons be transferred from group I " to group P to 
keep the epidemic under control. 



» 
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So far we have not §aid precisely what we mean by 
"under control. " Let us recall how the epidemic begins. 
The disease is introduced into t>he. community by a small . 
numbdr of people. So 1(0) is small, P(0) = 0, and 
S(0) * N. 

i 

(A fc word about the symbol * : When We say one expres- 
sion is a good approximation to another, we almost always 
are thinking of the percentage error, rather than the 
actuaU»8ize of the error. For example, -we might welj. 
write 1001 » 1000, but would be very unlikely to write 
2*1, even though 1001 - 1000 ■ 2 - 1 « a l.) 

The more rapidly I(t) grows, the worH^ the- epidemic 
becomes A Let us adopt as our definition of M under con- 
tY<^"th3 / t I(t) s-tops growing (i.e., I 1 (t) '< 0) after 
some 'time. * 



3.2 The Threshold Remoyal Rate 

Can control be achieved in our model? We shall pre- 
sent some calculations, and leave it to you to finish them. 

Dividing Equation 2 by Equation 3: , 



W s -b S(t) 



% 

In S(t) = P(t) * c. 



Putting t ■ 0 and recalling that P(0) = 0 yje get 



and so* 



In S(0) = c 



In S(t) = In S(0) - Ji. P-(t), 
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v Writing S 0 3 S(0) and solvirig for S(t): 

(4) S(t) « S 0 e" k i P(t)/k 2 . 

Mow it's your turn. 

Exercise : 

2. (a) Use Equations 1, 2, and 3 to show that 
I»(t^- (^S(t) - k 2 )I(t). 

(b) Show that S(t) is a decreasing function for all t. 

(c) Using (b), show that, if t > 0-and k 2 2. k,S , thjn 
k,S(t) < k 2 . 

* * 

(d) .Using (a) and^Jc), show that, if t > 0 and k 2 ~> kjSj, 

then I 1 (t) < 0. 



Recall that k 2 is the proportionality constant which 
tells us how fast persons are removed from group I to group 
P (the one we can influence by quarantine, etc.), and k, 
is the one which tells t us how fast the epidemic is spread- 
ing. Exercise 2 shows that we can keep the epi'demic under 
Control if we can establish k 2 > SqICj. This critical" value 
S 0 kj is called the threshold removal rate* It varies* 
directly with k, and with S 0 . ^ut S 0 « N. So we have the 
not very surprising result that the threshold -removal rate 
varies directly with the^ate at which the epidemic spreads 
and with the population.. 

Exer cises \~r* r - ' 

• } * 

3. A yet simpler (and less realistic) mo^el of an epidemic would be 

one without any provision for removal.- An infective remaiYib an' 

infective. * If N = S(t) + I(t), and if we make the same assum^ 

tions as before concerning contact between infective and suscep- 

tibles, we get: 

S'(t) = -k!S(t)(N - S(t)). 
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\ . Writing 3(0) = S 0 , find an expression for S(t). 

*• *(HinVc^witidififerentiate 

:^ t . i.' * 

Q0>y using theUdentity 

: 

1, 1 - v , 
K _ u v - u u(v - u) #; * 

4. Forsthe S(t) obtained in Exercise 3, evaluate lim fr(t) . 
* * What does this imply about the size of the epidemic? 

5. For the Stt) obtained in Exercise'3,v fip^he time t when the 
*jatfe of the spread <of the epidemic is at its maximum. 

i * * 



A. 



t. , ' 4. A MILD EPIDEMIC 

4,1 Extent of the Epidemic 

Now let us ask the question: suppose, k 2 is almost 
but not quite equal to Sokj, so we do aft quite "control" 
the epidemic. For instance, suppose 0.95 Sokj, < k 2 « < S 0 kft 
What portion of tjie community wil,l eventually catch the 
disease? For t > 0, we have remarked that the dumber who 0 
have caughVthe disease by time t is I(t) + P(t). So' if 
the -epidemic lasts for timeXrXi . e . , S'(t). = 0 for t > T) , 
the number we are looking for ik I(t) + .P*(t/. Let us call 
this number the extent,* and writie it E.^^te/ 

4.2*. An Eqftation for the Extent ^ / . * 



To find E,^we shall begin fty observing that P(t) is 
defined foV c aZ^t > 0, not ju§t for 0 < t < T. Figure 5 
'shows the^ graph of a typical: step, function P(t). starting 
* at Ifime t*=*~-p and extending well beyond t » T. .It makes 
it clear that 1>y .some time T**, later than T but not too 
much later, the slope of tfye smooth approximation must be 
close to zero. That is ♦ m > 

(5) ( P'(T*) # * 0. • 
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Figure 5.' Smooth approaimatiofNt© P(t) 
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Equation 3 immediately tells us I(T*) * 0 . But the ewm 
I(tj^+ P{t) does, not change after t = T*, and so 

*(6)* *E = I(t) + *P{t>) *= I(T*) + ^(T*X * P(T*). t 

/ We have assumed the total population does~ not change 
during or shortly after the epidemic. Specifically, let 
/us take this to mean during the time interval. 0 < t ± T*. 
Then, using Equations 3, 1, and 4 in that order, 

P'(t) * kiKt) = ki(N - P(t> - S(t))« 

y(N - P(t) - s 0 e; k i pft)/ ^)^ . 

throughout this interval'. Setting t T* and using- (53 
and (6) 

(7) U 0 * k 2 (N"- 5 - Soe" klE/ki ). 



4,3 An Approximation for e T ^ l£ ^ 2 * * ' * 

: 

The appearrance of both a linear and an exponential' 
term in (7) makes it very difficulty if not impossible^ 
vto solVe- for E. There is a way to circumvent this diffi- 
culty, provided k l E/k 2 is small. Recall that>for'any 
positive x and any'positive integer n " 

e nX = 1 - * + £ - • + r.i> x n * 



with an error of at most 



x n+1 



Setting x .= k(-E/k 2 



and n ■ 2, we obtain the approximation » ' ' - 

f k E 

with an eYj)or of at mogt ^' 1 



lfM)3 
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4.4 EstMnatihft the Extent^ * i 

before we can use (8) we must of course;assUre our- 
selves that this error term is small enough for *our pur- 
poses,.. Recall that S 0 ~ N; that' is, initially virtually 
everyone is susceptible. If we make th# very modest ' 
assumption that "virtually everyone" .means "pver 99^\Ci« 
other, words, the per.sdns who initially ' introduce* the "dis- 
ease constitute less than one percent of the population), 
then we can s^iow. that E < W>. With this restriction on"' 
E the maximum %rror in using ,(8) to estimate e" k V E / k 2 
works oufto be'less than one-half of one percent bf the 
true value. 



c rt takes a lot of messy algebra. to prove these asser- 
tions, and right now that would distract us- from the main 
argument. So we shall leave that algebra for the appendix, 
and proceed with our estimation. 

^ Replacing e " k > E/k * in (7)**by the' estimate given in 
<8.), and also dividing (*7)*by k 2 , gives us 
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0 « N — B - S. 



TT 7 TT" 



^2 - k 2 

Since, again, So ■ N, we can also replace So by N, obtain- 
ing 



0 =' N* E - N 

T7 * 1 * 2k: 



k,B iJc, 2 B 2 

R 2 K 2 



Hki« B 



2k 2 



E = Nk* 



2k. 



Nk r 



kjN 



- 1 



k 2 1 
N - — 



since k a = s 0 k r E Nk i • 



Exercises: 



6. (a) Assume kj = 10~ 6 , k 2 = .95,\and Find the approximate 

* value for* the extent E of tn^ epidemic. 

f (p) Do the same for k. = »99. 




• 4,5 The Relative Rem oval 

Sometimes k 2 /kj i^called the relative removal rate. 
Its threshold value is S 0 , which approximately equals N. 
With this* terminology, (10) says that in a mild epidemic, 
that is, one for which the removal rate is very near its 
threshold', the total number of persons infected sooner or 
later is approximately 26, where d is the amount by which 
the "relative ^removal rate falls short of its ^threshold 
'(« * N - k 2 /k a )/ , . * 

12 
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5/ APPENDJX 

In thi£ appendix we shall* justify the assertion 
made in the first paragraph of Section 4.4. Specifi- 
cally, if . ' m j 

(11) 0.95S o k, < k 2 < S 9 \ 

(the epidemic is nearly but not quite "controlled 1 ') and if 

(12) 0.99N < S fl <*>!*. ' 

*(over^99l of the population is initi^lly'sirsceptible) , 



then - E-<-Mh --■ - 

Ke can make a -very rough estimate of E graphically. 
Writing . 

r 

f(x) = N - x - S 0 e" k i x/k 2 

^we see that E is the positive root of f(x) = 0; that is, 
J^the x-coordinate of the point where the graph of f crosses 
the positive x-axis. 

To get a rough idea what, this graph looks like, we 
first compute , 

f (0) = N - S Q > 0 

(note f(0) is small since S Q * N) and <• 

f(N) . -S 0 e" k i N/k 2 <;0 - ■ . 

(since the exponential "function is always positive), 
thus showing that the graph ^crosses the x-axis between 
0 and We leave it to you (see Exercise 7) to show 
that f u (x\ < 0 for all x, and that therefore the graph is 
concave downward and connot cross the 'positive x-axis 
„ more, than once. 




Exercise: 



7. If the function f is defined by 



; x f (x) = N - x - S 0 e 



-k 1 x/k 2 



for all real x, show that f"(x) < 0 for all x. 
Carabine this with the fact f(0) > Ojjto show the Igraph 
of f crosses the positive x-axis at raost once. 

- x_ . 

V 1 

It follows that if we q an find any positive number, 

M for which f(M) < 0, then we" can concknie JO < f E < M° (see 
Figure 6). We shall now find such a M. v | 

Equation 1-1 can be rewritten * i 



0/95 < 



5^7 



< 1 



or, taking reciprocals and reversing the inequalities,* 



1 < S ° kl < 1 
1 < ~TT OS 




Figure 6. Graph of f(x) p N - x 



kix/k 2 



(See 5-2 for additional /information.) 



v 14 



From (12) we get 

N < irte 



i 



and hence, multiplying by k 4 '/k 2 and using (.13), 



(14) 



0.99 



< (0.99)(0.95) 



With this inequality and, a calculator, let^tw calculate 



f.(}> 



N , 
T 



S A e 



*(k x H/k 2 ) 




N 



< J N - 0.99Ne"* (k i N/k 2V 



N(0.75 - o.99e^ 1/ (°-")( 0 - 95 ^ 
* N(0.75 - 0.7<y -■ -0.01N 0 / 



Thus j is an example of a point M such that f(M) 0. 
Therefore, 0 < E < j . 

^Remember (Section 4.3) that the error in our es*ti- 

*■ • k E/k * 1 f k E) 3 

mate of e» 1 ' 2 is less than • As we remarked in 

* 2 
Section 4.4, with E < W this works out to less than 

one-Aalf of orte percent of the true value. We'll leave , 

the computational details to you (see Exercise 8). 



Exercise : 

8. (a) Show that if 0 < E < N/4, an4 if kj and k 2 are ^restricted 
as in the text, then the error in using (8) to estimate 
e -k l E/k 2 . s less than 0 0032 

(b) Show that, under the same cpnditions, the value- of 

-k E/k ) 
e 12 15 greater than* 0.76, so that the error of 

t 

•part (a) is less than one-half of one percent. 
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6- * ANSWERS T6 EXERCISES * . 

\ The model is (as pointed out in Section 2.1) only ifeasonable 
when short epidemics are analyzed. ..The raode'l does not consider 
interference from other possible outbreaks or epidemics that 
might occur during our time interval. 

There is also no accommodation ' for the spread of disease by 
' infected objects' '(towels, water, infected air, bathrooms, etc.). 
2. a) ' Given ^ 

(2) >(t) = -k,S(t)I(t) 
and t (3r P'(t) = kji(t) h < 

it follows that . * 

• • 

, I(t) * N - S(t) - P(t)- 

P(t) =» -S'(t) - P'(t) 

= kjSftjlft) - H 2 I(t) 

C / = 4 (k,S(t) - k 2 )I(t). 

b) ^ Given t >^0, we want to show that S'(t) < 0. This is 

* trivialrbecause-t > OVimplies^ I(t) > 0; the disease is in- 
* fectious. * • ^ t 

c) Given: t > 'Q and k 2 > k^Sj (or we could write k 2 > kjS(O)).' 

At this point the epidemic has not started, but the instant 

* • * 

, t S(t) decreases implying that S(t) < S 0 . Substituting 

S(t) jfor ST'Hre can write k 2 > kjS(t). * 

d) Given: I» (t) = (kjS(t) - k 2 )I*(t) , t > 0' 
V an4 ^ k 2 > k^S, (which implies k r S(t) A 2 ) 

Vt folldW^hat k,S(t) - k 2 < 0 ' 
|t and i( t ) > 0 

hence 'l'(t) * 0. 



4\/0. / This impli^that eventually everyone becomes ill. 



NET 



0 ifS 0 >|N.\ 



6. aX 100,000 
b) 20,000 



7. 



f"(x> = 



k,E 



k 2 l "° 



8. a) In (8) the error 4s .at most - 

6 



We have seen that 



1 



k 2 (0.99)(0.95) 



(Equation U) . 



Let E = j . Then the J?£ot f 



s at most 



1 

6 , 



1l £ 

k 2 * 4 



a I • 1_ 
" 6 * 64 

= .00277 



k^ 



6 64 



1 



(0.99)(0.95) 



b) If 



0 < E < £ 
4 



. 0 < x — < .. 1 » I JL 



then 



*2 ~ 4 k 2 
From Equation (14) we then have 



Therefore, 



k 2 < 4 ' < 4 ',(0.99) (0.95) 



k*E 



k 2 . 4 ,(0,997(0. 95) " : u -f 



- -0.26581. 
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(rounded down, not up 4 , to be sure" this inequality is 
1 * preserved), and , " 



k E 

e' > e ^°- 2 f 581 * 0.76658 



(rounding down* again).. 
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SPECIAL .ASSISTANCE SUPPLEMENT 



IS-I) 

' When a, B > 0, the functions 
(14) f(x) = -ae" Bx 

have graphs like those shown below in Figure 7. 



0 . ,1 2 3 4 




f(x) Fi * ure 8 - • 4 

♦ 9 
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Subtracting x from -e" x pulls the graph of «e" x away from 
the positive x-axis; Figure 9 shows tiie graph of f (x) = -x - c~ x 
as a smooth curve through the points (0,1), (1,-1 -i), 



(2,-2-i),... 



* 4 



0 * 



* -4 



\ ■ 



h Figure 9. 

; v 

The graph "of : * 

(15)^ f(x)** -x - ae" 6x # a.B > 0 
( behaves'^bas'ically the same way. 3 

The addition oF a constant N > 0 to the formula for f (x) 
* tra^lates^ the curvd vertically upward to produce a curve like 
the graph of . s 

-k 1 w/k 2 



f(w) = J? - w - Xj« 



in Figure 6. 
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1. RADIOACTIVE TRACER TECHNIQUE 

In the human bloodstreanupotaSsium ions (K + ) are 
constantly\moving into and out of the red blood cells"*' 
(erythrocytes); that is, the- surfaces of the erythro- 
y>. cytes are permeable to K ions. Ions move from the 

'plasma into the red cells at ^a certain rate, while 
other ions jfithin the "tells move out into the plasma at 
a certain rate. The determination of these two rates 
.(that is, of the permeability of the cells surfaces to 
K ions in b^oth directions is of great help to -both 
physiologists and doctors in their efforts to understand 
the structure and behavjior of these cells, and thus 
ultimately to combat blcJod* diseases . 

A technique to determine .these permeabilities works 
^ as follows. A fixed' quantity S of radioadtive K l+2+ ions* 
is introduced ii\to the blood. Initially, all" the%e ions 
are in the plasma. The amount remaining in the plasma at 
various subsequent times is determined by taking blood 
plasma samples and measuring the radioactivity present. 
Our ♦problem is to determine the permeabilities from these 
raw d^ta. 

' 2. A CLOSED TWO COMPARTMENT MODEL 
2.1 N otation- . 

. — : — \ * 

We shall establish a mathematical model depicting 
those aspects oA the situation which interest us. Since 
. ' A > \e have no need tV distinguish one red cell from another, 
we shall represent the bloodstream/schematically by two 
boxes, one {fir cells,. the other for plasma (Figure 1). 
If t is the elapsed time since the introduction of the 
K* 2+ ions, we shall denote the amount of K> 2+ ions in the 
-boxes by C(t) and P(t) respectively* Thus, C(0) « 0 
' • and P(0) '—"S. 



(I) 



(2) 



\ 




plasm/ 
P(t) 


' corpuscles 
C(t) . 

* 

1 


■ ->" 




T 2 P(t) 





Figure 1. Two Compartment Model of Bloodstream 



\ 



2.2 Assumptions 

We shall Assume that this two compartment system is 
closed; tha't .is, there is no loss of K** 2+ from the system. 
In our notation this says * 

) • * 

(1) ■ C(t) + P(t) = S. « 

We shall also assume that at time t the number of ions 
moving from* Box 1 to Box 2 (upper arrow) per unit time 
(the trans fer rate) is proportional to C(t), while the 
transfer rate from 'Box 2^ to Box 1 is proportional to 
P(t). If the respective constants of proportionality 
\the coefficients of transfer) are positive numbers 
and k 2 , our assumption says 

(2) " P'(t) -k 2 P(t) +. k!.C(t). 

The units of Y\ and k 2 are reciprocals of time* (for « 
example, min" 1 or hr" 1 ). In Equations (1) and (2) , S 
is predetermined by the experimenter and P(t) is observed 
empirically, so C(t) can bt easily computed. We must 
figure out how to determine k x and k 2 % " 
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/ % 3. THE* FORMULA FOR P(t) 



V ilv - We begin -- by f iT ) din S P W explicitly. It is easy 

'£} to solve for C(t) in Equation 1 and, substitute the resujt 

!s\ „* n ?° Equation 2: / 

;f { _J P f (0 88 jk 2 P(t) + MS -*P(t» 

5/ 88 *|S - (ki ♦ k 2 ) P(t).' 

For the moment let! us assume -that P'(t) is never zero. 
-We can then aivffle Equation (3) by it* right side: 

K ' M - nfr^fel P(t) " *• 



Since P(0) ■ S, the denominator if negative *(= -k 2 S) 
when t * 0. We can conclude that it- is always negative-" 
it is never zero jand therefore cannot change sign. Let 
us multiply this last equation by -1 and then * 
'anti.different.iate both sides: 



V 'i\l) * k,s dt = - f ldt ■ -t ♦ c. 



tki .+ k?) Ht) - k 1 S 
To antidifferentiate the left side, write 
. ' u - (in + k 2 ) P(t) - kjS 

r * St *"C k «-.Ma> p, <«V 

% Thus,, we have 

1 fl du' «. «. 
T x + k 2 |u3t dt8 - t + c » . * 

and sxjice u > 0 we get 

' fc, i k 2 ln (( k i +> k 2)-P(t) - kjS) - -t + C. 
• To .evaluate C, we set t «- 0 and use P(0) = Si 

' ' i. 



In k 2 S *\C. 



Thus, 



k 1 hr 2 (ln ((k * + k 2> *m - k i s ) m k 2 s) * -t 



^ ^ |,t ' | 1 ;i"'; 1 ' 1 ;-' 

« 

To obtain Equation 4 we had to assume that P f (ti J* 0 
But this apparent restriction turns ou$ to be no 
restriction at all. For it is now a routine computation 
to show that the function P(t) given by Equation 4 
actually satisfies- Equation 3 for 'all t (see Erterci'se 2) 
and is^thus the.function we seek. ^Incidentally , it is 
also easy to confirm that P"(t).is never zero (see 
Exercise 2 again): V 

Exercises - * 

t. a) Obtain an expression for C(t) from Equations 1 and 4. 
• , b) Obtain the following expression for C(t) by first us<lng 
Equations 1 and 2: * 



I) - kiC(t) V 



, C« (t) » k 2 P(t) - kiC(t) = k 2 (S - C(t)) 

Show that this* result agrees wljth the. answer to^fa)'. * - 

2. <a) Compute P«(t) If P(t)*is given by Equation i 4. 

b) Show that, for this P(t), P«(t) - kiS - (ki + .k 2 ) 

c) Show that, for this P(t), P'(t) Is never zero. A 



J^V'#'^' V DETERMINING kr AND ko 

? u > 4>1 ■> Computational ^Preliminaries t , 



To determine lei and k 2 from EquatibhaJ^notice* thate^ 



1 V 7 ; ;., as ' 1 approached infinity /the express ioif^Karentheses 
*^r- approaches.!, since'V ( * \ + *2** approaches^! 



Therefore, *■ 

%T'. ; * - ts) v iim P(tT^ kis ; 

t-H» - *1 *2 • 




Le ! ^ s ^ s ^ 1 this rvalue* Q. /Dividing Equation 1 4 by Q x we' 
obtain v ^ v / 7*/ * . , \ 



.4 ; * \ " obtain 



(6)lh 6£^i - l)V-Cki * k 2 )t + In > ■* * 
♦ * . . • * • " • -I 

-v j The expression on the left si'de " ofj-Equation 6" is 

" a new function" of t. Let us -give, it 1 a Same: ' 

ay ■> • get) -in^ifsi-.i). 

Wow Equation 6 tells ^us that, according* to our model, 



W - » 8.W "-(Ic! +.k 2 )t + Jn k ■ 



pit 



SO * he S^Ph. of ,g ft), theoretically at least, is a straight 
' " - .... Xine.with slope ^(k, ♦ k 2 -) and y- intercept In |t. ? sb if ' 
; we coul-i. compute g (tO 'from the experimental data^iislng. , 
• Equation and then plot the points Jtr gOOl, We oould . 
, . ^cqoapfish two thing's": - ■ ' 

■ . . . " •-.-„ '•.•...'•.*•. ■ > V 

^v-;. ( ^ If the . k POirits come close to lyih'g on - a straight' . 

ISs- • ■ • line :v we^ould use 'this fact to co1ifirm;the Accuracy of . 

p& v . fC.'*^ 1 *;-;^' a11 ' "-'is the;jmodel which,, by Equations • 

fe-C'K.;-;- ^t*fP»^^A ; t^^he w pb*ints will "lie 'on a iin»: : 

felir ; • : r • : - r J^' • ^?^ht line; h'ns ^lop£ m and y- 
I'r- >c ; ".' V;' ' /'intixiept -b;/ i^'," could: write # #/> - • • 

«;.-.. (k^^k^ .r, •- \, ••. • .' 



lineC- 

» 



and then do a little algebra to find k*i and- k 2 (see 
Exercise 4) . . j ' n # ^ 



(a) Evaluate D = 1 im C(t).' 
v (b) Using (a),, evaluate £ . 



Solve Equations 9 for H x and k 2 in/terms of m arid'b.- • 



1.2 /Betl 



rmining Q 



There is bne tatch, though. Equation ,7 contains the 
1 4 *** ki S ' ■ 

k t .* k 2 ^ definition. , 



symbol Q. Now Q = l im P(t) 



But this is a dead eh'd, since we don't know kj and k 2 
in fact:, we are trying to determine 'them. • 

Luckily, tliere is a way-)out . " 'The, statement 
^ = me y^' ,in ex P e . r inental* terms? that the 

observed amoun/ ?(t\ of ions in ^ he plasina appfoa^Jies 

Q as an equilibrium amount/ Experimental evidence ' * 
confirms that it is feasible to^ continue monitorirtg the * ' 
plasma until P(t)* does^jiQt changej or change's very little 
with further: passage of time. *Wefcan then take Ais \ 
nearly* constant 1 value to be Q. ' ♦ ' '* j ~ 

4»3 An Example „ . ' 9 

• * Tosillastrate this method of determining ki and kji '* 
le^us considex the data in Table 1 for a hypothetical % / 
permeabilrity s1fe.udy. • \ 9 - . 



TABLE 1 



9 & — ( B 

Hypotheti cal Data for Permeabj 1 1 ty. Study 





' 0 


4 500 


,1000 


'1500 1 


2000 


2500 


3000 


[3500 


4G0O 


9 

4500 


5000 


P(t) r (mg) 


- 5-Q0 r 


2.96 


?.oi' 




i-rt 


1.01 


0.97 


0.92 


0.87- 


0.85 


0.85 



It appears that P(t) is settling down to an equi- 
librium of about 0.85 mg, after t = about 4500 minutes. 
We shall take Q = 0.85 mg. 



We can now use Equation 7 to calculate*- the values 



of g(t): 



We get 



g(t) = in ( Z^l* - 1) 



In ( 



P(t) 



t 


0 ' 


* 500 


1000 | 1500' 


2000 


' 2500 


3000" 


3500 


4000 


9(t) 


J. 50 


0.91 


0.31 }-0.28 


-1.08 


-1.67 


-1.96 


-2.5G, 


-3.75 



It JLs impossible to Compute • g (4S*00) or g(5060) on the, 
b^is ©I our 3ata, -because ^ln ({jfff 1) = In 0 is not 

\ dejfined. This difficulty arises because our data are 

given to only^ two decimal* places . *With more precise 
j measurements we might have found, for example, that 

* P(450Q)*=. 0,854 and P(5000) 0.851, ' • 

>^> 

ha 



Jhe points (t ? g(t)) are plotted in Figure 2. They 
e closo'tq kyi^ig "on a^ straight, line.. In Figure 2 we « 



have drawn in, *by eye, wjiat appears to be- the "best- 
fitting" straight "line. In doing this- it is wise to 
u^ a transparent straight edge," so our view of the 
points is not blocked. -We . have tri^d to. draw the line% 
so som£ 'of* the 4)oint§ are slightly above it, some 
slightly b^ow 1t, # and none too far fronr&t/ 

. The line 'we have drawn appears to have its V 4 

* • •? * 

y-intercept at-a^bout b,- 1.55. It. also passes pretty 

nearly* through "(500 „ 0;9) and (400, -3.S). Therefore its' 
slope ^m is about ^ I ggg * ^ n L26 x J0"\ 

* ^ % We cotfld now use jpcyiatlons 9 'to* find kj and k 2 , -But 
if y.ou have "done " Exercise- 3 you haye* discovered ' tha£* ' '* 
Equations- 9 can^be rewritten' 



* t ^ Figure- 2. 
Plug'ging in, b- = l y 55 and m * 1.26 x 10* we get 



.3 



kj = 2.21 x 10" min " . / / 

k 2 =/l.'04 x 10* min' " . , ■> 

4.4 Some Comments on the "Beist^Fitting" _L*ine » * 

(a)- 1^ is not important Whether the lifte'pafses 
through any of the given points. In fact it would be a 
mistake simply to draw the line ^determined by two -of the 
points. We might be, unfortunate enough, to pi£k two 
points 'which are inaccurate because of experimental" 
error or rouijiioff error. ~ . 

1 (b) * In finding the s^ope. of the line; use^two 
poinff^on the line itself, ratfier than two .of* the given' 
p'oTjrts^,7*As .mentj.oned/in 4 . 4 (a) fome of the* given 
F^oint^s are^bounrf* to be Vlightly off. The line in effect 
"average's out". these errors. ' * % . 

( r * (c) Al a $o in finding the Slope, use. two pbints 
fairly far. apart, -If they are close together, then, the 



denominator in the slope formula will be small, and a. 
slight error in reading the coordinates of one point may 
result in a huge error in the slope. 

5. EXERCISES " < ♦ ' , 



5. In a permeability study the function g(t) = fn - 1) has 



been computed and plotted in Figure 3- Determine k x and k 2 for 
this experiment. * \ \ 

* * - - 

2 



'cn 
O 



-I 



-2 



iit 



rr 



tit: 



4-L 



■4-H 



r- u -r 



■M-t-i 



ffc 



e+3 



4-h 



' 1 



-H- 



t ; 



4p 



500 ^o<ro ]>500 2 00o 

» time* (in minutes) 
' Fiqure 3. 



6', Determine kj*and k 2 from the following data for a hypothetical 
permeability' study. / • ^ 







3 


4 


? 


6* 


7 


8 


3' 




2.05 


KB 


.1,28 


1.20 




J .13 


K12. 


1,10 
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Consider a two compartment model, not*closed, in which P(t) is 

ma/ntained at a constant leveMby continuous adjustment from the 

ouAsMe, for example).. If P(j) = P Q the determining equation for 

C(t) becomes * , t . • 

i — ± * 
C r (t) ».k,P 0 - kiC(t). . 

(a) Find C(t) if 0(0)*= 0. 

(b) Find D = lim C(t): 

(c) If Q d n m p( t ), find jL. 

In the models presented J,n £<ercj se 6, find fche.time at which Jb* 
C'(t) tsj^s oo a maximum, and compute that maximum, . 
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6, ANSWERS TO EXERCISES 



''Jt 



kj 



2. (a)^ P'(t) = 



k,S 



kj + k 2 



- (k, ♦ k 2 ) e - (k l + k 2 )t 



^) kjS - (kj +, k 2 ) P(t)' * 

= k*? - (kj + kaJ^rc-^ 



Ui + k 2 ) 



! e "(kl + k 2 )t" 



* kj 



kl s - kjs i +5f 



ki p "(ki + k 2 )t 







'' positive* 




3. 


Ca) Je,S ' 
, kr + k 2 










T 


kj-^SU-T 
l + e b • 




5. 






6. 


kj.* d. 14' fir* 1 




7. 


(a) j» 0 <1 - 


/ - 




• k? B 
(b) kf p o ^ 


■r •••> 







-k 2 S e " (k l + >2)t , 
(c) All three of the numbers k 2 , S, arrd e~ (kl + k 2 }t ^ 



1 + .eb 



f ■ IJaxfmuif^ccurs at t - 0 (endpolm maximum) a*f equals k,P„ 
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Return to: 
EDC/UMAP ' 
55 Chapel St." 
Newton, MA 02160' 



Student: ' If you have trouble with a specific part of this unit, please^ fill 
out this for* and take it to .your instructor for assistance. The information . 
•you give will •'help the author to revise the ' unit. 

Your Name • * . Unit No. 



Page 
*' 0 JJpoer 
"0 Mid die 


OR. 


Section • 
Paragraph . 


OR 


0 Lower 
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0 


1 
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Model Exam 
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Problem No. 



Instructor : Please indicate your resolution of the difficulty in this box. 



o 



Corrected errors* in materials ♦ List ^corrections *here: 



1 , \7J Gave student better .explanation, example, or procedure than; in unit. 

Give brief ' outline of your addition here: ' ". 



fc. 



■ y 
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Assisted s^jident in acquiring, general- learning and problem-solving 
' skills (not using examples from this unit.) 



" ft 
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Instructor's Signature 



•V 



. * • •. ' Pleases use reverse if necessary. 'CO 
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<} 

Date 
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Check tHe. choice for each question thkfc comes closest to your personal opinion. 

1- How useful was the amount of detail in the unit? 

p* , " , » . ^ — - * 

* Not enough detail to understand the unit 



_Unit would have been clearer with iiore detail 4 
^Appropriate amount of detail. . \ **P 

JJnit. was- occasionally too detailed, but this was. not distracting 
J^oo much- detail} I was oftep distracted 



How helpful were the problem' answfers? 
• Sample solution^ were too brief; I could not do the intermediate steps 



Efficient "information was^ given to solve the. problems 
^Sample solution^ were toc(^det£iled ; I .didn't need them 



3. Except for fulfilling the prerequisites , how much did you use other sources (for 
^exam ple, instructor, frienjs, or otjfer books) in order to understand the unit? 

^ A Lot . Somewhat A Little Not at all 

~ » * , * 



4. How long w^s this uni'f in, comparison^ to the amoiftxt of time you generally spend on 
a lesson (lecture and homework assignment) in a typical math or. science course? 

_ ° 

Much^ Somewhat [ About 0 Somewhat 0 Much - 
^Longer b Longer ' the Same / Shorter <* S horter 

*■ — \ 3 * 

5. We're dny of the following parts of the unit confusing or distracting? (Check 
' as many lifs apply.) . jy* 

P rerequisites %* e 1 

Statement of skills and concept?, (objectives) 4 % 

Paragraph headings _ . 

E xamples ~ 

* J _^ ^Special Assistance Supplement (if present^ 

, Other , please explain - f 



Were afty of the following parts of the unit particularly helpful? (Check as inany 
*s apply . ) , 

P rerequisites • * 
Statement of skills and concepts (objectives) 



^Examples 
.Problems 



-/ ; P aragraph headings* ^ # • ~* -> ^ 

.Table of -Contents .... * ^ ' ' 

S pecial Assistance Supplement (if present) ■> 
Other, please explain t . * 

* Please describe .anything in the junit that -you*did: not particularly like. 
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1. GROWTH MODELS 

The behavior of the economy of a given society 
over time is of interest not only to the economist, 
but also to all citizens living in that society. 
After all, the increasing abundance or t scarcity of , 
jobs, goods and services depends crucially on how 
fast the economy is growing relative to the ^mderlying 
population. Thus it is not surprising that words 
like growth, stagnation- (and* more recently stagflation), 
etc, have become commonplace even in the everyday , * 
world of the nightly news broadcast. 

For their part mathematical economists have 

developed'many different models, called' growth models* 

to describe the expansionary processes in an economy. 

\ * 

In this unit we will be studying a particular 
model of growth, applicable* to planned economies in* 
Khich sl\1 means of production are socially owned.. 
It was developed by the Russian economist G. A. Feldman 
(1928) in connection with planning for the centrally 
controlled Soviet economy. 

The purpose or this model i$ to describe the be- 
havior over' time of a two-sector economy in which sectoral 
investment allocations are controlled by a central authority 
according tp an overall economic plan. 

The impetus for the construction of Feldman's 
model came in 1927 when the Soviet- Union embarked on a 
sequence of 5 -year plans for the^ifxpansion of its 
economy. The devastation, caused -by the First World War 
and the following Civil War nad, to a large extent, been 
overcome and prewar production' levels had been restored. 
The first publication of £he model occurred in the 
November 1928 issue of the-Soviet State Planning 
Commission's journal Planovoc Khoziastvo (The Planned 
Economy) * k » . I 
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. v * ~ Altndugh3the model 's direct impact onYSoviet plWing 
policy formulatio* is questionable, there is np dqubt/^t 
its 'flavor 1 is wejl iiKJceeping with the initial sequence 
of Soviet 5-year plans with their strong emphasis, on 7 
building up the heavy industry sector of the, Soviet ' 
economy. >^ ' 

*>As»an indication of the model f s durability » we~ 
note that its main features were duplicated in ''the 1950 f s 
(apparently independently of Feldman's wtfrk) by the Indian 
ec6nomisf Mahalonobis in his work on a planning model' for 
the Indian economy^ 



2., DEFINITIONS 



'J 



2.1 Rate of Output and National Income 1 
" As is usually done in Ma/rxian economics we divide 
the economy .into two sectors/ the producer goods sector 
(Sector I), producing goods /to be invested in production 
in both sectors, (like drop/ forges , milling >macl&nes and 
tractors), and the consumer! goods sector (Sector II), 
producing gopds, to be consumed nonproductively by the 
populatipn (like bagels , 'television sets,' and footballs). 
We assume>ahat both sectors a.re; producing a continuous 
stre-am \f goods, / ' " 4 



I Before we can model/ the quantitative behavior of 

^ the outputs in "these. sectors over time, rt we must first 
define a quantitative measure of how much* each sector,' 
is producing at any gi,ren time. This we proceed to do 
, as follows: let J(t,jf) stand for theTiet output /measured, 
say, in dollars) of. Sector I' between time t and, time, u 
(measured, say, in years with an appropriate point in 
^ftime chosen as t ■ 0), 

•s f ' * - i 

The "net*; here means the output remaining aftjer all 
o the »>wear and tear 1 * in the investment goods being used • 
in the two sectors "has btfen macfe good from the "gross 11 
output of Sector I. .We use the net output , jsince only 
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this part of the output of Sector »I can be used to expand 

the economy ^over time. * The rest t>f .the output of Sectoral, 

i.e. the- difference between the "gross" and "net" Outputs, 

is required just to keep' the economy at the level it had 

already reached, ' ' * 

• 

If we agjee 'that ail of , the net output of Sector I 

* goes for investment^ it seems appropriate %o call the 

* quantity 

i' (t , u ) = iiUil 
' . av^ ' J u-t 
< • 

the (average) annual* rate'of investment or the (average) 
annual* rate of output in Sector J. over the time interval 
from t.to.u. (Note; 'the units here are dollars/year.) 
The ( instantaneous ) annual rate of investment is then 
defined as * 

i(t) * lim I (t,t+At) = lim «J(t,t+At) . ' 
4 ay At-0 — At ' 

'The quantity I(t) measures what the annual rafe of 

output in Sector I is at ,a given instant t of time. 

Sjmila.rly we can, and do, define C(t), t he 
(instantaneous) an'nual rate of output in Senator XI. 
The sum- of the rates of output in the two sec'tors 
is v then galled the national income: y(t) = I(t) + C(t). 

2 . 2 &- An 'Analogy * • , 

* * Xpu may have nitt'iced some similarity between 

the.procedure^we ijsed to -define national income-.and 

the one usfed e to define instantaneous velocity in 

first semester Calaulus or Physics courses. Your 

study of velocity probably began by def inin^k * - 

"distance" or^ "position" function s(t) which described 4 

the location of an object at time t usually as Tf^ # 

the direct distance from some arbitrary *ppinjfc ^called 

* it 
zero. It then defirfed the average velocity over the 

time interval frpm ( t to u as ' 

s(u)-- sCt) . \ * a / • 

x u - t - - ' 
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This average velocity ,and our average annual rate 
of output each amount to a ratio of two changes or 
differences. 'The numerators look quite different, but 
they really aren't/ In studying velocity it was pos- 
sible to define a function s o'f time whose change from 
time t to time, u could be calculated by subtraction 
and was just what you needed. There does not seem to 
be a convenient way to do the analagous thing here. 
Nonetheless, saying "J(t,u) is the net output between 
time t and time u" is very similar to saying "styf), - 
s(£) %s the net motion (change in position) between 
time t and time u." 

Thus natioital .income can be thought of as the 
derivative of a (somewhat fictitious) "cumulative 
total output" function just as velocity is the deriva- 
tive of the "cumulative total distance traveled" func- 
tion. Table 1 further illustrates this analogy,. 

# 

TABLE 1 * \ 



Cumulative 
total over 
a time, 
interval" 


Jypical 
units 


Average rate 
over a time 
interval 


Typical 
units 


Instantaneous 
rate at an 
instant of time 


T 

Typical 
units 


distance 
traveled' 


meters 


average 
velocity 


meters 
sec ■ 


velocity 


meters 
*sec 


total net 
output in 
Sector I r 
P-A) . 


dollars 


average rate 

of net 
investment 


dollars 


I (t) = .rate"* - 

of net 
investment 


dollars 


year 


year 


total- out- 
put in 

Sector II * 
(=B) 


dollars 


average rate 
of* output 
in Sector II 


dollars 


C(t) = rate 
of output 
in Sector II 

* 


dollars 


year 


, year 


total ' 
output 
(*A+B) 


dollars 


average 
national 
income 


dollars 


Y(t)=I(t)+C(t) 
= national , 
income 


dollars 


year 


year 



♦ > 
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THE MODEL 



3.1 Assumptions 

« With the definition of I(t) and C(t) in hand we 
are' ready to start building Feldman ! s mpde*l. 

We let Ip(t) and I c (t) stand for the net rates of 
investment in Sectors I and II respectively. , The en- 
tire oiltput I(t) of Sector I is to be invested in Sectors 
I and II. The central planning authority decides how to 
spl^t the investment pie betwe'qn the two sectors and 
allots a constant (i.e., independent of time) , ^positive 
fraction s o*f the output of Sector I to. Sector I (see 
Figure 1) . 



> 

* / 


* 

sl(t) 


Producer Goods 




Sector^ I 





(1-s) Kt) 



Consumer Goods 




^ Sector II 


C(t) . * 







J 



Figure^l.* The net output of Sector I is produced at the 

rate af I(t) dollars per year. It is channeled 

* back into producer goods at the frate of sl(t) 

dollars per year, and invested in consumer 
goods at thyrate of (l-s)I(t) dollars per year. 

#1 s . - 

° Nqwv we need ttTexplore the r|lationship between I(t), 
C(tif)<Ip(t) and I c (t), all of which are assumed to be 
djifferentiable functions of t. We" let'*K p jft) and K c (t) 
stand for -the capital stock, l.eT*, the t quantity of invest- 
ment goods invested at, time t in the Sectors I and IT re- 
• speStively. By -assumption, there is a stable relationship 
between the capital stock and the corresponding output in 
each sector,' i.e., capital "stock is proportional to output; 



U) xhy~ = v p r constant or K p (t) = v p I(t) 

C'Ctj = V C * Cons tant ;or K c (t) ^ v c I(t) 

The constant^ v p and v c in equations in (jL)and (ii) 
are called tjre- capital ^output ratios (for obvious reasons), 
For the sako^of simplicity, we assumeXthe constants have 
a common value v: 

v p = v c ^= V . . , 

If we differentiate both sides of the equations in (I) and 
(ii) and note that the rate -of increase in the capital 
stocje k is precisely the corresponding rate of investment, 
we obtain: . r * 

I p (t) - K p '(t) - v • I'(t) 

» and \^ 

i c ('t) - * c ' (t) = v • ecu 

t 

We can suiwuarize the* discussion above in "the fol- 
lowing two assumptions of the peldman. model : 

** . * 

Assumption 1: There exists s (0 < s <"l) 
y 1 such that I p (t) *' sl(f) and * . 

r c w s (i^) nth 

Assumption 2: , I(t) andjc($j are different! able 
1 » * ' / 

and I'(t) = | I p (t) and C'(t> ? /l c (t) . ' V. 

The crucial (and distinguishingVassumption of the 
model is Assumption 1. The claim iZ makes would almost 
certainly be false in the absence M controlled alloca-- 
* tion according to a plan. ' It «al<?o focuses on the most 
important parameter of the model, namely s. For, as 
we snail see in the next section-, it is the choice of 
s that decides" how fast" the economy grows, both 
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overall (Y{t)) and i^ terms of a its sectors (I(t) and 
CU), 

To avoid any confu\ion let unstress again; in the 
model presented above, s^and v are vtime - independent 
constants .(parameters), wHereas I p (t) ^ l Q (t) , I(t) and 
C(t) are increasin g functions of the triable t^time). 

Finally we note that Assumption 1 can be modified 
to .allow the possibility of % * 0 or s - 1,* in which 
cases one of the s_ectors does not grow at all. $oth the 
statement and~derivation of t^e resi^Us of the model in 
^the case s = 0 or s = 1 are considerably simpler thari 
*in the case 0<s<l. . We ask you to g\ through one such 
( derivation in Exercise 12, \ 

S.*2 Derivation of Results * . r 

We suppose, at" some starting point in time (con- 
veniently taken ?to be t = 0) the values of I«, C, and 
\^ff known: * ICO) « I Q , C(0) = C Q and Y(0) = Y g . 

x . Now the equations of Assumptions 1 and 2 can 
^be combined and rewritten as: , * *- 

en p(t) = f Kt) . • . 

and _ ' ' 

,P \ 

.Dividing -both sides of Equation 1 by I(t) and 
integrating we obtain 

t 

1 fdt 



c;( t ) -i^s. i( t ). 





For notational convenience we Refine k - ^ and obtain 



I(t) =.I Q e 



k£ 



•Substituting Equatlpn (3) into Equation' (2) 
' and keeping in mind that k = ^, we obtain aften some 
algebraic manipulation . , ' ' . 



C'(t) -.i^A. e 



= ll& 5. i e * 
-%1-s x , kt 

Jtow we~pan integrate eas'ilytfo get an expression for C(t) 



C»(t)dt.*,±§5. i ( 
0 ' \ 



u kt 
ke Kt dt 



C (t) - c 0 -»i^i 0 



I-s 



.kt 



T4) «*. C(t)'- ♦ ^2- I Q (e J 

Adding Equations (3) and (4), yields 
Y(t) - I(t) i C(t) 



1). 



.kt 



.kt 



since + Cq s Y Q , 'and finally • 

I 



1-S x ^ ,_kt 



(5) 



The rela^ve rate^ of .growth of a diFferentiable quan- 
%ity*i(t).is f'(t)/f(t). '"It measures *not how fast the ' 
♦ quantity is changing in absolute- terms, but rather how 
fast it is changing Relative to its owiv'size*. "*From 
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Equations (1), (2), (4) and (S) we can obtain -expressions 
for the relative rates of groyth for I (t) , C(t) and Y(t) : 



■(6) 

(7) • 

(8) - 



Y££y- = f = k\ - . (directly from (1)) 



C'(t) 



ro s i 

17?: it? 



e* kt + l 



(Exercise 1^ 



Y'(t) a 



1 


c 


y s ; 1 

ro j 


* (Exerci'se 
e xt + 1 











Exercise 1 

(a) Prove Equation (7) 
> * 



•(b) 



A, C ! (t) s 



Exercise "2 
(a) 
(b) 



Prove Equation (8) , 



* 1 



Show ^[!^ 



+ — as t ■* °° . 



3. 3 \ The Average Propensity* to Save * 
* * * . \ 

Another important) quantity is q (t) ,**de fined as the* 

ratio 'of investment (rate of output of Sector I) to .the 

national incpme . 



Thus a(t) measures the fraction of the- tota^ output 
of ^he economy, at a given point in time, which is saved 
(invested^ rather than consumed . So its name, althoug 
cumbersome, should not conte as ^a surprise: ct(t) is* 
called the average propensity to save . 



If'we let a 



0 



o(0) 



the ratio of I A to Y n then 
« % 0 0, 
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f.- i 



«(t) - — r-T 



Y 0 



. . . V : i) e -« + 4 • . 

y * • ^ . - 

• If s > a Q one can shoV (Exercise 3) that^ <*(t) is*an 

increasing .function of time. In any case, s*ince 
- kt * 

r> e- + *0 as t * », .we see that 

* * 
» lim a(t) » s. 

That-.is, th^> average propensity to £ave approaches the 
fraction of investment devoted to Sector I as time- 
goes on. , • m i* 



? Exercise 3- 4 

• Calculate a'(t). Use your* answer to show that 
.(a) ^aft) is an increasing function if s > a fl : 

4 * 

Cb) -a(t} is a depreasing function "It s < o Q %- 
• Exercise 4 * , , 

^ • ' Compute a(t) for t = 5 and t = 10* if 'a Q = 0.05, 
s » X).5, and v = 5. ' / N ' •* * ? . \ 



t . - 4. NUMERICAL EXAMPLES 

t ,In*this section we provide 'we results of calcula- 
tions of. the relative rates of growth of I(t),£(-t), 
and Y(/t) and of the average propensity to save a(t) 
for some reaSortabW^ values of the parameters of the 
model; namely, v « 1, 5; a Q = 0.1* 0.3; and s = 0.3, ,0.7. 

It"is certainly not unusual 'for an- economy to be 
reinvesting 101 « 0.1) orHo% (a Q = 0.3) of its ^ 

' - \ -io" 




total output. As fo^ the reasonableness -of v - 1 or *5 
(when time, is measured in years '[see Exe.rcise 5j), We 
can refer to empirical determinations of the marginal 
capital coefficient/ In fact, (with time measured- in 
years), Feldman estimated the capital-output ratio v' 
to be 2.4 in the years 1928-33. Leontief (1939) fo^nd 
comparable* ratios to range from 0.076 to 7,1 in 
various branches o£ the American economy. 

We have^tced one high (Table 2) and one low 
(Table 3) value \of s in the computation? to^llustrate 
the effect of favoring one or the other of the two 
sectors for investment. 



T . • . . < 11 
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TABLE 2 



Sector I Favored f^r Investment (s = 0.7) 





( 


f s 


. a ( 






' V 


= S 


0 


0.3 


1 


t ' 

[years ^ 


r/i 


Y'/Y 

(%) 


cvc, 
(%) 


a 


171 

• (%) 


Y'/Y 

(%) 


C7c 
• (%) 


a 




1 


14 


2.3 


0.8 ' 


0. 11 


14 # 


6.S " 






% 


S 


14 « 


3.S * 


1.3 


0.18 




0.** 


4.4 






20 


14 


"10.3 


6\3 


0.S1 


14 


13.0 ' 


11.0 


0.6S 


• 


SO 


14 V 

1 


13.9 


13.8 


0.70 


14 


14.0 


.13.9 


0.70 










f 
















v = 1 


r; 


= 0.1 . 




V 


= 1 


a 0 =D.3 




1 


70 


17. S 


6.4* 


0.18 


70 


42.1 ■ 


21.8 


0.4*2 




s 


70 


S?.3 


43.6* 


0.S9 


70 


67.7 


61.7 


0.67 




20 


70 


70 


70 


0.70 


- 70 


70 




0/70 




SO 


70 


70 


70 


0.70 


70 


70 


70 


0.70 






i 


c 
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TABLE 3 




\ 










Sector II Favored for 


Investment (s 


= 0.3) 






* 




v = S 


a o 


= 0.1 




> • 

V 


= s 


a Q p 0.3 


- 


t 


iVi 

r (%) 


Y'/Y 

(%) 


,C7C 

(%) 


• a 


r/i ' 
'(%) 


Y7Y 

(%) 


C7c 

(%) 


a 




1 


6 


2.1 


1.6 


0.10 


6 




6 


0,30 




s 


6 


*2.4 


1.9 


0.12 


6 


6 I 


6 


0.30 


t 


20 


6 


- 3.7 


3.2 


0.19 


6 


6 


6 


0.30 




«so . 


- 6 


- S.S 


S.3 


0.27 


6 




6 


0.30 


























v = 1 


a o 


= 0.L' 




V 


= 1 


ot Q = 0.3 




1 


3b 


12.1 


9.6 


0.12 


30 


30 


30 


0.30 




s 


30 
*• 


20.7 


18.3 


0.21 


30 


"30 


'so ' 


0.30 




20 


30 


29.9 


29. 8, 


0.30 


30 


30 


30 


.0.30 




SO 


30 


30 ' 


30 


0.30 


30 


30 

• 


' 30 


0.30 
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Exercise 5 -* \^ 

(a t l What are the units of v? 

(b) ^ How does t v get affected if we switch from dollars 

to dimes as the measure of output? 

(c) How does v get affected if we switch from years po 
months as *the measure of time?. 

Exercise 6 f * 

— ■ ■■■ ■ # 

t » rt* vi f 

A^surae^Q - 0.2, v =* 3, Evaluate ~ , ^at t 
% for d 

(a) ' s - 0..1 . . 

(b) s = 0.6 

Cc5 s = o> 



5'. CONCLWING REMARKS 

In concluding, we makfe three remarks* on this model : v 

1) I 1 /I is given by ^, and does not depend on a Q ; ' 

i.e., it is'qujlte possible to have a very fast 

growing sector I, while £he initial average 
> • * • 

propensity to save is low. This seems to have 

happened in the Soviet Union. 

2) The relative* rat e,s of. growth of the * natjLonal income 
and the consumer goods sector* eventually approach 
the growth of the producer goods sector (Exercises 

, 1 and 2) . ' ' ^ » / ■ 

'3) The average propensity to save rises to s (provided 
<Xq < s) . The' empirical verification of this predic- 
tion of the modtel seems to have been a source o!k . 
controversy among economists. Different analyses 
of the Soviet economy have led to widely differing 
estimates: from a' regaining essentially constant 
at .'23 in the years 1928 *to 1937, to a increasing 
5 from .17 to .37 in the same time period* (See 
Exercise 11.) * „ 
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EXERCISES 



Exercise 7 
Given ; 



C Q - 2I 0 , s = 0.75 and v< = 3 



(a) How many years must pass before the value of I(t). 
catches up with th^ value of C(t), 

(b) Sketcfi the graphs of I(t) and C(t) 



(Assume I n ■ 1.) 



Exercise ^8' 



* Under what circumstances will - 
' * . * C'(t) 9 Y'(t) _ lUt) 

for all values of t? • 

(Hint: First see when ^ifl . i^l , and then 

( 



Exercise 



What is the relationship_between y^ffi - and a(t)? 
Exgrcise 10 * 

In Feldman's model is it possible -to have 

(Hint: First see when <j$l < '.and * ': 

then when fjgl < ^01 0 ' ' • 

Exercise 11^ 1 
* 

Discuss possible reasons for t^e mistaken estimation, 
of a(tj if I(t) and* Y{t) are measured in current irices 
and the Inflation rates in the two sectors" are different 
Exercise 12 

# Suppose. s- = 0 and derive a formula for C(t). 



14 ' 
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B^ercise 13 

It takes not only capital but also labor to produce 
output v What are the implicit assumptions of the Feldman 
Model 'about the supply of lab<?r? * 

Exercise 14 



What. are the implicit assumptions of the Feldman 
model about international tradfe? 

Exercise 15 <• 

Suppose we drop the assumption v p - v c = V, i.e., 
Vp and v£ are now two different constants & How does 
this' modify Equation (4)? 

• 

. 7. REFERENCE 
• * • 
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8. ANSWERS TO EXERCISES 



* 

t t , ks ( s / tt p - D^ 1 

*• ° l*J = TT - " rr* — is positive if 

tes/« 0 - l)e Kt + i]« 

s > ct 0 , and negatiw if s < a n . 

„4. a(5) = 0.077 , * 

a(l f 0; = 0.116 

5. (a) Unit is time; e.g., years, 
(b) Stays same. . 

(cr) *v(months) = 12 v(years) 

6. (a) I VI - 3.331 C'/C - 7.21 Y'/% = '6. SI 
(b) ' IVI = 201 C'/C = 3.91 Y'/Y = 7.61 

" ; (c) I r /I t 301 C'/C = l.l| Y'/Y = 8.41 

7. (a)^ 3.66/years 



8. a Q = s 



9. a(t) - v • 01 



10v No. 



12. I'(t) = i I(t) = 0 
.'.I (t) = I Q = constant 

c(t) = Ls-s i( t ) = b. 



C(t), = c 0 ♦ ^t 



IS'. C(t);- C 0 + ^^l 0 (ev-p.!) 
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PREFACE 

In your mathematical 1 ! studies up to this pointy you 
have often been ^called upon to make a mathematical 
"description" of a situation. This description usually . 
consisted ^f an Equation of some sort. The familiar 
"wprd problems'' in an algebra course gave rise, to one 
or more linear equations, or perhaps a quadratic or 
exponential equation. y 

When the! situation to be described contains a non- 
uniform rate of change, the equation will contain a 
derivative. It is then called a differential equation. 

It is the purpose of this module to show how to 
describe .certain physical situations by means of differ- 
ential equations and how to solve these equations by 
simple graphical or numerical techniques. 

• The differential equations (DE's) 'treated in this 
moduli will^be ordinary (containing no partial deriva- 
tives) order (will contain- only first derivatives), 
and of firet degree (the derivative will not be.^raised 
to any power higher than one)$ ^ ' ' 

I would like to express ^appreciation to Mary Jane 
Neuendorf fer ktid tfilliaa U. Walton\pf Project CMC fox 
their extensive help in preparing this module, and to 
the many peopyeUrho reviewed the draft and offered 
valuable suggestions. 
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' ' • ' » ' ■ Chapter 1 * ' .» 

THE OPTICAL FILTER PROBLEM 

I threw down my pencil in frustration. Taking my 
action as a signal, Polly and Herb did^Jhe same. 

. Professor Arclet. didn't notice.. He haxf do.zed 1 of f ; \ 
his habit When* he gaye\a test,- and he was snoriifg iightly 

•Bow, forming silent equations with, his lips. * 
* .. * * * * 

Herb- held the book poised above the floor. .His 
chubby face showed both discouragement and apprehension. 
He looked at mje with a questipning smile. -That book. 
Another reason fyr our discontent. Heavy. Complicated. > 
Eighteen dollars. I knew it was unkind, but I gave Herb 
the nod, * * 

FiVe pounds of calculus came crashing .to the flooV. 
r Arclet sprang to his, feet, banging , his Jcnees painfully on 
the desk. He staggered to the .blackboard and began to 
lecture, a' 'continuation of- the derivat;i'o* begun in h*s 
sleep. J , - 

••Professor Arclet," 1 interrupted, "We cari't do 
this- test." . J ' j 

'The exam is whai he calls a prq-teat. It was the one" 
.on differential equations. Arclet always gave Us" a pre-" ' 
test When we got to a new topic, and thoselof us that ■ 
-passed didn't 'have to come to 'class- until ,the next, topic 
caae up. 'i neyeV passed any of these tests, but tkey did 
give me soine ideas of .what was cooing, and what l.was • 
supposed to -be able to-do later on. Sometimes he used the 
same test aftpr *fe. finished' the topic. 



"Well," he said, '•rubbing his sore knees and absently 

staring opt. at the Vermont landscape,, "that means we need • 

to go .over this material. *. ♦ , . • 

^ . ^ j * * » 

There was V groan from Polly; who had freckles and 
wavy auburn hair. "But professor," she said, "differ- 
ential equations seem So, so— nowhere! 'We spend all that 
time and effort learning this stuff, and -I 1 11 bet we never 
use it in a million years." , Polly's face was flushed and 
her voice 'shook. , V * » 

"Yeah, yeah," from Herb and me. 

Encouraged by this support, Polly" went on. , "What is 
the stuff .good foT? >Can you%give us a single. example of 
^_where soraeone'used <k differential equation to solve a real, 
important problem?" / * * 

! , Polly sat down to our applause. . . 

* Arole.t was silent foj 4 longtime, his thick, bushy^ 
brows drawn -together in* cpncentration. This was' our first 
open rebellion, and ^wondered how^e would handle it. 

He shuf-fled around to the fifrn* of the desk and sat ' 
on it. We krtew this was symbolic* -removing the barriers ' 
between student? and teacher. This usually meant that he 
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was going to tell us not to think of ourselves as students 
versus teracRer, but as'a groups of pec^le^seelcing the truth 
together. Instead, he just satjthefe, kicking his heels* 
'against the side .of the* desk. 

, "Let me tell you^a lit-t^e story," he began.after a 
miriute or two of thumping. "It's 3 " about a young friend of 
mine, named Denis DropmoVe."* 

V settled cSmftortably "into /my seat, ready to doze. 

1 /'Young- Denis," Arclet said,- "was a proposal wxiter . • 
at the deadly Nightshade Sunglass Company, tack in the * 
•sixties. 1 One day, thj company # *eceivea a request from 
^the^ Department of %fense to prepare a Quotation* for a 



* badly needed; item/ Some i^eirtioned^ foreign power had m ' 
developed a new weapon, -the ballistic Laser for Infantry 
Neutralization and Destruction, code narae^BLIflD. Denis's 
company was asked to pr.oduce a protective eyecover to bet- 
worn by soldiers: Sunglasses*, Head Attached, Defensive, ^ 
Expendable, Shatterproof, code name SHADES., 



\ 




— / 

"The *)OD\Kad sent along a small sample of the -lens 
^material to, be 4ised'/ Only %he thickness could be changed. 
Another requirement was that the glasses reduce the ihten- 
saty off the laser beari to JO, percent 6^its^value. ' 

* w "Testing* the' sample in the 4aboratory, the company 
^■lsco^vered that th# 1mm thick materia? would reduce the .* 
. iriteasity of *a laseir beam by IS percent, 'Well, good, \ • 
^said/Al, the cpmpa^iy's softspoken" chief 1 engineer . ' r If one, 
.?millimeter*Sf that material^will remove 15 percent, then 
-six millimeters WilP remove the requir^jj' 90 percent ^, Denis 
write^nto the proposal that we will make /the lenses,, 
six millimeters thick. 1 . * : ~ - 

"Denis chuckled. ' 'Al, \ he ,said, loud enough for 
eyeryone ^n *he office to hear, 'yo\i doiPt know^ any thing 
about optical filters, -do you?' * \ * 

+* - ^ "Al smiled* lip at Ddnis. 1 You 're. fired, ' he said'""' * 
softly," ■ ** '■ 

e jThe clas£*be|i rang'. Arclet's stofy grabbed us, but * 
rso^did' lunch; We aU ran out. * " ' ^ . 



^Chapter 2 
THE SAGGING BEAM PROBLEM 



► Tod<ay, Polly, Herb and I were not la^e for class for^ 
the fiyst time in*the semester. Professor ArcLet wa.s 
Seated at the edge of the desk, ^exactly as we had left him 
yesterday-. 

"Denis Dropmore sat drinking beer at home," Arclet 
bejan, "dreading the return of his wife. How could he 
fcreak the n*ws? At 5:30, Angelica Dropmore came home, r' 
She looked distressed. Denis took his young wife by the- 
hand. 'Honey, 1 fte stammered. 'Today I . . /;' 

"'Denis,* I got fired,' wailed Angelica. . ' 

_ "Denis was what^ypu'd call 1 flabbergasted. Angelica 
poured out he*r story. 

"'I- was at the bank, as usual, today ,* tidying up^iny 
desk after the last customer left, and couldn't help hear- 
ing Nichols and Dymes , those two" creepy vice presidents, 
having. a' battle. You know the ones--always trying .to 
outdo each other. o • 

"'Nichols shouted, "Under the safe:." 

"'Dymes screamed, "No, you twit, it goes .under the 
midspan." * ' 4. * \ 

"Angelica raced oh," 'It "alidad to do with that heavy 
safe coming" in ^next week.XThey found a location directly 
over a> thirty -foot -long steV; beam, but couldn't get it 
closer than six feet from the end* support. -Old Mr. Usury, 
the bahk president, was afraid that bending of the^eam- 
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•under that massive load would cause cracking of the * 
e^xpensive marble floor tiles.. He tpld the two VP'*s 
» to install a third, vertical coiunjn under the beam at 
tjie .point where tlie sag would be greatest if the third 
column -wasn't* there:' Nichols wanted to put it directly 
under the safe, while Dymes said no way, that the great 
est sag would be at the midpoint of the'beam. 

Remember, 1 Denis, I studied structures in archi- , 
tecture school, so I Jcnew better. 

v,,, The v column should* go between the safe and the 
mid^oiftt, 1 I blurted. *J / • 

; ,M Two pairs" of eyes turned/ on me. ' Disbelief and 
'* contempt, I felt it ii^ their stare. Nichols and Dykes 
spoke in unison, agVeeing with .each othe^r for-the first 
time since they joined the bank. * 

0 M, 'You*re fired,!' they said."' 





Chapter 3 . 

C 

9 

THE FISH POND PROBLEM 



. "Denis hate.d t£ burden his wife with further bad 
news;" Arclet said, "but hated more to keep anything^ rom 
her, so he told Angelica the events leading jup to his own 
faring. However, af ter^publicly apologizing to Al for . 
his brash words, Denis had been told that ' he .could have 
his job back if, he brought ii>, by Monday morning, a com- 
putation of the correct lens thickness. That* was fine,, 
but even though he knew that Al's answer coul* not be 
right, he, himself, -was not sure. how o to go about finding' 
the correct* answer. 

"'Tears ^formed in Angelica's eyes but she brushed 
them away. 'Well, ajt least we have plenty to eat s For 
a while, anyway. 1 - 

"She was referring to the fish. * 




✓-"Ten'days previous, old Cy Seepage had put a dam 
across a 'stream running through his property in order 
to trap some trout. Incensed by the blocking of their 
favorite trout stream, Denis and Mike Mossy had, the 
following day, \opened the drain. Nater flowed from the cJ 
pond, carrying with it trout. Denis and Mike, collected 
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th'em in a wire basket hidden by the dense, hillside 
brush. They had left the basket in place day and night, 
occasionally 'emptying it and storing the trout in Mike's. ' 
freezer. Every day, old Seepage* checked the level. of 
the>pond,' planning to net a ll the troufas soon as the 
pond was full. The level continued to rise, for water 
flowed in faster thair.it could be drained out. Seepage * 
had no reason to suspect foul play." 

"The Dropmore's telephone rang. 

"It's for you'.' Angelica handed the receiver to Denis, 
'It's Mike Mossy. \ She wrinkled her nose. 

# 

"'Hi, Mike, what's up?' 

- J "'The pond is. Nearly full. Seepage says he's going 
up tomorrow to catch and count the fish', and half the town 
is going. up to watch.* ■ • 

V - • / 

"• That's great, Mike. ,A11 have to do is divvy up 
the fish we swiped.' " •» * 

■ *- - : " ' , 

"'Er-thafs whj/I called. I've got a proposition. 
Let's each' of us give a guess at the- fish left in the 
pond. The one coming the closest gets all the fish. 
How about, it:?' . ' 

"Denis was hesitant. Mike wasn't too bright, but he" 
war lucky; always winning contests,' Denis, was finally 
shamed into agreeing, ' - 

"'How could you be so stuprn, • Angelica cried, when 
she heard what Denis 'had done.' 'Don't you remember that 
both our families are, coming to dinner tomorrow night, 
and I was going to serve' that trout? .Ten guests coming * 
and no main dish, and.no money left either. What are we 
going to do?'" % 

** 



Chapter 4 

'MODELLING THE OPTICAL FILTER PROBLEM 



By now, Herb and Polly and I were all involved with 
Denis and. Angelica Dropmore. We begged Arclet to continue 
the story even after the bell rang. He just shook his . 
shaggy head, so off we went, discussing' their problems 
continuously until class the next day. 1 ' 

Arclet came into the classroom with a cardboard carton 
and, dumped the- contents onto the desk. It was a funny col- 
lection of junk'. There was a desk lamp, measuring spoons, 
wood blocks, a bag of peas, a hacksaw blade, and lots of 
other things.* Arclet stood beside the heap, looking arch. 

"Aren't you going to continue the story today?" we 
wanted to know. 

- "This is the continuation of the story," he said 
mysteriously, indicating the pile of junk with a "grand 
sweep of his drm. He was a pain 'when he put on the 
theatrics. 

He cleared his throat. "To continue, we left Denis 
on a Friday night, with a highly troubled mind. At eight 
o'clock, he received another 'phone call. J It' was Mr. ' 
Usury, the president of the bank. He was very kind, and 
said he suspected that what Angelica said about the beam ' 
*■ was correct. If she could come in by 8 AM on Monday with 
> something- that would convihce him that she was right, Mrl 
Usury would not only rehire her, but wbuld give her a raise 
as well, 

"The challenge was now clear. Angelica and Denis 
were faced* with three problems, one of which had to be 
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solved by three' o'clock* the following day, and the other 
two by Monday morning. . * ■ • N , 

"When Angelica came down Saturday morning^ she found 
Denis asleep at the kitchen table. Heaped on file table 
and the counters 'were ■ the very items you see before you 
now, and* sheets of papeV covered with* scrawled figures 
and graphs JL ay strewn about the floor. He -had devised 
ways $o model the tffcree problems. * 

"Is that r what you're going to £how us now?" I asked 
Professor Arclet. v , • 

«. Vtiot exactly," 'Arclet said. "That is what you're 
going io show me." ' 

Not' only did Artlet make me perform aij experiment 
right-on the sfot, but he make me write it up as* well, 
complete with .objective, steps, conclusion- -the whdle 
bit. Much of the»following description is lifted right 
out of my notebook. 

« 

Title : THE OPTICAL FILTER EXPERIMENT 

, Objective: & learn about filters. 

» 

At this point? Arclet objected *to my objective as 
to^rarfowT" I crossed it out and wrote: 




Objective: To see hop 'different equations can 
arise^from a physical problem. % 

x ' m ■ : , ,-< 

•Now Arclet objected to my .objective as being too 
objective. * / \ 

"Make your objective more sub j active," he urged. 

"What I. mean," he said in response tosmy vacant 
look,, "is that .you should phrase, your Objective in terms 
of what;, you will be able' to da after the lesson, that 
you couldn't do before." V \ 



So I wrote: 



Objective : After this lesson,* I^should be able to 

write a differential equation to describe 
a physical problem, (if IVs real easy.) 

Arclet sniffed at this, but let me go on to the 
step. ■ 



Materials: 



1. Desk lamp with 60 watt bulb. 

2. Photographic exposure meter wltb a 
scale graduated In EV (exposure value). 

3. About AO sheets of translucent tracing 
paper, or typist's onion-skin. These 
will be our filters. 

A. A sheet of glass large enough to cover 
the front on the, lamp Tiou*lng. 



Procedure: 



Step J: 'Swivel the lamp upward, so t$at It points 
at the celling. Place the glass over Ft,' 
and tape It In place. PJace one sheet of 
paper (filter) on the glass. 




gur* 1: Apparatus for^the Optic*! -Experiment 
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Stec 2: Hold-lng the meter right againSt the filter 
c take a first reading. This will be the 
starting value,, so it should produce a 
needle deflection close to the top »>f the 
, scale. If it's too high, cut d^wrfCfhe 

light by using more filters or a smaller 
bulb, (f too low, use a larger bulb; 

Step 3: Once you get this large deflection, .tape 
down the filters and mark the outline of 
- the meter^on the -top one so the meter 
can always be returned to the same loca- 
tion. Rotate the meter dial untilj*the 
number ten of the EV^scale is opposite 
the needle. — • • 

Turn off the lamp when not^aking readings 
so the paper won't overheat. » 

Step k: Add filters, singly or in groups, each 
^ time noting the meter reading and the 
* total numberof "filters added-. Do not 
coun^ the filters taped to the lamp, as 
these were only used to adjust the start- 
i ing Intensity, and to control the color of 
the light reaching the meter. 

# Step 5* Keep adding filters until the light reach- 
iog the meter Is so' weak that no further 
readings can be obtained* 



I wound up with a two-column table of data look- 
ing' something like this:* 



^ Experimental Data : 



x 



t 

Number of 
Filters 


* / L 

light level 
; (EV) 


0 


10. *0 


1* 


7.5 ■ 


2 


5.5 

* 
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Plot of the Data: 



10. 

Ui 









■ i » • 
U-L- 


i ; — 






- 








. 1 i ! 


— i 












1 ' "■ 

1 — i — , — 


! 
i 










- 
















— i 


i — 




— 


— - 












i i — 

: i 




















i : i . 
















; 
















i — 










1 


— 
















H 


-i 






I 


























1 


i— 








i 




























































k 














1 — r— 
1 

































5 10 IS 

A/aMBER OF FILTERS 



Figure 2: Plot of Light Transmission w . 

Vs. Number of Filters JF 



After plotting the graph, I attempted to explain to 
the others who were looking over my shoulder-. 

"Now, it's pretty ^clear that not every fijter is / 
absorbing the same amount of light, for if it was, the 
curve would ibe a 'Straight line sloping downward. ^It's 
plain that each filter is absorbing less light than the 
one before it, because the reduction- in 'light per filter 
gets smaller, for each Additional layer." 

"That makes sense," responded.* Arclet . "When we say 
that a filter absorbs,; say, 25 percent, it'means 25 
percent of the light Teaching the filter. Since less 
light teaches the farther filters, they absorb less. 

"In other words: The amount absorbed depends upon 
how much is present* 

. ,l In other wordsj The raze of decrease of^the light 
IpVfA is proportional to the light level. 
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< "In other words: The rate gf change of' L ie 
proportional to I 9 if we let L represent* tk& light 

level* n % % 

* * • 

# . "How can you talk about a rate?" asked'Herb. ."The 
filter, isn't moving." * , 

"Quite sp," said Arclet, "but do'we need movement 
to have a rate of change? What, about interest ratt.-ehf" 
"Oh, yeah," said Herb. * * 
Arclet went to the board and wrote, : 

[The rate of change of £] [is] [Proportional to] [£] 

"Now for an equation," he 'said. "There are four 
expre^sion^ each enclosed in* brackets, and four of us: 
Each of us will write the proper' mathematical symbol below 

. one of the expressions. I'll g0 first." Under [Is] he 
chalked in a huge equals sign, and, with sickening coyness,. 

; handed' the chalk to Herb., "Next," he said. 

Herb dashed to jthe board and, under U] , wrote* "£" 
and tossed the chalk to Polly* . 

Under [P.ropqrtional to]^Pg04y^wrote the proportionality 
♦syabol but just -stood rf^r* withja dissatisfied wrinkle- 
in her nose. "How can weUave^hjJtf/an equals sign and a 
proportionality sign in tif^ame equation?" Her quesl 
was directed at the blaa/bbard*, *ut Arclet answered instead? 



ind a / 
is t ion ( 



"You cannot have bofch. Instead of the proportionality 
symbol'; how about a 

"Of course," said Polly, and she wrote 'a constant of 
pxoportionality, k* The equation now read, ' 

t tThe rate of v change of hi *\kl ' 

- and the chal'k was passed to me. '« 

* • • ... . ' * 

'* 

15 
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Remembering that rate of change always meant a 
derivative, I wrote, , 1 ' v f, . ' 



dl v — 

Under the remaining expression. 

"Oh, I see," said Arclet mockingly, "The derivative 
of L. with respect to I-don't-know-whar; -Very, resourceful . 
Smarting, I wrote, 

' si*. " m i 

r 

where t.j total thickness of filters Ofuraber of filter 
sheets) . ^ 

" Just one thing missing now K " m said Arclet and,' after 

waiting a moment for one of us to supply :,£he missing ghing 
went to the bfcard and wrote a minus ^igqSuth a big 
flourish. '."The slopes are all. negative, right?" * 
The equation was then, ' 




.V 



To see if this equation really held for my data, I 
sketched in vt£e tan'gent to" the curve at several points. 
I then measured the slop # e of those tangents "(riW^ oyer . 
run). At each point where I measured the slope, I read 
tH£ ordinate L K and* divided the slope, (dL/dt) , by the 
ordinate £o obtain fc. * 



dh/di 
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1 got yifPproximateiy the sqjne. value for k at all of 
the ^iiifv taken. *. * 

' ^ "Not* bad, not bad at .all ,"' Arclet said, as I returned 

to my seat. "Now take another look at that, Equation* Is 

the.re anything different about it, or. i<s it like, all the 

other equations we've seen so far?" < 

« " 

"It's got. a derivative in it," Herb said. 

* "Right J" Arclet was* encouraged by the ale^t response, 
"And any equation^ containing a derivative is called.. 
He waved his arms like a conductor, trying to get^j£ 
fo sing put the remainder 0 f his, sentence. We sat 
there silent. . • 

,f A DIFFERENTIAL EQUATION, M he hissed through his 
teeth. ;,, Class dismissed. tr 



>J T*»JillllL- ' 
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* Chapter 5 ' * 

MODELLING THE SAGGING BEAM 

* "Tcdav, we will do another experiment,; 1 Arclet 
said after we'd settled irtto our seats. "To 'begin/ 1 

"But the story," we hoVlec^ \ - 

• "Are you here f or^learning u or fofr entertainment?" 
the^rofeSsor growled. (We knew we toad g%ne too far/) 

^len, more softlys "WeJ 11 get back to the story. But 

first it J s important for „ you to, understand the experi-' 

ments that Denis Dropmore stayed- up all night*to perfoVm 
". p» y \ 

in his kitchen. , 

"Now, who will do this one? Polly, my dear, Come * 
up-here." . * " ^ * 

Polly went up to tfie deslc^and performed* the experi- 
n\ent with some assistance From Arclet. The following 
outline is -almost word-for-word from the no£es I took 
that day-; , i ' ' M . " ^ A 

■ ' "\ A • - > * k 

Titles TH£ BEAM EXPERIMENT • 

. Objective : To* be able to make a model of a beam deflection' 

— *0 9% 

prob%m. 

Materials : 1 . Hacksaw blade » 

2. Sheet^f wood", at* least 8i*by^11 inches, 
* • 'and at least i Inch thick k ' *• 

♦ ^ 3- Two nal Is o 

* * 4. Sheet 'of rectangular graph pager 

5.* Ruler , 



Procedure: * 
Step 1: 

Step 2:' 
Step 3: 



Step 4: 
\ 

Step 5: 



Step 6: 



Tape the graph paper to the board, so that 
the edge of the paper Is even with, the edge 
of the board. 

Draw a line down the 'middle of the paper, 
parallel to the long edge. 

' Drive 'the two nails into this mid-line. 
Let them extend about % # lncb* The distance 
between the nails should*. be '10 inches. 

Prop the boaVd into a vertical position, 
lay the hacksaw blade across the nails. 

One and a half Inches from one of the nails, 
pushfigainst the blade wfth the eraser end 
of af pencil. Get an assistant to* trace the 
x shape of the blade on the graph paper. 

Changing the pressure against the blade, 
draw several such deflection cuVves. 




\ 



Figure 3: Apparatus for the Beam Deflection. Experiment 
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Arclet pointed out that many physical problems required 
special knowledge of the field in order to even write the 
equation, and that this was one such case.' 

"Fortunately, I happen ^o have a broad engineering 
background v in addition to my mathematical training ,", said 
Arclet. "Now write this down. 11 He tilted his head back 
and defocussed his eyes, as if he were reading something 
off the inside of his skull. "The second derivative of 
tKe vertical displacement"!/ with respect to the position x 
along the beam is directly proportional to the bending 
moment M and inversely proportional to the product of the 
modulus of elasticity E of the material' and the moment of 
inertia I of the beam cross-section." 

This entire sentence was delolyjered without taking a 
breath .and in a fl^at monotone, like a third-grader reciting 
a poem on Parent Vs Day. \ 

As with the optica 7 ^ f ilteVproblem, we wrote the 
expression in brackets and,, with less difficulty this time, 
wrote the symbols beneath* 



The second derivative 
of y with respect 
to x 



is directly 
proportional 
to 



Y 



or, 



dx 2 





M 




J 






M 


0 



and inversely' 
proportional 
to FI 

1 

EI • 



d 2 y m 



dx 2 



EI 



Arclef then pointed ovk that E was, in fact, the constat 
of proportionality in this equation, so that a separate k 



Our 


final e 




M 


dx 2 


EI 
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.* . 'fpon't woriV abqut this equation now," Arclet said. 
"We'll go Jhto it in' detail when I explain later how we 
solved' this problem numerically. For now, it's enough 
thaf$ you ?find the point of maximum deflection on your , 
hacjcsaw curve^and verify that it lies between the load 

O 4 « 

and fhe midpoint, as Angelica had claimed. M 



O 




1TTT Curves $co*\ Bww\ De4lecf7ov\ tr^pe^^ev^t 




'^te:|.;vl:r:,J. ^ |>f ::l.::.|r: 1: 



- V The, curves -Polly drew of the deflected ,hacksaw blade*. 



' ; , - £V*l*f e . n ^ Qtir second lesson in writing differ- 

, \ .e&t^tt; equations ^f or different, physical problems,, Arclet 
V- g^ve 1 us the 'following qydz.. 



G 




Quiz #1 

1. If the rate at which heat is lost by a body Is propor- 
tional to the temperature difference between the body 
and its surroundings, write a differential equation 

to describe this situation. 

2. An object moves In a medium offering resistance * 
proportional to the square root of Its velocity at any * 
instant. Describe this situation with a differential 
equation. 

3. State in words *what a differential equation is., 



r or answers see page 67. 
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« Chapter 6 

* PROFESSOR ARCLET TO TrfE RESCUE 

Today, Arclet was seated on the edge of his desk, 
a signal that the story was to resume. We took our seats 
in a hurry. "^Ny j t 

"You will recall, "^he beg^i, "that we left Denis 
Dropmore asleep in a kitchen filled with the debris of 
the two experiments we have just performed. Sweet 
Angelica, fearing that her husband^s mind might have 
s*napped under the stress of their troubles, woke him 
gently/ 

"•Careful of the fish! 1 he shouted, pushing her arm 
away from the little pile of green, split peas with which 
he had unsuccessully tried to simulate the fish pond 
problem. # * 

"•There, there, Denis dear. Of course I won 1 t hurt 
your little fi^hies. Now, you come upstairs and lie down 
whiLe I make some 'phone 'calls . 1 

"•No time, no time, 1 he cried, and began rushing 
about the kitchen, rumpled and unshaveii, counting piles 
of split peas, an^ holding sheets of tracing paper up to 
his swollen, red eyes. 

•"•So much to do, arid so little time/ I feel that 
I'm so close to an answer but X 4pn f t know^what to do 
next. 1 „ - ^ 

,M Now, Denis, piease lie down, .and let me call the 
doctor. • * . ' 

- ,„ % "Denis stopped dead in his tracks.^;* ^* 



,u The doctor? The doctor! His eyes were wild, and 
he shook tge hacksaw blade menacingly, f 0f course! I 
need ttje Joetor. Angelica-,— you're a genius- 1 He ran to 
the phone and called the only doctor'of mathematics he 
knew, -Me, Arclet, of course, Realizing the urgency of 
his problems, I went over at once, 

"At Denis 1 house, we exchanged the briefest of 
greetings, and went right to work, I scanned the results 
of his experiments, and v^as astounded at their cleverness. 
Realize, of course, that he had- no equations written, as 
we have done, and I set about yriting these, I sent Denis 
off to take some measurements at the pond, 

"He returned in an hour with the measurements and 
estimates, and with the news that people had already, 
■begun to gather at the pond, Mik* Mossy was there, 
annoyed that Denis didn't have his guess ready yet, 

"Since the fish problem required a solution sooner 
than the other two, we tackled that first. We solved * 
it simply by drawing a tangent field* 

"Class dismissed." 



Chapfer 7 

TANGENT FIELDS, AND SOLUTIONS TO DE'S 

Arclet was so obvious. He had left us hanging yester- 
day, with that' reference to tangent' fields : today we were 
suposed s to rush into class and yell, what is a tangent * ' 
field? What is a, solution? Teach us", teaches. I hate 
to admit it, but I was a littlfc curious. Arclet came intp 
'the room. 

* *, 

"Lesson today, right? 1 ' we asked. He winked and 
began writing on the blackboard. 

Title: THE MEANING OF A SOLUTION TO A DIFFERENTIAL 
EQUATION. INTRODUCTION TO TANGENT FIELDS. 

Objective: &y the end of this leaaon, you should be able 
to:' " ' j • 

a> Verify -whether a particular equation^ graph, 
.or tapje of point pairs is a solution to a 
given 'differential equation, 
b) Draw a tangent fieV&l 

a) Use a tangent field to sketch a solution. • 

^ These preliminaries out of the way, Arclet began 
*to lecture. * 

"When you solve algebraic equations, what do you 
get, aside from a headache?" 

It looked as if we were in for one of Arclet f s 
lighter lectures, 

4 4 

VYouf get s.ome number, the root, which is the value 
of x at which the plot of the equation crosses the i-axis. 
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quadratic 



v * the/cutV;e/cros ses in .more than one place, as a 

? , v way, you get: a.;roQt; at each crossing. Fine. 

r V^^ : &>? -'^But^h^ differential equa- 

HU;4tiont . What sbrt>OL§, an^r;. can we expect? - A^iimler? 

One objective of this, unit 
■fr ,T to help^you tq ^understand What the" solution to a 
"i 4 : ] .differential equation ;means> and how:io recognize one, 

,;Af ter) coapieting this, unit t you should be able to examine . 
. what- i.s% claimed. ,t o" be the solution of a ^particular DE,.and 
%?X$ty impostor or not. We will 

* S°£r^P$ <s 5 ol ^i0ns 4 V Not- >y e t * ^ , 
f > "Look it this differential equation: 

Y 



Jfwe think in graphical terms, dy/da: is the.slppe of some 
c - curve.-y « rWe«,canrr^ m. For 

^ '-siipiipity; . \t rewrit^ tlie DE> using w. • 

sr- t^'-' m *t2y -x' 0 ' '? \ '' "' I 

: V"- 1 % * v-"- v ' - -1 .'„ , r , , '* t 

1*1^1' P*f*f ^"^^i J**/* 4 ^"'* £h at , will^satisfy this 

Take, for Example, m ■'l f \ 
V^Vi^ ^S^ahi^y^i. * Substituting/" we 'get^ i - .„ ; 



r . i f 2(»> > s - o . 

* "It should be° apparent that there" are infinitely jnany . 
such combinations, that will 'satisfy our equation* Just- 
lO ^:r/PWk t «iy valiies for * and y out of the air, and solve 
;*^J tW equation for mi \. « *- , ? %t \ 4 * 

" ~" "Suppose , that we pick the vajutfs 1 and 5 for x and y* 
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Theft 



* - 2y 
1 . 2(3) 
-5 



Thus the set of numbers -5, 1 and 3 will satisfy o<ur ] 
equation. Now cdmpute three qx four. more sets, of numbers 
that will also work. To keep from getting too spread 

out, take values of x and y between 0 and 10. 

' * ,* 

"Got the numbers? All right, now plot them on a 
sheet of graph paper, using x for the abscissa and y for 
the ordinate in theN, 4 usual way. Alongside each point, write' 
in the-vaiue of m. You shoul4 now have something looking 
more or less .like this. 
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gure 5: 



.S\ /ft , 

Slopes at Various Points 



"Draw a smooth curve through the points.. In fact, , 
try to draw in your curve so that it will have the slope 
required at. each point." * 

"It doesn't work, right? Don't woyry, I didn't 
expect it to. But now I ask you the following question. 
Instead of^chOosing x and y at random, can you think of , . 
some clever way of picking them so that when all the points 
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are- connected, the slopes at each pointjare precisely 
>hose, required by the equation^ Why should you go \o all 
this troubled Because ' * 

"WHEN YOU PO THAT, YOUR COLLECTION , OF ( JOINTS 

WILL BE A SOLUflON.OF THE DIFFERENTIAL * x 

EQUATION. » 



I wouldn't asjc you to do something that wasn't important. 

"At this point, take five minutes to struggle with 
\hat problem.' You can work by trial and error if yoji like. 
A good way to start would be to draw a shojrt line thro ugh 
each^of your points, with a slope of the proper value. 



V 




f igure 6: Part of the tjangent^f t«1d for % + Zy - x - 0 



'After five minutes, Arclet t>n>te on the b^ard, Tangent 

Fields. He then continued" to lecture. 
/ 

"If you had graphically indicated the slopes at your 
points, and had adcied more points your plJt, you would 
arrive at a diagram knoun as a tangent fi eld , or^a direction 
field* IT we sketch in a ci^rye whose clc~e is always 
Ihe same, as that of the surrounding tann.ent field, we 
know that any point on the curve wi^l satisfy the differ- 
ential equation. The 'collection of points on the 0 curve 
is therefore one solution of the DE. 

"Now g<> ahead and complete your own tangent 
field, and sketch in several possible solutions." 

- He gave us a few' minutes to do that before * 
continuing. 

Different Forms o5 a Solution 

"The curves that you sketched previously are 
solutions to the differential* equation . For the*, 
"moment, don*t worry that there are an infinite % 
number of such curves. We get to that problem . • 

.later-. The information in these ^cjurves can be . * * <4L 
presented in three dif f erent Vays . . . * 

a) Graphically, as they are now. * * x 

b) Numerically, as a talkie, of x>y pairs. 

c) Analytically, ^as the equation .of the * 

4 curve 'connecting >he points. < < J *• - # « 

» * - - V 4 " ' ** ' 

"In^our ne'xt lesson w£ will solve a^ DE graphically, 
and obtain a curve. Later, tile 1 11 do a numerical * 
solution and arrive at a, table of point pairs for 
an answer. An analytic solution is obtaine3 A by « „ 

manipulating the differential equation. We 1 re not 
'going to do that." ' " • / 
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;'X I raised ay hands to clap, but Arclet's icy stare ♦ 
stopped me cold. " - , 

Boundary «: Conditions J 6 

; "Why so. many solutions? At this stage in your 
. math career you should not be too surprised to get an 
infinite "number of possible* answer^ when you reverse 
a mathematical operation. Iftfck what happened wljen you * 
tooJc the inverse of a trig function; or an anti-derivative 
Involving a diffgrei^fial equation, integration must 
be performed, and it. is the unknown constant of integra-'v 
tion which makes the result indefinite. As with the — 
-anti -derivative, some additional information is needed 
before the unknown constants can be evaluated. These 
additional facts are called boundary conditions <*" 

"For example, supposel told, you that, in our 
previous problem ,^tfie slope had to be -7 when a: .was 
zero. Wouldn 1 1 jjp merely^draw in the one curve having 
the required slope whereit crossed ^the. #-axis, and 
ignore Jthe rest? Or, if you Khew that the curve had 
to pass through the point- (4, 5), this) information 
would exclude all possible curves bux one. 

/ ' X 

"If you, have access to a computer,", you might 
enjoy writing s program to plot tangent ^elds. . Have 
the program print, numbers giving, the slope, at the , 
proper coordinates N on the paper. If» yqu'hav,e a plotter * 
available, ^ou even can haVe it .draw in the . tangent 
^lines. Here is a tangent field produced by such a . 
program ' . 



U4> 



Figure 7. A Computer-Generated Tangept Field 
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bNe t where k - 0.5, an'd b 



0.5. 



\ 

Arclet* concluded the lesson \by giving us the 
following quiz. » • . • ' 
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Chapter 8 



Quiz #2 



1. Verify that the equation 

y - j (In x + 1) ' 
Is a solution to the differential equation 

by plotting It for values of a: from 1 to 5, and measur- 
ing slopes of at least three points. 

2 t Given the differential equation 

plot the tangent field, letting x range from 1 to '5 * 
and y range from 0 to 5. Sketch In a possible 
Oblution. \* ' 



,«4 



Por^te 



nswers see page 68- 



THE FISH POND PROBLEM SOLVED WITH *A TANGENT FIELD 

"Time was of the essence," Arclet said, continuing 
the tale. "Mike Mossy had just called, 'reminding Denis 
of their agreement; if one of thera-failetf to make a 
guess by the^tirae the fish were counted, he forfeits his 
share automatically, I worked furiously over the graph 
while D^nis.fed roe the measurements he had taken 'and 
Angelica did the computations on a pocket calculator. 

"Here what Denis found at the pond. 

. a) Water flowed into the pond at an estimated rate 
of three gallons per minute. 

Each ten gallons flowing in brought with it two 
fish. 4 ' * 

c) .^Greedy old Cy; Seepage had gotteji the Fish *and 
Game Department ^o stock the pond free of 
charge by telling them he was going to allow 
* public fishing; They had put in one thousand 
"trout the same day tha*t the dam was finished. 
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Figure 8, Ffeh Pond Flow Rates 
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d) Since Denis and Mike had opened the drain after 
one day of fillirtg,' the pond would then have 
confined 4320 gallons. 

e) Water ran out 'the drain at the rate of 2 gal/min. 

f) The pond was completely filled nine days after 
■ the .opening of the drain, at which time it 

contained 17,280 gallons. 

"Denis an3\Angfil±ca-. watched with* awe and gratitude, 
as I developed my clever solution. £ 

"Let Q represent the number of fish in the pond at 
any instant, ^hen,, 

dQ 

ZTF * ra * e o£ chan 8 e of the number of fish 

Also, let 

C * the concentration o£ fish; 

that/Js f the number] of fisij per gallon of water, at any 
instant. We assume that the fish are evenly distributed 
^tl^roughout the pond. Since water containing 0.2 fish per 



gal^ 



on was .entering at the rate of three gallons' per 
minute, fish were coming in at the rate of 3 x 0.2 =0.6 



|h per minute. ♦ / 



"Since water containing C fish per gallon°was leaving 
at the rate of two* gallons per minute, fish were being 
drained off at the rate of %C fish per minute* Here we 
are making 'the assumption that the drained water contains 
the *ame concentration of fish as in any other place in 
the pond. 
» . 

,"The rate of change of the number of fish in the pond 

is obvious!^ [fish in] minus [fish out] or, 

. §jr « 0.6 - 2C. 0 (equation 1) 
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The total amount of water, W, in the pond, if we start 
reckoning time at the opening of the drain, is, 

W » 4320 + Zt - It ■ 4320 + t 9 

where t is* in minutes. * . 

- "The concentration of fish is, then the total number 
of fish divided by the total amount of water, or 

* ' _ Q Q 

0 " V 9 4320 + i . (equation 2) 

By substituting equation 2 into equation 1, the rate of 
change of the number of fish becomes, 

Y 



dQ 



0.6 - 2C = 0.6 



,2Q 
4320 + t 9 



which was the differential equation describing the fish 
population in thg p^ond. 
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Figure 9a: The Start of Arclet's Graph 

"The next step was to plot a tangent field. I 'let 
thr! abiscissa be tine, in minutes, extending from zero 
to * 

9 days x 24 hr/day x 60 min/hr « 12,960 min / 



Ivrould also compute the number of fish in the pond at 
time. zero. It was, 

1000 '+ 4320 gal x 0 . 2* f*ish/gal = 1864 fish. 

* *' * ' 

(See'figure 9a.) - * « ^ 

"I next computed^slopes' at' various points on*\he 
graph, naturally taking nhem only where I expected the 
curve to^Be^ *I computed them from the -differential^* 
equation \ * »• 

slope = °' 6 - TrnF T 1 ! ■ 
No vr this was the graph: 
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" Figure 9b: Arclet's Tangent field for the*FIsh Problem 

and I could easily sketch Ufa the curve showing the growth 
..of the -fish populati6n. Where the curve intersected, the 
12,96D minute' fine was J the number of .fish at the hme the 
pond was full. ^" — * • * 
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Figure 9c: Arcle-t's Solution to the Fish Problem 



"Denis startled me by tearing the sheet from my 
hand before the ink was dry, and flew out the door. 
Angelica and I waited for several hours. Their dinner 
guests started to arrive, and we maxle feeble excuses 
for hi? absence! He returned home at six. We had a 
very fine dinner of boiled potatoes, broccoli with 
hollandaise sauce, pink chablis, and, oh yes, broiled 
trout A 

Arclet stopped speaking. * His eyes were^lazed 
as he smdcked his lips and recalled that fabulous 
dinner. Herb wanted to drop his calculus book again, 
but'I said no. It frightened me whenArclet got like 
this. • - 

We slipped out after the bell rang, leaving the 
professor sitting, on the edge of the desk, muttering 
something about 'lir^ sherbet. 
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• « Chapter 9 

SOLVING DIFFERENTIAL EQUATIONS GRAPHICALLY 

v • * <* • 

Today Arclet was a^l* business . He strode briskly 
into the room, always a bad -sign, and began *a rapid-fire 
barrage of questions ♦ What was wrong with his solution 
of the fi sh problem? Why draw infill those tangents 
when only a few are needed? How. could you limit your 
work to only those needed? And so on, ' 

I slumped in my chair, \ It was a good day to 'lie low! 

Confronted with a cpunter-barrage of silence, Ajxlet 
sullenly began the lecture -he had been softening us'up 
for* ije wrote the board, J t * 

Title s GRAPHICAL SOLUTION OF DIFFERENTIAL EQUATIONS - 

Objective : Upon completion 6f tliia Jleeaon, you should be 

able to solve a differential equation by a 

graphical application of Eudev's Method* . • 
" N • * 

\ In ^ 1 * a * UIiit Wj 5 solved a differential equation 

graphically {^plotting a t tanggnt fiefd," Arclet be^aii. ' 

"This is fine, but it takes "more work than 4s reaJ,lyT* 

needed,^ Jy the €nd of this unit you' should be afae *to - 

produce a graphical solution by thef fgich easiej Euler's" 

method* ^ ' * 

" ' % • - 

"I- will illustrated the method-'by" doing -an example. 

We wish to solve the equation* 



with the boundary condition that y = 1 when 1. Suppose, 
also, that we're interested only in the region between ^ 
x « 1 and x * 10. 

"We' start by plotting our starting vajue, (1,1) on 
rectangular Coordinate paper.. Now our differential 
equation tells us what the slope at any po^nt should be:. 



• ( slope . 

t y 

At (1,1) the slope. should be 

m*= 12 - 1 

m * y ■ 1 , 

m so, through (1,1), we draw a short line with a Vjope of 1. 
Our construction should be: 

\ 




dx Y>> 



Figure 10a: Construction for a Graphical Solution * 



"Now we step along the x-axis to a new. value df. x_ 
The. size of the* step is important; large steps give inac^ * 
curate results and small steps are a lot* of- work. I'm. 
going to take steps'of'l forlfrhis ^monstfttibn. 

"At x *-2, the value of y as read from our previously 
drawn line, is 2. Mark this point on the graph. Through 1 
• this point, we draw another short- line having the required 
slope at this point, which is: 



^ F I gure-40fchr— Con tfn u a t-ion-of-thtf-eonst r uc tt oiy- 



At'x » 3, y has the Value of 4, so the slope of the next 
segment *is* . 

'•■t: 2 - 25 - 

"The process is repeated as" far out as is needed, in 

this case until x = 10." ' t 

> • , - *■ % 

At this-point Arclet- paused dramatically. . Then he* 
flipped on his ever-ready overhead projector, pulled out 
a transparency, and said triumphantly, "Viola, ^the- final 
solution!" (see figure 11) 




Quiz #3 

1. Given the differential equation 

ax Zx * 

>fith the boundary conditions, y « 1 when x « 1, solve 
graphically by Eulefr's method, taking intervals, delta 
x of 1," and find the value of y when a? equals 5. 



answer see page 70., 



Chapter 10 . " * . 

/ v 

A GRAPHICAL SOLUTION OF THE FILTER PROBLEM \ 

* *> 

A • ^Wasn't it strange how-Arclet's story was requiring so 
much work from us to be able to, understand it Befoue he 
could begin nil' narrative the following day, I nad 'my harfd 
up * 



"Professor, i*s this a true story? 11 " ~. v ' 

t 

s ey< 



^Arclet looked over my head, focussing his eyes on 
infinity, and sighed. "My boy, what is«£ruth?" £| 

v * * * /I 

"Heavy," breathed Polly. . /• 

> ♦ ' **» 

"Bull," muttered Herb. * ^ ' _ ' » 

"To continue,"* said the pjrofe*ssor, "Denis v ^ngelica, 

and I were groggy from/fche huge meal, the cha#3^s ? and the % 

cognac J had thoughtfully remembered \o brj.ng,*.so we'didn't 

tackle the remaining^ two problems after tfteguests left 

but waited until the next morning, 

• - « 

"We decided tfte**next most urgent. pfobTem was* to deter- 

♦ . • * ?■ 

mime" tjie* thickness _of the sunglass lens, 

"Angelica had persuaded .Mr . Usury to let^Ber into the 
bank, so while Denis and L worked'on* the lens problem, 
she went off to take measurements on the steel beam* 

"We had decide^'^o <fo 0 a graphical solution. On a 
sheet of rectangular paper, we took lens thicknessjt^in^ 
millimeters along the •abscissa, and light ^transmission L 
along the ordinate, graduated from zero to'100*. Wej 
marked out two known points : $ transmission # 6f 100 percent 
cat a thickness of zero, *and a transmission of 85 percent 




si**' ?* th ? c Mi B ?? of lam* Vie connected the'tso points,- and* 
<>**tff. ; *> <: o.easure,d; pie, slope frf the line.; Our gVaph was then: 

... • -V 
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4 6 . t » 12. 
Figure 12: Start of the Graph! ca fool ut Ion 

"If you recall from our filter experiment; the 
differential, equation We arrived at was: 



slope > , 



where I is the light transmitted, t is the filter thick- r 
ness, and k is a constant. We had to find k before we • 
could- proceed with the construction. We solved DE for k. 

* * ».-»*** 



slope 
£ » 



':>;. , at the fiVsit' point on our graph. 



V 



Here. ; = 100 and the* slone is '-15, so - 

-15 



k = 



= 0.15 



and we used that 'number throughout' -the construction. 

"Through our second point (1,85*) we drew a line with 
a s.lope of; [ • * . 

. , Slope = -0.15(85,) =■ -12.75 , 

and saw. that it intersected the line t = 2' at a transmis- 
sion of 72.25* W£ Icept -repeating the computations" and 
the construction until our." curve dropped below the ten ' 
percent transmission level. ' Here is the complete graph. 
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. <^ s '4 4 6 8 Itf 12 14^ |6 
Figure 13. Graphical Solution of 4fe* -Ja, for * « 0.15 
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\. - *'As. $ou can. see, the lens thickness required was 

"Denis could, now prove to his boss Al that the 
required lens thickness was. not 6mm, save his job, save 
the eyesight of our fighting men, and keep his family 
together « * * 

"At that moment Angelica came home with measurements 
of the st.eel beam. Her display of gratitude and affec- 
tion more than repaid me for my efforts." 1 

• "We'll get'hi* for this," I whispered to Herb. 



• ^Untt 83: SOLVING DIFFERENTIAL EQUATIONS ^NUMERICALLY 

Chapter 11 Solving Differential Equations Numerically . . . . 

Quiz M * 

-Chapter 12^ _ A JumeJI c_aJ__Solu.t Loh„to .tbe^agging_fleam_P_rob lem . 
"Chapter 13 That Exam Again . . ~ 

That Exam 



• I x " Chapter 11 

SOLVING DIFFERENTIAL EQUATIONS NUMERICALLY 

. . 5 

"Yesterday you expressed some doubts about the 
veracity of my little story/' Arclet began with some 
sternness. "Well, this should settle the matter." He 
tossed the •phone book Tnto ay lap. f*opened.it where 
the-page- had-be^n-turned downy-andgazed af^the name— — r- 
circled in red. W ' 

Denis Dropmore, it said. ArClet would live to 
regret showing me 'that 'phone book, 

H I would really like to continue ' the'story," Arclet 
said, "buV before you can understand what happened next, 
it will be necessary first to have a lesson on the * 
numerical solution of differential equations," 

I sailed. It was all so clear now. 

Arclet walked to the board and began to write.- 

Title:. NUMERICAL SOLUTION ^F DIFFERENTIAL- EQUATIONS* . 

Objective: Upon completion of this lesson, ygu should 
be able to carry out a numerical solution 
of a differential equation by Euler's Method, 
by hand/ on a calculator, or by computer, 

•Then Arclet began;the lecture. /< 

"If you've understood the graphical solution of a 
ififferenfial equation by Euler^s method r*iSe~way to per- 
form a numerical computation should be obvious." Ke will 
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use the same equation we had for the graphical solution 
as an example: 

% - f- sl °p e m 

with the starting value, y * 1, when x » 1, We .now 

increment x by some cTiosen amount Ax.* Let's tate 1 for 

our step size now, and later we'll see the effect of 

making it smaller, m < 4 

* » 
"The slope of the curve, at our first point is, 

slope = — = i- * 1 _ 



and the change in ordinate going from the first to the 

* * 

second point will be the slope times the change in x, or 

J J 

w • Ly = mAxj/lU) - 1 

-* 

so the ordinate at the second* point will be 

'« 

v y+Ay«l + l=°2 s » 

■ i 

The second point' on our curve is thus* f2,2). 

"Repeating the computation, the slope at the second 
point is, • s ' „ 

* , x 2 2 2 * * % 

* slope = — = _ = 2 , 

K 3 2 • 
< 

The change iri ordinate is, 

: i/ = mbx = 2(1) = 2 . - v 

So the ordinate at ihe third point is 

* * ** ** 

* - y + *A# « 2 + '2 ? 4 . * 

The third'point is thus (3,4), * 

"We continue the computation in a* similar way as 
far* as we need to go, obtaining each ordinate y n frora 
the ordinate y n - x of the preceding point, and 'the slope 

0 ' 50 
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r » n *i ft the preceding point by . the equation: 

\ ,f A computer really cones an handy here.* The follow 
iag short program, written in BASIC,, is designed^to 
perform this computation, taking a step size Ax of 1. 

10 PRIflT "X"% "Y" , "SLOPE" 
Y 20 S»l 

30 Y»l * - ] . , 

35 FOR X-l TO*10*STEP S " * • . 
40 M-X+2/Y < 
50* PRINT X, Y, M 

6a'Y»M*S+Y , J 
* 70 NEXT X x 
^J&tiO END v fr- 

"This program. will also print the slopes, an 



Step size»1.0 



3 



extra dividend/ 


Here is a RUN: 






X 


Y 


SLOPE 




1 


1 


1 




2 


. 2 


2 




-3 


" 4 


2.% 


* 


<4 


6.25 s * 


2.56 




5 


8.31 


2,837 




6 


11.647. 


*3»090 




' 7 


14.738 


3.324 




8 


18.063 


3.543 




9 


21.606 


3.748 




- 10 


25.355 


3.94.3. 



"How that we have a program, it is art easy matter 
to determine -the effect of changing *the Vtep size. 
Notice that the program was written so that/only line 
2Q need b*e changed.^ Talcing a' step size of /0.5 results 
in the following table; * / 

:* ' X , .Y • SLOPE - - ' 

:1 -1 ' .• 1 . * 

• . 1.5 '1.5 ' i«5 c 

' 2, 2.25 • 1.777 Stepasize- 0.5 

• • • , 2.5 . 3.138 1.991 ' c 

... - ; ;„ . 3 , -4.134 2.17'6 . * 

- #*• * •*> ... 3.5. S.222» ■ 2.345 

"4 • 6. 395 2. SOI' , . c 

• * " 4.S, 7.646 2.648 * 



X 



5 


8.970 


2.786 


5.5 


10.364 


2.918 


6 


11.823 


3.044 


6/5 


13.345 


3.165 


7 ' 


14.9^8 


3.282 


7 S 


1 ft ^fiQ 




8 • 


18.267 


3.503- 


8.5 


20.018 


' 3.609 


9 


21.823 


3.711 


%% 


23.679 


v 3.811 


t 


2J.584 


3.908 




1 - 





4, 



"Now reduce the step size to 0.1, afc the same 
;timeTigging th^ program to print only at integral 
values of «, to save time and space. 

% • " . • J 

SLOPE * / ' 

hit ■ ■ ■ 

' 2.452 

* siooe step si ^o.i 

3.246 ' - 

. 3.470 - 
.* 3.680 . 
3.879 

reduce the step size by~ 





Y ' 




1 


2 


2.357 


3 


4\2S3 


4 . 


• 6.525 


5 


9.111 


6 


11.975 


7 


15.091 


8 * 


18.440 


9 


22-. 0G6 


10 


25.777 */ 






filial 


RUN,t>Jet^s 


■ of ma 


g^itiide.. 






X 


- Y 


1 


1 


2^ 


> 2.378 


3 


4.278 


4 


6 6.554 


5* 


; -9.143 


6 


12.-010 




•c 15.128 


A ' 


18:479 


9 


22.048 


10 


' 2$,.821 



SLOPE 
- 1 V 

\j.6ai 

* 2.103 
CT|»'441 „St£p size* 0.01 

t!997 
3,2-38 
, 3.463 
3.673 
3.872 



>ize= 0, 

7 



u,, At thi^ point, it might be ^interesting to compare 
'the final value of *y obtained u^ing the various step 
sizes. Let's make another 'table. 




t 


step size N 


ordinate at 
x - 10 


« 


- , 1.0 
, 0.5 v 
0.1 
0.01 . 


^F.355 
25.585 • 
25.777 
25.871 




-Theoretical * 
* 


25.8*26. 



"Not too bad. Even with our coarsest step, our 
final value is less than two percent different from" 
.the theoretically correct' answer, and with-our smallest 
$*ep§; our' answer is correct to four significant 
figures." < „ 

{ Polly interrupted. "Where did that theoretical 
value coae from?" 

"I was hoping you would be curious.* Remember 
when I spoke about analytic solutions? 3 ; Now would 
be a good time- for you to scan that chapter in your 
text, 'and we 1 1,1 get to it in a week or so." ' 

€ 

. "I'll do that tonight* professor," said , Polly. . 

Herb muttered ^something inaudible. 

At the end of the lesson, Arclet#gave the ' *\ 
following take -home quiz. ' * 



168 



( i 



f * 



if 



Quiz Ik 



1. Given the differential equation 

. dx x - 2x z y 
and the boundary condition 



y 88 3 when, a: = 1 

do a numerical^ solution by Euler's method, taking 
intervals no larger than 1, and find the value of 
y when x equals 5. . You may program this problem 
on the computer if you wish.. 



For answer see page 71. 
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• 4 

X n A NUMERICAL SOLUTION TO THE SAGGING BEAM PROBLEM 

Today Arclet began, "It was now Sunday night, and 
we had just finished a supper of 'pan fried trout served 
with % a savoxy made of soy sauce, chablis, and*a touch of 
dill.' We 'attacked the last and most difficult of the 
three' problems. 

."Angelica had 4°ne an admirable job, of collect} 
.the data for which I had seri^ her. We began by makinfe' 
«a diagram and carefully listing all we knew about the 
problem. 




•Figure f*». ' Free Body Diagram of the Beam 



"I drew in coordinate axes as shown, and indicated 

the forces acting*j>n the Seam." We felt -it reasonable 

to .assume that all other Toads , including" the weight- of 

the;be , am itself, would 'be negligible in comparison to 

the weight; of the huge safe. ^ ,* f * 
* .. '•' . ♦* 

"We found the dimensions of ttie beam cross-section 

in orie'of Angelica's architecture books, and ^t gave the 

foment of inertia. / * 

I - 12,000 inch* 



*We also found the modulus of 'elasticity of steel, 

JB * 30,000,000 psi 

"I took as my starting point the differential equa- 
tion for the deflection curve of a beam, which we wrote 
in class the otKer da^: 

lx* EI 

whereby is the vertical deflection at any point x along 
the beam, M is the bending moment acting on the beam, 
and £ and I are as defined above. 

. "Before solving the DE,, it was necessary to know 
how the* bending moment M varied with position x. I 
considered a section of Ueam lyi4qj t,o the left of the 
sate, : <^ \ , . 




and asked, 'What moment M would be Required to keep £hat 
section of beam dn equilibrium (l^ep it from rotating)?* 

"It wo^ild have to be J 

M * Rx . , • 

where R is the vertical'reaction .at the left-end of the 
• heam, and x is .the distance to % the section* being eonsidered. 
We can find R without too much trouble, "by asking what, 
/force R is needed to keep the, entire beam' frpm rotating . 
about the right hand support point. * 
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* Figure 16, Taking Moments About the Right End 

"For equilibrium, we must have 
30* = 6?, 

where P. is the weight of the safe. The bending moment is * 
then, 

M =■ Rx = | x , 
where * is in feet, so our differential equation 'is 

"Angelica had found the weight of the safe from the 
manufacturer's catalog, and we added some weight for 
contents, and took P » 10,000 pounds. Putting in the known 
values, I computed, 

. * . ( io. ooo ltn 

5Er 5(30,000,000 lb/in z ) (12,000- in") 



£ X 10 " y fl ^,^^■>)■ 

5(3 x 10')(1.2 xl0») in" " [TF x 10 J in 

and converting to dimensions in feet 

t*7'[h* H'^VPiT^ = 8 x 10 " 7 ft 

* So out DE becomes, v 

. •' ' \ 

, gf [«x 10' 7 ) x 
wi'th all distances expressed in- feet. 

i7o ' ; 
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"It *e now take the integral of both side's, we get 
g|= (8 x iO~ 7 J a£ + c -*~*4ope, 

an equation giving us the^slope at any s distance x along 

the beam. This is the differential equation we must solve. 

m 

"I de.cided to do a numerical solution by Euler's 
Method. For a starting value, I knew that y had* to be 
zero when x was zero. 

"At this point, fienis, who had followed my' prev^hts 
solutions with much care, interrupted me. 

"•How can you find - the slopes from that equation 
when you don't know the value of C?' 

!t, We have another boundary condition we haven. 1 1 
Used yet, 1 I replied. 'We know that the deflection y 
also has to be zero at\^e other support, where <x = * 30. 
All we have to do is keep guessing at C until* we find 
one that gives a displacement of .zer t o at the right end. 
We ought to get it in about fiFteen tries.* 

"•But don't you have to do the entire CQmjjutation, 
with each guess at C, before you know whether it is good 
or not? 1 Denis asked in apprehension . 1 . 

"That's right." . 

"•And how long will each computation* take? 1 
"About an. hour. ' * . ' 

"There was a long silence. Mt was now a lit^Le past 
midnight. 

"It was time to play my trump card. m I went to my 
car and returned with a portable teletype and £n acoustic 
•^coupler. A quick 'phone call to'th* computation center 
brought the full power of the college's massive computer 
into the? Dropniore 's kitchen, and the little teletype 
was soon spitting oYit columns of figures. 
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"Here, in BASIC, is the program I wrote: 



5 PRINT ,f X", "Y", "SLOPE" 
10 C » 11.7E-5 
20 S * .5 

30 FOR X * 0 to 30 STEPfS 
40 M « 4E-7 * X*2 - C > 
50 PRINT X, 1E5*Y, 1E5*M * 
60 Y 38 Y + M*S 
- • 70 NEXT X 

"To get a first guess for C that would riot be too 
wild, I assumed that the slope would be zero when x was 
about* 20 feet. Solving the DE for C, . • • , 

0 - 4 x 10" 7 (20) 2 + C 

Q « -16 x 10~ 5 

"I entered, this guess for C on line 10, ran the* 
program' and observed the final value of y. I then 
changed C in a way that kept reducing the final deflection^ 

"The vaJLue of C that I finally used (11.7 x 10" 5 ) 
gave a deflection at the end of 0.05 x 1Q" 5 ft, $s compared 
* to the maximiftu deflection of 136 x 10" 5 elsewhere on the 
beam^ * * 

- "Here is a copy of the final RUN. 
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X 




SLOPE 


0 




-11.7 


.5 


/ - 5.85 


. -11.69 


1 


-11.695 


-11.66 


1.5 


-17.525 


-11.61 


2 


-23.33 ' 


-11.54 


.2.5 


-29.1 


-11.45 


3 


* -34.825 . 
-*0.495 


-11.34 


3.5. 


-11.21 


4 


-46.1 ' 


-11.06 


4.- 5 


-51.63 


. -10.89 


5 


-57.075 


-10.7 


5.5 


-62.425 


-10.49 


6 


-67-67 


-10.26 


6.5 


-72.8 


-10.01 


7 


-77.805 


- 9.74 


7.5 


-82.675 


- 9.45 


a 


-87.4 


- a. 14 
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Maximum 
Deflection 
Here 



X 


Tsx 10 5 


SLOPE 


8.5 
9 ^ 


-91.97 


- 8.81 

- 8. ,46 \ 


-96.375 


9.5 


-100.605 


- 8\09 


10 


-104.65 


- 7.7 


10.5 


-108.5 


- 7.29 


11 


-112.145 


- 6.86 


11.5 


-115.575 


- 6.41 


12 


-118.78 


- 5.94 


12.5 


-121.75 


-> 5.45 


13 


-124.475 


- 4.94 


13.5 


-126.945 


- 4.41 


14 


-129.15 


- 3.86 


14.5 


-131.08 


* - 3.29 


15 


-132.725 


- 2.7 . 


15.5 


-134.075 


- 2.09 


16 


-135.12 


- 1.46 


16.5 


-135.85* 


- .81 


17 


-136.255 


.14 


17.5 


(-136.325) 


.55 


18 


-136.05 


1.26 , 


18.5 


-135.42 


1.99 


19 


-134.425 


2.74 


19.5 


-133.055 


3.51 


20 


-131.3 


4.3 


20.5 


-129.15 


5.11 


21 


-126.595 


5.94 


21.5 


-123.625 


6.79 


22 


-120.23 


7.66 


22-5 


-116.4 — 


8.55 


23 


-112.125 


9.46 


23.5 


-107.395 


10.39 


24 


-102.2 


11.34 


.24.5 


- 96.53 


12.31 


25 


- 90.375 


13.3 


25.5 


- 83.725, 


14.31 


26 


- 76.57 




26.5 


- 68.9 


~ 16.39 


27 


- 60.705 


17.46 


27.5 


- 51.975 


18.^5 


28 


- 42.7 


19.66 


28.5 


- 32.87 


20.79 


29 * 


- 22.475 


21.94 *' 


29.5 


- 11.505 


23.11- 


30 • 


5.00000000E-02 


24.3 



"Notice that the maximum deflecfion occurs just 
over 17 feet from the end of the beam, ^between the safe 
and tne midspan, as Angelica had predicted." ^ . 
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Arcletrose slowly to fiis feet. "And that, dear 
class, is the end of the story. Denis* and Angelica 
bath regained t.heir jobs, and were more highly regarded 
than ever before. Convinced of the power of mathematics-- 
in the Hands of an expert- -they are both studying calculus 
in their spare time. I visit them often, giving them 
small tips to facilitate their studies, and the affectidn 
and respect they shower ujfon *me is almost embarrassing. 

"Have a -good weekend, and remember that we have a 
test on differential equations on Monday." He strutted 

• ' . A- - • 

* We sat there silent for a long time; heads hung down 
and shaking slowly from sidf to side. 

Finally Herb spoke, "Mat, didja call?" 

If • "Yeah* They never heard of Arclet!" • ' 

"Chea£ <trick, ""Polly said. 

"We've lTeen had. Taken, conned,' duped, used." 

« ^ 

"Should we go through with the plan?" 

••Let's vote," « 

Three thumbs, pointed downward. *That night, I phoned 

♦ 

my brother irt Seattle. 
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Chapter 13 
1 THAT EXAM AGAIN 

Today, Arclet seemed a different person. His hair 
was brushed, he wore a well-pressed suit with a clean 
shirt, and his shoes were shined. He strode into the 
classroom , carrying a smart little attache case with metal 
trim instead of his battered, old briefcase. 

"Who cares \bout differential equations?" he asked. 

""I'll tell you »who caresf. Boeing Aircraft Company on the 

West Coast cares, that's who." He "waved a telegram at us, 

♦» 

"Listen." 




HAVE URGENT PROBLEM INVOLVING DIFF EQN. YOUR 
NAME RER TO US. HOPE 'YOU CAN COME ASAP. YOUR 
USUAL FEE PLUS TRAVEL EXPENSES PAID. 

"I've arranged for a*subsfitute teacher, who will 
be here starting tomorrow. Seq you in a week. v " 

Halfway out the door, he stopped short/ "I nearly 
forgot.* Finish this *test today and leave it with my 
'substitute." 

It was exactly the same exam that Arclet had used « 
as a pre-test for this topic. In fact, the same test 
that had prompted his wkrole ^ridiculous story. 

We Kad come full circle. 
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. THAT EXAM 

1, Given the differential equation 

r 

'. - •" #u?y' 

^ 3xy 2 
determine whether the equation . 
, y m x V l - In x * 

is a solution. 

2. *0nly one of the curves graphed below is a possible 

solution to the differential equation? 



* ~ y 



Check enough points to determine which curve represents 
a solution* 



:ifct 



.ttc 



tttt«f: 
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Given the differential equation 

verify whether the following table of polrfts represents 
a solution by checking at least 3 points. 



0.5 






1.0 




-2.000 


1.5 




■»1.184 


2.0 




-0.793* 


2.5 




-0.568 


3.0 




.-O.^ 


3.5 




-0.323 


k.O' 


•* 


-0.250 


k.5 




-O.^* 


5.0 




-0.153 



2 ' 



Plot a tangent ffeld for the di fferential ^guation 



X 2 -*2i/ 
o dx m x 

where 1 » x 3 5 and 0 £ z/ £.5. 

Sketch In 2 possible* solutions. 

Given the differential equation 



x 2 - y 



with the boundary conditions 
y ■ 1 when x « 1 , 



find a solution graphically, taking Interval^, del^j x 
of 1 and find the value of y when x * 7. 



• 17 



^^••r^^H^^- - ' "V"' '• • 



*" L 
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6» Given the differentia} equation 




* 



w %v and the boundary concfitlons 

J i s. *• - * v * • * " ■ ** * 

, V • . y m a y* 1 ** «. * o » * - v , - *■ 

I*' *f * numerical solution^ takhfc intervals, delta * of 1, \ 

.and find "the' value of t/ wheQ'X*« 10 * 

7* Write a dlWrentfal equation to descYUbe the following" . ^ ' « • 

\ 4' situations. We define <your symbols/ V lV\ J 



_APP£NDU£j| V 

A Answers to Quiz #T . . . . „ 

B Answers *to Quiz #2 
C Answer to Quiz #3' 
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a) A body *all& In a medium offering resistance pVopor^ 
tional to the speed -at any -.instant. * } ^ 9 

b) A part (eft' moves in a "horizontal line a^Ted upon by 
an attractive force which varies Inversely as the 

• . cube of the distance: from a fixed point. * • 
' *c) The ratejpf flcfr frog,* i tank of uniform cross- 
section is proportional to the square root .of the - _ 
> \ UguU depth. . \ , fc . ' * ' 
. *d) ftadium decomposes at a ratr proportional to. the 
*\ present amount 1 ' >r^" ^ / 

Poke a hole in Arclet J s fish story by cajculattng^ho^ ^ 
ajahy f|sh t Mike's freezer wou^W have to hold. 1 



D Answer to Quiz #4 v . . /71 
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E Answers to That Exam 
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H V ?1l " ? 



*- 1. Let: «■ temperature at any # time t ? c 

* temperature of surroundings 

♦ m 
then, s • 

^ ^ which is known as Newton's law of cooling. 

2. Let i? be the resisting force of* the medfum on the object* 
The statement of the problem gives. 

(1) . ' R m -k/Z. \ m * 

Assuming fl^tb t>e the only force acting on the object, Newton's 
second law gives ~ ' ^ 1 

<2H- ma = f 

where, m is the object 's^mass and a is the object's acceleration. 
„ f Substituting (2) Into^(l) gives 
x . te) " " ma » r fa>. ^ ' 

Acceleration .Is defined as the rate of change of velocity, or 



which when substituted into (3)' Qives » • 



a differentiated equation describing, the sysjdm. 

If additional ^orces^ acted 'on the" body/ equation (2j would have 
to be modified^ account for them. One addlt Jonal'force, F. 
it m driving the object' through the. medium would, lead to the follow- 
* I* In* differential eqoa^tlon:,^ . . 

M . Jv -See'te^t page W. . l . :' ' " « • . * 

J '• 




, Since the agreement is excellent, we may-assume that the equation 
y *4 (lr. x + 1) 



is a solution of the differential equation 
&L m x + 2y m 

in the domain shown by the graph. , 



2. . Calculation of ^ for ^ - ' % «2-4' 

N ■ . * ax ax x x 























1 to 


ft 


J t 


0 


0 


2 


\ 


» J ' 


i 


. 1 


4 




• 2 . 


V : 




• 3 


3 


-1 


2 


0 


* 0 ' 


2,0 




1 


0.5 


.1.5 




2 


. no 


1.0 




3 


. 1.5 


0.5 




d 


0 


2.0 




l . 


• .25 


K75 




2 


■ .50. 


1.50 




3 


.75 


.1.25 



Tangent i Field, for 



.etc 



Ith two possible solutions. 



■r*. ~- ^ 



* J /- 



'7 i 






» 










<N ^ 






f 


'— 1 






> 




L i 












f- 


• • 

-i— **- - - ' 


* i 






f '■ 




i < / 














Answer to Quiz #4 



1«* Given: 



x - 2x*y 



tn 9 



% . and Euler's Method: y n - i/ n - x + nr n ^ Ax, 

-Y- ' ». ' 

^ 4 we let ir 0 ■ 3i when - l»and Ax - 1 for this calculation. 
^' It Is easiest to make a table for calculation as follows: 



'•1 

2 
k 



1 «rf 

*Z 2.%00 

fc 3 2:26a 

*• 2.200 

5 2.167 



m mAx 



S 



2'— 19.20 

2.260 
3 - 40.68 



-0.600 : 



2.200 2.200 . „,„ 

4 - 7o.4o "SOo " 0 • 033, 



t/ + mtx 



3*00 - 0.600 
2.400 - 0,140 

2.260 - 0 r 060 

* - - » 
2.200 - 0^3- 



*; * ^ The value of ^obtalnejl^when x -,5 will depend upon the step size 
^ used^as shown li\ the table — — * 



* step size 


«/<5> 








2.167 • 


0.1 


2.521 




.2.5*9., 


0.001 ' 


2.551 



p 



Answers to That Exam 

1. Differentiating the proposed solution, 
„ ^ a 2*- 3 Zn a» 
3(1 + In x)*l\ 
Substituting y and its derivative Into the PE^we get, 

'2 - 3 In x _ 3Cx/l - in x) 173 } 3 - x 3 - 
3(1 + Zn x) 2/3 ^ 3x{x(l - Znx)* 73 } 2 



/ • 3x 3 - 3x 3 Zn x ^ x 3 

, - a 2 - 3 Zn .x 

- • 3(1- In x) 2/3 " 

t . which verifies the proposed solution 4 . *' 

2. Curve 3 is the only curve shown that can be a solution 9 . ' 

\ To determine If a curve Is a possible solution, plcjc a convenient 

point on the*cuQ£e,,# measure the slope of. the tangent line at that" 

* * . / 

point aYid compare the measured slope wfth the derivative calcu 

r lated by substituting the coordinates of the point Into the 

equation for the derivative, 

dy • 1 



dx x - y. • • 

Curve It At the point (-1,1) the measured slope = 0, b 




» The measure^slope of 0 rs-not undefined, therefore e 
the solution ^nnot be^curve 1. 
Curve, 2: At (9.50, 1,96). the measured slope = 0 (at least < |0*02|), 



but • du 1 1 

* • 1- 9.50 -<i. 9 6 '7iy°- 13 ' . * t 

and 0,13 Is not < |0,02|, therefore the.S9lution Is no.t 

likely to be curve 2. ' . * . 

CheckJng another point on curve 2 we have at (0,0) the- 0 

measured slope » t.1, but # < * ^ * 

du % \ J J < <> • ' 

fa m q - q "'q* m undefined, and/ 1.1 Is not 

undefined, therefore the^olutton cannot- be etirve 2., * 

■ \ i 185- * • 



Curve**: At (2,0) the measured slope = 0.5_<and - 

Therefore the solution .could be curved. 
Curve 4: At (2,0)* the measured slope = 0 a,nd ' ' 



* Therefore the solution cannot be curve W 4V 

Since the solution Is supposed to be one of the curves, and it 
cannot be curve-1, 2, or 4, bur could be'curve 3, We might assume 
that It must be curve 3. It would be wise, however,* to check 

another pointer two on curve 3. ' % 

1 . 6 

of 

Curve y. At (4, 0.8) the measured slope = 0.32 and 
A m 1 ' — K V 

• • . • £ "Pol" IT" °- 312 5. . - 

which is in close agreement with th^measured slope'. 
At (7.5, 1.6) the measured-slope = 0.175, and 



• r ^"7.5 \ 1?6 ' 53" °-tf95\ . 
# agafrv In close agreement with the measured slope. > 

• « N 

• * '.. . * 
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3. The gfven points' are^ plotted^below, and the slopes at x = r, 2, 
and 3 are^ measured. These compare well with the slopes computed 
from the DE, as shown in the table. 








y 


Measured 
Slope 


3 - xu 

>2x z ' 




I 


-2.000 

* 


. 2.40 


2.50 


f 


2 


-0.793 


• 0.55 


6.57 . 




3 


-0,423 * 


. 0.22 


0,24 
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k m Calculation of slopes at different* points for 

» o -» •* 

• t 

ax x 







x* '- 1y 

t i 




1 


0 


1 - 0' 


1 




2 


' 1 - * . 






4 


"8. 


-7 


<• 


P 

2 

t 


9 r'0< 
•9 t 4. 
, 9 - 8 


3 

5/3 * 
"3 . 



This'table is continued and the Jesuits plotted below. 



V 




A numerical solution wi thyintervals of 1 unit yields the 
following table of point pairs. 



X 


y 


A 

u 


1 .00 


1 
1 


1 . 00 


2 


m cn 


} 




k - 


2.66 




^3.26 


' e 


3.86 


i 


M8 * 


8 


5.08 ' 


9- 


5.69 » 


10 


6.31 



\ 



(For method see answers to Quiz §k, page 75.) 

Let IT-be the weight of.the body, and R be the resisting force. 
Then by Newton's second law. of motion 
£ F - ma 1 ' 
■ f 

y + A 



but 



/? ■ -fa; where v ■ instataneous velocity, so 



W - kv - m 



do 
IB 



since the mass m ■ W/$ /> 

cfy . kg 

Let* be the distance from the fixed point, and F the attractive 
force. Then,, 

k 



F - 



ma 



v 1 2 
where Is the particle mass and a is the acceleration, , 



then\ 



k d*x 
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» * * * > 

7c. Let R * the' rate, of ftow from thetank, and y ? the l iquid depth v 
Then, from the statement^ of the peob'lem, 

P - Ci /y . 

But, the rate of flow must also be proportional to 1 the'rate at 
which the water level is changiog. , Thus, 

Setting the two expressions for ratev equal to each other we have~ 

1 Ct 

Dividing by -C 2 , and letting = k we have finally # 

8.. The number of fish collected in Mike's freezer would equal the 
difference between the number of fish that came into the pond 

and the number left in the pond when Cy Seepage counted them, 

i 

or A (fish in] m [fish into] m [fish left] . 

(freezer)*, [ pond j " [ in pond J • 
i 

^Arxlet, used the tangent field method and calculated the fish left 

in the pond to'l>e about 3500 fish (see figure 9c)., 

* /' 
The fish into the pond would.be the number that viere carried in \ ' 

by the water plus the number put in by the Fish and Game 

Department, that is j * " 

(12:960 min)(0.6v fish/min) + 1000 fish « 7776 fish + 1000 fish 

- 8776 fish x 
^ 4 * * 8800 fish 

Then the total number of fish caught in the basket would be 
8800 fish - 3500 fish = 530G-fish 

Some freezer! •> 

^ - Some dinner I s 

\ 
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* Prerequisite Skills : r t f , (¥ ^ * 

1. Ability to integrate * V dt. ' , \ 

JO f(t) 

2. Ability to use first and second derivatives as aids in graphing 
functions. 

3. Knowledge. of lim'e" C (k positive). / . ♦ T 
Output Skills ; 

1. Know equations governing radioactive chains ("parents and « 
daughters").. • 

2. Know the meaning o£ transient equilibrium and secular equilibrium. 

3. Knew the approximations relevant to transient and secular equi- 
librium* and knoj/ the circumstances under wh'ich they are applicable 

Other Related Units : o> 
/ Kinetics of Single Reactant Reactions (U232) 
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* RADIOACTIVE CHAINS: PARENTS AND DAUGHTERS 



1 • INTRODUCTION 
1.1 Radioactive Decay 

, ; Radioactive decay is a first order reaction. This 
means that if a radioactive substance is not being re- 
plenished in any way, then its amount (number of atoms) 
N(t) decreases at a rate proportional to that amount:*. 
(1) N'(t) = -AN(t), " . ' 

.Where A is a positive constant known as the disintegration 
constaht 'or decay constant. 

The elementary consequences of Equation (1) are , 
discussed in many elementary, calculus textbo.oks. In our 
Unit 232 (Kinetics of Single Reactant Reactions f we* 
' discuss first order reactions in greater detail. In that 1 
unit.you can 'discover how experimenters determined empir> 
ically that radioactive decay is a first order 'process , 
and what this suggests about the^mechanism of radioactivity. 
1«2 Chains ^ , 

^When a radioactive substance A decays into a substance 
B, A'and B are referred to as the parent and the daughter. 
It may happen that B^ itself is radioactive and is the 
parent. of a new daughter C; and so on. In fact, this is 
a very common situation.* There are three , chains like this/ 
beginning respectively Wf ith U 238 , U 235 , and Th 23? , whose 
lengths ar^e 19^^17, and li. They do not overlap, 'and 
together accouift 5or all naturally occurring radioactive 
substances beyontl Thallium (atomic number £1) Jn the- . . 
periodic table. Esich of these chains entfs with a stable 
(non- radioactive) form of lead. ' * 
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SETTING UP .THE EQUATIONS 



2.1 Notation and Assumptions \ 

We shall consider the relationship' between one 
parent A and her radioactive daughterly. We shall write 
Nj(t) and N 2 (t) for their amounts^at time t, and \ l and 
* 2 for their decay constants! Figure Arnay help you 



N t (t) 



N 2 (t) 







B 








► 


► 





Figure 1. Schematic representation of a radioactive chain. 

remember this notation. Since yt he rate of decay' of B into 
C depends only on the amount of'B present, and not on the 
amount of C, we do not care whe'ttier C is Stable or .radio-*** 
-active. * 

, Now imagine that at'the instant' t = 0 we have a 
freshly prepared .amount N Q of A, and none of B. Th-at is, 
NLj(O) = N Q and N 2 (0) = 0. Imagine that the chain of % 
reactions hi Figure* 1 then proceeds without external 
interference. ^ - * 

2.2 The Equations '. ^ . 

* — & 1 1 

Since A^is not being replenished, Equation, tl) applies 

directly, and we have* 

Nj(t^ - -A^ft). . ^ 

If B were not being replenished, Equation (1) would 
apply again, and N 2 (t) would be changing at the i^'te < 
-A 2 N 2 (t). But B iq being replenished. * Each, atom of A 
which decays becomes an atom of B, and this' is happening 
at the rate A^ft). So altogether we have 

N*(t) = A^Ct) - A 2 N 2 (t).- 
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(3) 



We, are confronted with the following system of 
equations : 

(20 • Ni(t) = -X.N.(t) 

• v. • \ ■ * ' 

^.(t) = XjN^t) - X 2 N 2 (t) . 
NjCO) =~H Q . ^ 

N 2 (0) '= 0. - > 1 

, 3. SOLVING THE EQUATIONS 



3,1 .Solving for N ^t) 

It is fairly straightforward to solve Equation (2) 
for Nj(t)* This was done in Uni v t 232 and is probably done 
in your calculus textbook. We just diyide through by ~H 1 (t) 
and then- integrate from 0 to t: 

P.Ni(t) 



dt = - 



ft 



Xj dt. 



" Jo NlW 
This leads to "the equation: ■ 

t ln.CN^t)) - In (NjCO)) 



or 



since N 1 (0) = *The usual absolute value signs are not 

antu 



needed, because tfie .quantifies involved are positive 
.Finally, ' 4 s . . * 

(4) 



Exercise 1. 



\ Find a relationship between X.j and -the half life of A (the half 
life is the time t* at which N^t*) » JiN Q ) . 
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5.2 Solving for N 2 (t) 

Finding N 2 (t) is a bit more tricky. Applying 
Equation (4) to Equation (3) we get 

(5) ' N'(t) = XjNpe^l* - X 2 N 2 (t). 

Equation (5) probably looks quite different from any 
you have seen before. Let f s try to make a shrewd guess 
what kind of solution it has. It^^s^that the derivative 
of N 2 (t) is the sum of two terms, f^N 0 e" X l t and -X 2 N 2 (t). 
With luck, this nfight' remind us of the product rule: 

(6) if N 2 (t) *- u(t) --v(t) 

then N 2 (t) = u(t) • v ! (t) + v(t) • u.»(t). 

Can we' pick u(t) ajid v(t) so the terms in Equation ($) { 
match up with the terms in Equation (5)? In;other wo f rds, 
can we pick u(t) and v(t) so that 

(7) * u(t) • v'(t) = X-Nne^l* ' - 
and 

(8) ' j ,,v(t) • u'(t) = -X 2 N 2 (t)? ' \ 

Since N 2 (t) = u(t), "%(t), Equation (8) caft 6e rewritten " 
v(t) • u'(t) = -X 7 u(t)v'(t) 

and we are, in busines^! The v(t) factors cancel out, leaving 

us with 

* • / 

* u'(t) = -X 2 u(t) ' + * > 

which looks very much like Equation (2) and can be solved 



in the same way. First, 
Theri, writing R = u(0), ^ 



U(t) 
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Putting this in,to Equation (7) gives 
Re" X 2V(t) = X 1 N 0 e" X l t 

. ' v / (t ) = tp. ,CV X 1>*. • • .... 

If X 1 ■* X 2 we 'feel confident you can complete this 
solution yourself (see Exercise 2). ^* 

If Xj f X 2 , then X^ - X^ f 0 and we can write 

» • # • 

^ where K is the constant of integration. Then 

N 9 (t) = u(t)' • v(t) = ^l^L e' X l t + KRe' X 2 t . 
z A 2' A i 

Using the fact that: N 2 (0) = 0,'we get • » 

. x i N o • m 

a 2 -a 1 

; \ ^ . 

KR = - T±jL t 
A 2 A l 

X 1 N 0 , -x n t -. / -x~t, 



Exercise 2 . . # 4 

T~ 

*ind N 2 (t) if X L « X 2 . 

• # 
Exercise 3. 

Assuming C is stable, .find the ^ime at vhfch the total radio- 
activity {i.e., the total number of disintegrations of A-atoms and 
B-atoms per unit time) is greltest. 

- ■ 1 ■ 1 j 2 \ 

3.3 Graphs of N 1 (t) and N ^ -ft) 9 . • * 

It is easy to confirm from, Equation' (4) that N|(t) < 0 
^ for, all t, that N ! , f (t}->0 fens' all t, and that lim N,(t) = 0. 

^ With a little more work (see EocerciSe 4) it cap be con- 
« \ firmed that, lim N*(t) s 0 and that 



a' ♦ 



t-V<» 
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(10) 

and 

(11) / 
where 





> 


0 


if t 


< 




No(t) ■ 




0 


if t 




u 




< 


0 


if t 


> 


u 




r < 


0 


if t 


« 

< 


2t 0 


N 2 (t)i 




0 


if t 


= 


2t 0 




> 


0 


if t 


> 


2t Q> 


lnX x 




lnX 2 







V 

Therefore, the .graphs of N^(t) and N 2 (t) have the shape's 
shown in figures 2 and 3 .• # 




0 • ' t 

Figure 2. Typical graph of Nj(t) (amount of A as a function of time). 



MP) 




,Fi'gure 3. Typical graph of N 2 (t) (amount of 



* 0 V < 

: B as a function of 
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time) . 



Exercise A % , 

a. Show that lim N 2 (t) » 0. 

4 . b. Confirm Equation (10), 
N c. Confirm Equation (11). 

Exercise 5 ' . / 

Find 4 the time at which the greatest amount of B will be present, 
Exercise 6 



*or the chain B* 210 + Po 210 + Pb 206 , \ { « 1.37 x nT 1 day" 1 and 
X 2 - 5.1 x 10" 3 day"* 1 ^ 



aS. Use- Exercise 5 to. determine when the amount of'Po 210 will be 
* • greatest. . k 

b. If filially there are 10~ 8 grams ef Si 210 , how many*grams 

of Po wil1 t ^ere be when it is at its maximum amount?. « 

i . - ^ » — 

4". EQUILIBRIUM 

4.1 What is Equil ibrium A * > * 

^ In a continuing proqess such as the one we are d£s- 
cursing, it is natural to ask about "equilibrium" of the 
process. Webster's Seventh New Collegiate Dictionary (1965) 
defines equilibrium as "A static or dynamic state of 
balance between opposing' forces or action's. " To a scientist, 
"sfate of balance" means that certain measurable quantities 
remain constant. But in practice scientists frequently use 
the word "equilibrium" when the measurements under consider- 
ation are nearly constant rather than actually constant. * 
There is good reason for this. The process 'are discussing 
illustrates that reason nicely. We have been devoting our 
attention to Nj(0 and N 2 (t). For either of thes^e actually ■ 
to be constant over any time interval, its derivative would 
have to be zero throughout that interval. But N^(t) < 0 for 
> all t, and N^Ct) = 0 for only one t. So,, strictly speaking, 
it is impossible for either N^t) or N 2 (t) to be "in equili- 
brium." - > ■ * 

* * • < . 

But N^t) and N 2 (t)^ involve negative exponential func- 

tions. In fact, functions involving negative exponentials 
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occur fairly commonly in the description of physical and 
chemical processes. The most basic negative exponential 
function is & Z $ and if you know anything at all about 'it " 
you know that it approaches zero very fast, so that, al- 
though it is never^constant, it is before long practically 
equal to zero and therefore practically constant, 'This' 
characteristic ("never constant, but practically constant") 
carries*through to many of the more complicated functions 
involving negative exponentials. The quantities "they de- 
scribe never actually reach their limiting values, but 
usually came (and remain) extremely close to them within 
a reasonable iength of time- -perhaps even so close .that the 
difference is not measurable v Scientists often apply the 
words "equilibrium" or "steady state" to this situation. 

There is one more "thing we should say about the word 
"equilibrium* 4 before we move on to discussing specific 
cases of it. The functions N^t) and tt^(t) are examples 
of functions which come and remain extremely close 'to a 
cpnstant value (zero). - But it would be wrong to say that 
they are therefore in equilibrium, even allowing for the 
stretching of the definition which we j-ast discussedj The 
reason li*s in the other part of the definition: "oppoftftg 
forces or actions." There are no* "opposing- forces or 
actions." Rather than "having two things happen which cancel 
each other out," we have nothing happening at aM (in the 
limit). In pjain English, virtually all of A will have 
decayed' into B and then into C, so that there^will be vir- 
tually none of A or ^ left." It takes.no fancy , mathematics 
to see t*his.. If C is radioactive it will eventually decay, 
and so on, so that the limiting situation is ,that only the 
stable substance at the end of the chain remain. 

There are, however, two situations involving radio- • 
active chains to which the word "equilibrium" is usually 
applied. The first of these is known a£ transient equili-^ 
brium. Another^ known as secular equilibrium , can be 
regarded as a special case of -the first. 
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4.2 S<5me Comments on Approximation 

. Before we get into the mathematics of transient and 
* secular equilibrium, it will be wise to take a moment to 

discuss just what we mean by a "good approximation . 11 When 
scientists "say ^two numbers r and s are approximately equal, 
they almost always mean that the difference between r and 
s is small compared to either 'of .the numbers. For example 
they mighfsay 1002 * 1000 (depending on the context), but 
would almost never say 2=1. ' 

Saying r - s is small compared to (for example) s» 
. means = I - i is small, or I is near 1. In the numer- 
ical examples we just gave, ^£ = 1 .002 , which is very 
near 1, but y ■ 2, which is much further from I. 

This interpretation of approximation can be.^pplied to 
functions too. Let f s look specifically at negative expo- 
nential functions. If P and'Q are* any non-zero constants, 
» ■ and if 'a and b are constant's such that 0<a<b, then * 

pi) ' Pe ' at + f' bt , i + 2 e (-b)t ■ \ 

pe -ax p 

as t ■* «, since a - b is -negative. Therefore, for t large 
enough, Pe" at + Qe~ bt cah be approximated by Pe" at . We 
shall use this fact in the next section. • J 

4.3 Transient Equilibrium 

It has often been observed that in many chains in- 
volving a parent A -and a radioactive daughter B, after a 
^ while both parent and daughter appear to be decaying at the 

same rate, in the sense that in any given time dnterval (say 
from tj to t 2 , with t 1 large enough) parent £nd daughter t 

each lose the same fraction of their initial amount ~»/v , 

N 1 (t 2 )/N 1 (t 1 ) * N 2 (t 2 )/N 2 (t 1 ). This phenomenon is called 
transient equilibrium. Let us try to explain it mathe- * 
•matically. 

We can rewrite the equation of the preceding paragraph 
N 2 (t 2 )/N 1 Ct 2 ) = N 2 Ct 1 )/N 1 (t 1 ) . In -other words/ the observed 
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result is that N^t^/N^t) i<> a Constant . (This is why 
the word "equilibrium" is used in describing this phenom- 
enon.) Why* should this be so? Let ! s investigate this 
quotient, starting with the formulas for Nj(t) and N 2 (t) 
given in Equations (4) and (9) respectively. . 

We know from Equation (9) that 



A 2~ A 1 



If A i < A' 2 , Section 4.2 tells us that for t large enough 



(13) N 2 (tOL e~ X l Z . ' 

We also know, from Equation (4), that 

6 

N x (t) = NqB^I* 

(which is exact). Dividing Equation (13) by Equation (4), 
we get 

N 2 (t), \ x 

(14) NTTtJ * 



On the other hand, if A 2 < Aj then this does not go 
through as neatly. Equation (13) has to be replaced by 



A 1 N 0 -A 0 t 



N 2^ * T^T- e " A2 
L 12 ' 

and then Equation (14) becomes 

N 2 (t) m A i e (A r A 2 )ty 

This is a positive exponential. It does not have a finite 
limit., „ * 

So. th e mathematics tells us that transient equilibrium 
should be observed when * 1 <A 2 , but not otherwise. Sure 
enough, this is exactly what happens. 

* • • 

Another way of , looking at transient equilibrium is to 

compare the approximation given in Equation (13) directly 
with Equation (4). The exponent -A^ is the same in both 
cases. So for t -large enough B behaves as if it had the 
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same decay ^constant (and therefore the same half life*) as \ 
A. Sincp/x'j < X 2 this apparent half "life Is longer than B ! s 
natural half life, an observation which should appeal to * 
'your, common sense even w.ith no mathematics a£ all. After 
all, two things are TTappe'nin^' to* *B It 4s decaying at its y 
natural rate, and it is being replenished at a certain rate. 
Therefore you would expect its actual rate of disappearance 
to be somewhat slower than if 'it "were not being replenished . 
This c,ommdn 'sense observation may .help you remember which , 
way^tfre inequality \^ < A 2 goes for transient equilibrium. 



' 4.4 Secular Equilibrium 

, We have said that transient equilibrium occurs when \^ 
is Smaller than A 2 . Now let us suppose that 'Xj "is very 
small, . and very much smaller than * 2 - (Scientists, write 
X 1 <<X* to mean * 1 is very much smaller 'than * 2 .") This is 
. actually a very common occurrence. For- exawple , in the 
chain « ' 

,Ra 226 has, a half life of about .t*620 years*. "The decay con^- 
.tant far the first step .is =■ 4.^8 x io" yr * = ? , 
1.17 x I0" 6 day -1 . In contrast , Rn 22 has. a half lif> of 3.85 
days, so that \^ = 0.181 day" 1 = K81 x 10" 1 day" 1 ^ 

\ We know that whenever 0 < X, < e'^1* decrease.? to 

zero more slowly than e 2 1 *. If X<,<<X 9 , v the difference 

* ■ - \ t 

in these rates is so great that long after e 2 has become 
< tiny enough to neglect* in Equation (9), we oan still say 

e * 1. (When t = 200.^ays in the example given above, 
, eC X l t * 0.99977,, and e" X '2* * 2 x 10" 16 .) Then we Would have 

and * 



N 0 



(16) NjCt) « N Q , 

both approximately constant. 

— c 

*For the definition of half "life and its relation to the decay constant, 
see Exercise 1. For more about half life^see UKAP Unit 232 (Kinetics 
^of Single Reactaqt Reactions). 
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Not only are N, (t) and N 7 (t) decrea»ng*at the same 
rate-j but this rate is so slow that they are in fact vir- 
tually constant. This situation is known as secular equi- 
librium.. Again, as remarked inSection 4.1, we are 
stretching the term *a bit, since the quantities involved 
are not really constant. Here there is yet another abuse . 
of terminology, in that the "virtually* constant" values of 
Nj(t) and N 2 (t) are not their limiting values. Eventually 
(although perhaps none of us will live long enough to see 
-it) these amounts *wiLl begin to decay noticeably, and 

ultimately they will approach zero. 

m i — 

Rememb§r*we also remarked in Section 4.1 that for ail 

their abijse of the term, scientists dd*agree that it is 

wrong to apply the wor,d "equilibrium" to a situation in 

which "nothing is happening" (there are no opposing forces 

or reactions). t This is not a problem here. Plenty is 

happening. New B nuclei are being formed, and 'old ones are 

decaying. The total number of B nuclei remains the same, 

but they are not at all the same nuclei. t (The total number 

of people in New'York City is about the same as forty years 

ago, but they are certainly not the very same people.) 

4.5 More Comments on Approximation 

One thing* about Section 4.4 may puzzle you. Adding 
the approximations given in Equations (16) and (15) we 
get \ % « 

' N l( t)< + N 2 (t) - N Q ♦ ^ N Q > N Q . 

But the total number^ of atoms, including those of C and 

possibly later substances, must always equal N Q . How can 

this be? v # 

* * * 

What Has happened is that Nj(t) has decreased by a 

certain amount while N 2 (t) has increased by a lesser amount. 

But there was a lot* of A to begin with, so the decrease is 

small compared to the original amount, and NjUJI/Nq *^l. 

On the other hand, there was none of B to begin with, and 

even at secular equilibrium the^e^^ very little. Compared 



1 . ; • . 

to this. amount, the increase* is significant. It's as if 
Exxon Oil- (A) were to pay you (B) $10/000** and you were to 
use $10 of that money to bribe your math , teacher (C).. Ex*on 
will still have essentially the^ame amount of money *as be- 
fore,, and you will be much "better of ^financially ,^even 
thought the total ofiExxon's money anci your money, when, 
calculated precisely, ^tfill be less than before-** The 
numerical calculations in Exercise T may help you see what 
is going on. t ' , 

' Exercise 7 >' * 

22 2* 

/ a v Use Exercise 5 to show that the amount of Rn ' in the chain - 

of Section -4. 4 is greatest at about t = 66 days- . , 

# b. (Requires a calculator.) For the values of w t given below, 

X • ?9fi 299 

Is compute the precise amounts of Ra and Rn as given by 
Equations (4) and (9)i and v also the sum of these amounts. 
Then compute the approximations^ given by Equations (15) and , 
(.14), and also their 'Srum. Tabulate and compare these results. 
» Take tf Q « 10? atoms, and make all computations to the nearest 

" integer. - 
Use t (in days) - 1,2, 3,4,^,^,7,8,9, 10, 20,30, 40,50, 60,66, 70, 
* 60,^0,100,^00,300,1000. 



**ou will note tha't for this analogy it does^not matter whether your 
math teacher is stable or- not. 

\ ; ' Q » 13 

5? ERIC 208 



5 .* ANSWERS' TO 'EXERCISES ' 



-X.t* 1 
e 1 



-Xjt* = ln| - -In 2 



t* 



Writing X - ^ « X 2 : * 

XN ■ . 

v(t) = -g* t + K 
N 2 (t) = u(t)-v(t) - R&" Xt (^ t + k] = XN 0 te" Xt + RKe 
Since N 2 (0) = 0, we have RK - 0 and 

i 

N 2 (t) = XN 0 te" Xt . 

Set D*(t) - 0, where D(t) - X^^t) + X^fi^t) . 
D*(t) - XjNjd) + X 2 N*(t) . . 



- -X 2 N e'V + XlXz "° 



U 2 N X ? X 2 N 0) -X.t , X 1*2 N Q -X,f 



D'(t) - 0 when 



2 2 

x;x 0 \ , _ x,x; 



fx 2 . 12 -x 9 t - 

l x i + xpx7J e 1 ; jpr^ 2 

-2 2 

MVV + X i X 2 m e <x r x 2 )t 

X 1 X 2 



20^ 



* X 2 



i 



t - 



x r x 2 



In 



"1 



T (2X 2~ x i> 

u A 2 



lnX 2 + in(2X 2 -X 1 ) - 2 In X 2 



x r x 2 



4; a. Use e^l' + o'and e" X 2 t ->-0. 
X,N, 



b. N 2 (t) - 1 ^2-(-X 1 e~ X l t + X 2 e" X 2 t J - 0 



when 



. -X.t . -X«t 
Xje 1 .» X 2 e 2 



Since N 2 (0) ^0 
and lim N 2 (t) - 



from^Jiysical considerations), 
1 inequalities for N 2 (t) follow. 



when 




(Use the. fact that lrt jyij - 2ln~.) 

\ in A - - in A« 

5. This occurs when N£(t) = 0. That t^s^.at t Q -\ jj*— j — , by 

. Equation (10). 1 2 

6. a. By Problem 5, ; 



t „ In (1,37 x 10" 1 ) - In (5.1 x 1Q~ 3 ) 
° N 1.37x 10" 1 - 5 ,1 x 10~ 3 



24.95 days. 
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IS 



t». N 2 (t Q ) 



2 A l 1 



X,N 
X 



(1.37 x 1Q' 1 ) x i(T 8 



5.1x10" 



1.3lx 10" 



-1.37xl0" 1 x24.95 • 
e / -e 



■5.1xl(T 3 x24.95 



■ 8.81 x 10" grams. 
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7. Column A gives the exact amount of Ra 

* 222 
Column B gives the^ exact amount of Rn * . 

Column C gives the sum of columns A and B. 

226 

Column D gives the approximation of the amount of Ra 

o 

Column E gives the approximation of the amount of Rn 
Column F gives the sum of columns D and 
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t 


A 


B 


C 


D 


E 


V 


F 


. 1 


999,998,830 


1070* 


QQQ QQQ QQfl 
777 , 777 , 771/ 


10 9 


6464 


1,000,006,464 


2 


999,997,660 


1963 


• QQQ QQQ £9? 
777 , 777 > J 














3 


999,996,490 


2708 


999,999,198 














4 


999,995,320 


333b 


999,998,650 














5 


- 999.994. 150 


3849 


QQQ QQ7 QQQ 
777,77/ ,777 


t 








m 


6 


999,99w)80 


4282 


QQQ QQ7 

777,77/ jZOZ 














7 


^99,99^810 


4643 


999,996,453 














8 


999,990,640 


4945 


999,995,585 t 














9 


999,989,470 


5196 


999,994,666 




• 






> 




10 


999,988,300 


5406 


999,993,706 














20 


999,976,600 


6291 


999,982,891 














30 


999,964,901 


6436 


999,971,336 














40 


999,953,201 


6459 


999,959,660^ 














50 


999,941,502 


6463 


999,947,965 














60 


999,929,802 


6464 


999,936,266 






\ 


r 






66 


999,922,783 


6464 


999,929*247 














70 


,999,918,103 


"6464 


.999,924,567 














80 


999,906,404 


, 6464 


999,912,868 














90 


999,894,706 


6463 


999,901,169 










— 




100 


999,883,007 


6463 


999,889,470 












200 


999,766,927 


6463 


999,772,490 








300 


999,643*062 


6462 


999,655,523 


I 


1 


j 


1000 


998,830,684 


6457 


99^,837,141. 
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STUDENT FQRM 1 
Request for Help 



Return to: 

EDC/UMAP 

55 Chapel St, 

Newton, MAr-02160 



Student: If yoi^ have trouble with a specific part of this unit, please fill 
out th% form, and take it to your instructor for assistance. The information 
you give* will help the author to revise the unit. 

Your Name 



Page 


£ " 


1 




0 Upper 
O^iddle / 




Section 




. OR 


Paragraph 


V OR 


0 Lowef 






* 



Unit Ntf, 



Description of Difficulty: (Please be specific) 



Model Exam 
Problem No.^ 

Problem No. 



; r 



Instructor : Please indicate your resolution of the .difficulty iir^his box. 



o 



Corrected errors in materials. List corrections here: 



o 



Gave student better explanation^, example, or procedure than in unit/ k 
Give brief outline of your addition h6re: 



\ 



r () Assisted student in acquiring general learning and problem-solving 
.skills (not ufeing examples y from this unit.) 



-/ 



o 1 
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Instructor's Signature 



ERIC 



Please use ^reverse if necessary '. 



Return to: 

STUDENT FORM 2 EDC/UMAP ' " 

Unit Questionnaire V * 55 cha P els t- 



Name_ ^ 1 Unit No. ' Date 



Newton, JIA 021 60 



Institutidn . \ Course^No. 



Check the choice for each question that comep closest to your \ personal opinion. 

1. How useful was the amount of detlil in the unit ? ■ N , ^ 

Not enough detail to understand the unit, 

Unit would have been clearer with more detail 



^Appropriate amount of detail 



Unit was occasionally too detailed,' but this was npt distracting 
Too mudh detail; I was often distracted " * \* 



.2. How helpful were the problem ^answers ? 

Sample solutions were too brief; I coulti not do th^ihf ermjftliate steps 

* Sufficient' information was given to solve the problems 

Sample solutions were too detailed; I £idn f t need- tHem 

3. Exc'ept for fulfilling the prerequisites » how much did you use other sources (for 
example, instructor, friends, pr other books) in order to understand the unitff 

A Lot Somewhat » A Little • - Not at all 



How long was this unit in comparison^ to the amount of time* you generally spend on 
a lesson (lecture and homework assignment) in a typical math or science course* 

Much ■ Somewhat About * Sqmewhat , Much 
Longer Longer ' " the Same ; S horter S horter 



5. Were any of the" following parts oi; the unit confusing or distracting ? (Chack 
* as many as apply.) 

* ' \ * 

• P rerequisites * 

Statement of skills and concepts (objectives,) * - < 

" P aragraph headings . • ^ o 

Examples ^ 

■ ♦ Special: Assistance Supplement b (if present) K - ' 

Other, please explain ; ; , 



Were any of the following parts of the unit particularly helpful? (Check as many 
as apply.) , ^ ^* 

I P rerequisites' * - ^ 
^Statement of skills and -concepts (objectives) . • m 



^Examples, 
Problems 

mm g 

Paragraph headings J 
JTable of Content's 

^Special Assistance Supplement (if present) \ , ^ 
Other, please explai n * * ; - 



Please describe anything in the unit that' you did not particularly like. 



Please describe anything that you found particularly helpful, ^(Please use the back of 
this sheet if you need more space.) j 2 ^ z — — ; 
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